VOLUME  m 
SYSTEMS  PHASE 


CHAPTER  10 

AEROELASTlCnY/ 
AIRCRAFT  STRUCTURES 


1997011T  013 


USAFTEST PILOT  SCHOOL 
EDWARDS  AFB,CA 

Mo  \/,  90 


i) 


sc 


12.1  INTRODUCTICN 


Aeroelasticity  is  the  science  v^ch  deals  with  the  iraitual  interaction  of 
aerodynamic,  elastic,  and  inertial  forces  and  the  aircraft's  structural 
responses  to  these  .forces.  If  aircraft  could  be  designed  as  perfectly  rigid 
structures,  then  the  interaction  of  these  forces  wovild  not  be  important.  The 
wei^t  penalties  and  the  resultant  loss  of  aircraft  performance  make  this 
design  option  impractical.  This  chapter  is  intended  to  introduce  the  flight 
test  crew  to  the  subject  of  aeroelastic  phenonena  in  aircraft.  Paragraph  12.3 
deals  with  material  used  in  aircraft  structures,  including  stiructxnral  design 
considerations  and  use  of  cotposibe  materials  in  aircraft  structures. 
Paragraph  12.4  introduces  aircraft  structural  response  to  loads  and  Paragraph 
12.5  presents  static  and  (^mamic  aeroelastic  phenonena  due  to  the  interaction 
of  the  forces  mentioned  above. 

This  subject  area  has  taken  on  increased  importance  to  the  flight  test 
crew  in  recent  years.  New  high  strength-to^weight  and  high 
stiffness-to-wei^t  materials  are  being  used  in  aircraft.  New  aircraft  design 
and  fabrication  techniques  allow  the  structural  weight  to  be  minimized,  thus 
increasing  the  flexibility  of  today's  aircraft.  Active  feedback  control 
surface  actuation  to  alleviate  dynamic  loads  and  flutter  on  aircraft  is  being 
iKed  to  further  decrease  aircraft  structural  weight  and  thus  gain  increased 
performance.  New  conposite  aircraft  structinres  require  the  test  crew  to 
better  understand  materials  and  structural  response  to  loads.  This  chapter  is 
intended  to  introduce  the  test  team  to  this  important  subject  area. 
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12.2  ABBREVIATICNS  aND  SYMBOLS 
A,  a  airplitTjde  or  area 

ac  aerodynamic  center 

b  semi-chord  or  torsional  damping 

c  wing  chord  or  distance  from  neutral  axis  to  outermost  fiber 

eg  center  of  gravity 

^  variation  of  lift  coefficient  with  angle  of  attack 

max 

variation  of  the  lift  coefficient  with  angle  of  attack 
a 

variation  of  the  lift  coefficient  with  aileron  deflection 

d  linear  dairping  or  diameter 

E  modulus  of  elasticity  (Young's  Modulvis) 

e  strain  or  the  distance  between  the  aerodynamic  center  and  the 

elastic  axis 

EA  elastic  axis 

F,  f  force  or  frequency 

fj^  bending  stress 

f  conpressive  stress 

o 

f  normal  stress  or  natviral  frequency 

n 

fp  proportional  limit  stress 

f  shear  stress 

s 

f^  tensile  stress 

G  Modules  of  Rigidity 

g  gravity  or  flutter  damping 

h  vertical  displacement 

I  moment  of  inertia 

j  polar  moment  of  inertia 
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^ring  stiffness 
torsional  stiffness 
bending  stiffiiess 
lift  or  length 

lift  due  to  aileron  deflection 

variation  in  lift  due  to  wing  twist 

variation  in  lift  due  to  aileron  deflection 

Mach  or  itionent 
mass 

aerodynamic  moment 

moment  about  the  aerodynamic  center 

moment  about  the  aercx^mamic  center  due  to  aileron  deflection 
load 

dynamic  pressure 

dynamic  pressure  at  divergence  speed 
dynamic  pressure  at  aileron  reversal  speed 
radius 

wing  area  or  shear  force 
axial  force 

material  thickness  or  time 
divergence  speed 
flutter  speed 
aileron  reversal  speed 

argle  of  attack  or  coefficient  of  thermal  expansion 
frequency 
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“d  frequency 

ojp  flutter  frequency 

(0^  undaitped  natural  frequency 

(0^  undamped  natural  torsional  frequency 

a  stress 

Ojjgg  engineering  stresis 

0^  true  stress 

0^  ultimate  stress 

0y  yield  stress 

e  strain 

^ENG  engineering  strain 

true  strain 

7  shear  strain 

M  Poisson's  ratio 

S  change  in  length 

^  aileron  deflection 

a 

C  damping  ratio 

<|)  phase  angle 

12.3  AIRCRAFT  STRUCTURAL  MATERIALS 


12.3.1  Introduction  to  Design 

Ihe  primary  responsibility  of  the  engineer  is  the  design,  ccaistruction, 
and  maintenance  of  structures  and  machinery,  etc.  In  his  function  as  a 
designer,  he  makes  use  of  the  principles  of  thermodynamics,  electricity,  and 
the  statics  and  dynamics  of  solids  and  fluids,  but  he  is  ultimately 
by  the  materials  at  his  disposal.  In  the  past,  design  of  a  mechanism  or  ^stem 
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has  often  beai  a  function  separate  froti  the  consideration  of  the  TnaftaT-iq-)  of 
vdiich  the  mechanism  ves  to  be  constructed.  This  process  was  adequate  when 
there  were  a  very  limited  number  of  materials  available.  Now,  it  is  estimated 
that  a  designer  must  choose  from  as  many  as  75,000  alternative  materials.  In 
addition,  the  capdDility  of  designing  specific  materials  for  an  application 
exists.  This  has  a  far  reaching  effect  on  the  process  of  design;  now  the 
designer  must  consider  from  the  outset  the  materials  and  fahT-jr!qtion 
techniques  to  be  used. 

Perhaps  even  more  significant  is  the  growing  tendency  toward  development 
and  use  of  metamorphic  materials:  those  that  change  properties  as  the  service 
environment  changes.  Some  exairples  are  metals  that  form  metal  oxide  coatings 
in  a  corrosive  environment  and  inhibit  further  corrosion,  steels  that  are 
"self-healing"  in  order  to  prevent  crack  propagation,  ard  polymers  that  are 
formed  into  final  shape  "in  place."  These  materials  demand  a  unified  design 
approach  that  considers  the  material  as  a  dynamic  part  of  the  ^stem  rather 
than  passive  and  static. 

Because  of  the  intimate  interaction  between  the  part  and  the  naterial,  it 
is  necessary  that  designers  understand  basic  material  properties  and  mattarials 
engineers  understand  the  design  process.  At  least,  they  must  be  able  to 
communicate  in  a  mutually  corrprehensible  manner  since  they  must  work  together 
throuc^iout  the  design  process.  The  purpose  of  this  section  is  to  briefly 
introduce  the  design  process  and  some  material  properties  of  concern  to  the 
"designer."  Ihe  design  process  is  discussed  first,  th^  material  properties, 
fabrication,  and  finally  an  example  of  their  interaction. 

12.3.2  The  Design  Process 

The  fonrulation  of  aircraft  performance,  size,  aid  carriage  requirements 
by  using  ccannands  starts  the  design  process.  These  requirements  are  carried 
throu^  the  acquisition  process,  modified  through  cost  aid  current  (aid 
projected)  state-of-the-art  and  off-the-shelf  technology  caiprcmises.  The  end 
product  of  the  acquisition  cycle  is  the  contractual  specification  that  i  nog 
as  specifically  as  possible  the  performance  requirements  the  system  must  meet. 
It  is  this  specification  that  dictates  the  design  and  the  design  process. 

For  a  new  developnent  of  an  aircraft  or  weapon  system,  the  definition  of 
its  mission  dictates  the  design.  While  the  iterations  of  the  design  cycle  may 
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be  sunilair  among  cargo /  fighter,  and  heavy  bcniber  aircraft,  the  tolerances, 
structural  requiranents,  and  size  limitations  differ  widely.  Ihe  typical 
flic^t  profiles  dictate  the  loads  and  mission  load  cycles  to  be  considered  hy 
the  structural  designers.  Wing  and  control  surface  sizes  and  locations  are, 
in  part,  determined  by  the  performance  requirements,  and  these  in  turn 
influence  the  structural  load  paths  to  ensxire  proper  margins  of  safety. 
Static  ax]d  ^^namic  loads  are  analyzed  according  to  the  required  flight 
characteristics  and  internal  and  ejctemal  carriage  requironents.  Sizing, 
structural  design,  material  selection,  propulsion  interface  and  support, 
avionics  locations  and  functions,  control  ^stem  development,  aerodynamic 
considerations,  radar  and  IR  signatures,  and  human  factors  engineering  are 
analyzed,  developed,  coordinated,  ccnprcmised,  modified,  and  integrated 
throu^  innumerable  iterations.  Each  nut,  bolt,  washer,  and  rivet  is  analyzed 
and  examined  to  ensure  that  critical  load  p)aths  at  worst  case  conditions 
ejdiibit  satisfactory  margins  of  safety  under  static  and  <^amic  loadings. 
^9^®  12.1  depicts  a  sinplification  of  this  process  before  aiy  aluminum  is 
formed,  conposite  is  wound,  or  fastener  installed.  The  airfrane  design 
porocess,  from  an  aerodynamic  viesipoint,  must  also  contain  sufficient  structure 
to  permit  landing  and  taxi  loads. 

Qnoe  the  basic  analyses  have  been  cotpleted  and  preliminary  design  has 
given  way  to  initial  fabrication,  mock-ups  are  used  as  a  tool  to  determine 
actual  form  and  fit  of  hardware,  avionics,  control  ^sterns,  ojQpen  and 
paressurization  systems,  ducting,  hydraulics,  and  electrical  routings.  As  this 
sizing  progresses,  structural  coiponents  undergo  load  testing  to  determine  the 
accuracy  and  fidelity  of  analyses  and  to  update  simulations.  Fabrication 
to  installation  and  integration,  yielding  structures  that  finally  begin  to 
look  and  feel  like  aircraft.  Selection  of  materials  is  covered  elsesdiere  in 
the  chapter,  but  is  a  continuing  process  in  design  in  an  attenpt  to  provide 
required  shape,  strength,  size,  and  airflcws  at  minimum  weight.  The  design 
process  continues  throu^  fli^t  testing  the  entire  ^stem,  irrproving  on 
previous  attempts  to  upgrade  performance  with  minimum  cost  and  weight. 
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FIGURE  12.1.  DESIGN  PROCESS 
12.3.3  Material  Properties 

The  value  of  an  engineering  material  is  highly  relative  depending  on  the 
specific  circumstances  under  which  it  is  to  be  ijsed.  There  are  as  many 
different  materials  suitable  for  a  given  design  as  there  are  different  designs 
to  accomplish  a  certain  design  goal.  However,  there  is  a  "best"  material  for 
a  certain  design  when  a  sufficient  nimber  of  conditions  are  set  for  the 
material  to  meet. 

There  are  a  nun±>er  of  material  properties  that  are  basic  in  almost  all 
engineering  designs.  For  sane  applications,  more  than  one  property  may  be 
needed,  such  as  hi^  strength  associated  with  good  electrical  conductivity  and 
corrosion  resistance.  As  more  conditicms  are  introduced,  the  selection 
narrows  down.  Obvioi^ly,  an  infinite  nur±)er  of  conditions  such  as  the  three 
above  can  be  set  to  finally  limit  the  choice  of  available  materials. 
Unfortunately,  the  ideal  material  for  a  particular  application  having  a 
conbination  of  optimum  values  for  various  properties  does  not  always  exist  and 
a  ccrpranise  is  necessary  at  this  stage  to  reach  a  definite  decision.  This  is 
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only  possible  by  landerstanding  the  properties  vMch  describe  the  character  of 
the  material. 

Sane  of  these  properties  will  be  discussed  in  more  detail  on  the 
follcwing  pages.  These  properties  are:  mechanical,  physical,  chemical, 
thermal,  electrical,  and  optical.  In  addition,  availability,  cost, 
fabrication  techniques,  reliability,  and  maintainability  must  be  considered. 
Ihou^  not  often  thou^t  of  as  material  properties,  they  are  usually  based  on 
a  coribination  of  properties,  design,  and  economics. 

12.3.3.1  Mechanical  Properties.  Ihe  mechanical  properties  of  most  concern  at 
present  are  stiffiiess,  strength,  toughness,  ductility,  hardness,  and  fractxire 
toui^ess.  This  discussion  will  treat  only  tensile  loads,  but  the  same 
analysis  could  be  made  for  corpressive  and  shear  loads. 

Stiffness.  When  stiffness  is  considered,  there  are  two  definitions  that 
must  be  understood.  The  first  is  stiffness  as  it  refers  to  the  ratio  of 
stress  (a)  to  strain  (e)  in  the  elastic  region  of  the  stress-strain  diagram 
shown  in  Figure  12.2.  Stiffness  in  this  sense  is  Young's  Modulus  (E) ;  in 
torsion  the  analogous  parameter  is  G,  modulus  of  rigidity. 


FIGURE  12.2.  TYPICAL  STRESS-STRAIN  DIAGRAM 
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Stiffness  (E)  is  closely  related  to  the  bonding  between  atons  in  metals  and 
ceramics  and  the  bonding  (nunber  or  type)  between  chains  in  polymers.  For  a 
given  metal  or  ceramic,  there  is  little  that  can  be  done  to  significantly 
change  the  value  of  E  for  the  material.  In  the  case  of  polymers,  the  value  of 
E  can  be  changed  because  interchain  bonding  can  be  varied. 

Stiffness  (E)  is  not  necessarily  related  to  the  strength  of  a  material 
which  will  be  seen  when  strength  is  considered.  Figure  12.3  shows  vAiat  is 
meant  by  this  statement.  Material  A  has  the  highest  E  but  the  yield  strength 
is  less  than  for  Material  B  which  has  a  lower  E.  Ihe  meaning  of  yield 
strength  will  be  discussed  in  the  next  section. 


FIGURE  12.3.  STRESS-STRAIN  DIAGRAM,  DIFFERENT  MATERIALS 

Strength.  The  strength  of  a  material  is  a  mechanical  property  of  the 
material.  It  also  depends  on  vdiat  the  failure  criteria  are  defined  to  be. 
There  are  several  "strengths"  depending  on  vdiat  constitutes  failure:  (1) 

yield  strer^th  (o^) ;  (2)  ultimate  strength  (a^) ;  and  (3)  rupture  or  breaking 
strength.  The  cotmon  failure  criteria  are: 

1 .  Excessive  elastic  deformation  —  related  to  0^  since  when  stress  is 
equal  to  Oy,  plastic  deformation  begins. 
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2.  Excessive  plastic  defbnnation  -  an  exanple  of  this  is  creep  in  jet 
engine  tiorbine  blades. 

3.  Fractvire  -  separation  of  the  material  due  to  ductile  or  brittle 
fracture,  fetigue,  or  creep.  There  are  other  failure  criteria  such 
as  loss  of  appearance,  vAiich  often  must  be  considered.  The  three 
listed  above  are  of  major  concern  and  will  be  developed  further. 

Strength  data  for  a  given  material  are  usually  obtained  frcm  sirtple 
tensile  tests  of  standard  ^cimens  and  are  displayed  on  a  stress-strain 
diagram  as  shown  in  Figure  12.2.  Normally,  materials  display  a  region  of 
elastic  re^»nse  and  one  of  plastic  response.  The  response  is  taken  to  be 
elastic  (ccxnpletely  recoverable  under  ideal  conditions)  to  the  yield 
strength  (a^) .  Thereafter,  the  response  is  plastic,  and  the  deformation  is 
not  ccnplet^y  recoverable  vdien  the  load  is  removed  -  sane  permanait  set  being 
retained.  The  ultimate  strength  (a^)  at  Point  B  represents  the  ratio  of  the 
maximum  load  applied  to  the  specimen  to  the  original  specimen '  s 
cross-sectional  area  (A^) .  Fran  Point  B  to  C,  the  specimen  continues  to  be 
loaded.  However,  the  load  is  decreasing  in  magnitude  and  the  cross-sectional 
area  is  decreasing.  Consequently,  the  a/e  curve  decreases.  Point  C 
represents  the  complete  separation  of  the  material  (fracture) . 

Figure  12.2  showed  an  engineering  o/c  curve  in  vMch  all  loads  were 
divided  by  Aq.  Engineering  stress  is 


and  engineering  strain  (^gjqQ)  is  given  by 


^ENG 


(12.1) 


(12.2) 


However,  during  the  test,  the  cross-sectional  area  clearly  decreases  and  the 
length  increases.  If  instantaneoxis  values  of  cross-sectional  area  (A^)  and 
length  are  i:ised  rather  than  the  original  values,  true  stress  and  strain 
e^)  cure  obtained.  These  equations  are  '' 
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(12.3) 


and 

^TR 

but 

^ENG  and  can  be  related  by 

e^R  =  to  (1  + 

Figure  12.4  shows  a  ccnparison  of  true  and  engineering  stress-strain  diagrams. 
Ihe  true  o/e  curve  continues  to  increase  even  after  the  load  decreases  because 
the  instantaneous  cross-sectional  area  (A^)  is  decreasing  at  a  faster  rate  due 
to  necking  of  the  specimen.  This  discussion  will  consider  engineering  stress 
and  strain  since  the  engineer  is  concerned  with  the  original  areas,  lengths, 
etc. ,  in  the  design  phase. 

^^aterials  are  often  characterized  as  ductile  or  brittle  depending  on  the 
extent  of  their  plastic  behavior.  Figure  12.5  shows  hypothetical  materials  of 
varying  ductility.  Material  A  is  entirely  brittle  and  fractures  while  still 
responding  elastically.  Ceramics  display  this  type  of  behavior.  Material  B 
undergoes  some  plastic  deformation  before  fracture  and  is  characterized  as 
being  ductile.  The  strain  in  the  plastic  region  of  the  c/e  cur\^e  is  a  measure 
of  the  ductility  of  the  material.  If  the  material  is  greatly  strained  before 
fracture,  the  material  is  considered  to  be  ductile.  Ductile  and  brittle  are 
relative  terms;  thus.  Material  B  is  ductile  relative  to  Material  A  and  brittle 
relative  to  Material  C.  Uiere  are  three  additional  mechanical  properties  of 
great  concern  that  are  related  to  ductility:  hardness,  toughness,  and  fracture 
tou^mess. 
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FIGURE  12.5.  EFFECT  OF  DUCTILITy  ON  STRESS-STRAIN  CHARACTERISTICS 


Ifardness.  Ihere  are  a  number  of  ways  to  measure  the  hardness  of  a 
material:  abrasion,  r^und,  cutting,  indentation,  etc.  The  most  cannon 
method  for  metals  is  by  indentation.  Benetrators  of  specified  geometry  and 
material  are  forced  into  the  metal  by  known  loads  and  the  depth  of  penetration 
measured.  The  depth  of  penetration  is  then  a  measure  of  the  metal's  hardness. 
Since  the  penetrator  is  forced  into  the  metal,  the  metal  is  plastically 
deformed  indicating  that  hardness  and  ductility  are  directly  related.  If  the 
material  is  ductile,  it  can  be  easily  deformed  plasticeilly  and  consequently  is 
soft.  If  the  material  is  brittle,  it  cannot  be  easily  deformed  plastically 
and  is  therefore  a  hard  material.  Hardness,  then,  is  a  measure  of  the 
resistance  of  a  material  to  indentation. 

For  many  metals  (particularly  steels) ,  there  is  a  direct  relationship 
between  the  hardness  and  the  ultimate  or  tensile  strength  of  the  material. 
Uie  approximate  relationship  for  carbon  and  low  alloy  steels  is 

=  480  (BHN)  (12.6) 

vhere  BHN  is  the  Brinell  Hardness  Number  —  the  result  of  a  Brinell  Hardness 
Test.  Other  hardness  numbers  could  be  obtained  from  Rockwell  or  Vickers 
Hardness  Tests  and  then  related  to  the  BHN. 

Tou^mess.  lousiness  is  a  measure  of  the  ability  of  a  material  to  absorb 
energy  in  the  plastic  range.  One  way  to  measure  toughness  would  be  by  using 
the  area  under  the  a/e  diagram.  Thus,  in  Figure  12.5,  Material  C  would  have 
to  be  the  toui^est  and  Material  A  the  least  tou^. 

In  the  elastic  range,  there  is  a  property  related  to  toughness  vhich  is 
called  resilience.  Resilience  is  a  measure  of  the  capacity  to  absorb  energy 
in  the  elastic  range  and  is  the  area  under  the  elastic  portion  of  the  a/s 
curve.  Resilience  is  inportant  vhen  considering  such  things  as  leaf  springs 
for  automobiles.  Ihe  material  must  always  operate  within  the  elastic  range 
and  absorb  energy  in  that  range. 

There  are  other  properties  of  importance  such  as  electrical,  thermal, 
magnetic,  optical,  chanical,  and  physical.  Only  electrical,  thermal,  and 
chemical  propserties  will  be  briefly  mentioned  here. 

12.3.3.2  Electrical.  Some  materials  are  conductors  and  others  insulators; 
semiconductors  fall  in  between.  Ihe  electrical  conductivity  of  a  TnatP>-rial 
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depends  on  the  type  of  bonding.  If  the  valence  electrons  are  bound  as  they 
are  in  ionic  and  covalent  bonding,  then  the  material  is  an  ins\alator.  If  the 
valence  electrons  are  free  to  move  as  in  the  metallic  bond,  the  material  is  a 
conductor.  At  very  high  taiperatures  the  ions  in  ceramics  may  became 
electrical  conductors. 

12.3.3.3  Iherroal.  Thermal  conductiviiy  and  the  coefficient  of  thermal 
e}5)ansion  are  also  related  to  bonding.  Thermal  conductivity  requires  that 
electrons  be  free  to  move  so  the  same  materials  that  conduct  electricity  are 
normally  thermal  conductors.  Those  that  are  electrical  insulators  are  thermal 
insulators.  The  coefficient  of  thermal  ejq)ansion  (a)  is  inversely  related  to 
the  strength  of  the  atonic  bonds.  Therefore,  as  bond  strength  increases,  a 
decreases.  Many  designs  are  subjected  to  a  thermal  environment,  and  thermal 
stresses  are  set  up  that  result  in  failure  of  the  part.  Another  thermal 
effect  is  creep  or  accelerated  plastic  deformation  that  only  beccraes  a  problem 
at  about  one-half  the  material's  melting  point  on  an  absolute  scale. 

12.3.3.4  Chemical.  Chemical  properties,  as  used  here,  refer  to  the  corrosion 
and  oxidation  properties  of  the  material. 

12.3.4  Conclusion 

The  design  process  is  a  rather  long  and  exhaustive  process  because  of  the 
great  number  of  factors  to  be  considered  and  the  interplay  of  mechanics  and 
materials.  Good  design  is  possible  only  if  the  designer  is  familiar  vdth 
materials  and  every  materials  engineer  is  familiar  vdth  the  design  process. 

12.3.5  Ccnposite  Materials 

12.3.5.1  Introduction  to  Conposite  Materials.  Early  man  had  a  variety  of 
natural  materials  such  as  wood,  clay,  stone,  copper,  and  iron  avciilable. 
Althou^  he  did  not  understand  the  basis  for  a  material's  behavior,  he 
selected  materials  that  met  his  needs.  The  early  Egyptians  are  an  exartple  of 
men  vhci,  not  having  a  natural  material  available  to  meet  their  needs,  conbined 
two  materials  to  form  a  corposite  that  overcame  their  limitations.  They 
cotbined  mud  with  straw  to  produce  reinforced  bricks  vhich  they  xised  to  build 
their  cities. 
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Hiis  cycle  continues  and  even  more  so  in  today's  hi^ly  technological 
world.  Design  ideas  limited  by  materials  lead  to  the  developnent  of  new,  more 
advanced  materials.  The  need  for  higher  performance  materials  to  overcome 
present  design  limitations  has  created  a  new  class  of  materials  called 
conposites.  A  oonposite  material  is  a  combination  of  two  or  more 
constituents  from  the  three  general  classifications  of  materials:  metals, 
ceramics  and/or  polymers.  CJoirposite  materials  were  created  to  provide 
materials  with  improved  mechanical  properties  such  as  stiffness,  strength,  and 
high  temperature  stability.  The  development  of  composite  materials  is  a 
.revolutionary  advancement  in  materials  technology.  Until  recently,  materials 
selection  vas  a  substitution  affair  \«here  the  question  was,  "What  material  can 
I  use  in  place  of  alloy  X?"  The  material  selection  was  often  a  secondary 
decision  that  did  not  influence  the  design.  A  composite  material  allows  the 
materials  engineer  to  incorporate  the  required  prop)erties  vhich  greatly  affect 
the  design  and  should  be  considered  early  in  the  preliminary  design  phase. 

The  importance  of  composites  becoames  clear  if  one  examines  material 
properties  on  a  density-corrected  basis.  Table  12.1  illustrates  this  fact  by 
comparing  mechanical  properties  (density  corrected)  to  show  the  specific 
paroperties.  It  is  apparent  that  very  little  flexibility  in  material  selection 
actually  exists  for  conventional  materials.  Composite  materials  comhine  a 
reinforcing  material  with  high  specific  strength  and  specific  stiffness  with  a 
low  densi'ty'  matrix  material  that  results  in  a  system  that  is  synergistic 
(paroperties  of  the  vhole  are  greater  than  the  sum  of  the  individual 
constituents) . 

These  matericils  are  especially  useful  in  aerospace  applications  vhere 
veight  is  an  extremely  important  and  driving  design  parameter.  Development  of 
composites  has  been  primarily  financed  by  research  cind  develcptient  in  the 
aerospace  industry  conducted  by  both  private  industry  and  government,  A  large 
variety  of  laminated  composite  structures  have  been  designed  and  flight 
tested.  A  large  nimriber  of  materials  can  be  classified  as  composite  materials; 
however,  this  chapter  will  be  confined  to  the  more  advanced  fiber  reinforced 
composites  that  utilize  continuous  and  discontinuous  fibers  as  the  reinforce¬ 
ment  material. 
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Ti®LE  12.1 


MBOENICAL  PROPERTIES 


3 

Materials  p (lb/ in  ) 

E  (psixlO^) 

E/p  (inxlO*^) 

TS  (psixlO"^) 

Ductility  % 

Steel 

(Hot  Rolled) 
(Metal) 

.283 

30 

105 

65 

30 

Steel 

(Cold  Rolled) 
(Metal) 

.283 

30 

105 

85 

20 

Alumintm 

(Metal) 

.100 

10 

100 

63 

7 

Titanitm 

(Metal) 

.165 

17 

103 

165 

15 

Beryllium 

(Metal) 

.067 

42 

625 

83 

16 

Pine  Wood 
(Polymar) 

.0145 

1.3 

90 

7 

— 

Silica  Glass 
(Ceramic) 

.079 

10 

126 

— 

— 

Graphite 

Whiskers 

.060 

102 

1700 

2845 

— 

Boron/w 

Wire 

.095 

55 

580 

400 

— 

Boron  Epojqr 
(Ccarposite) 

.075 

20 

270 

133 

12.3.5.2  Continuous  Fiber  Iteinforceannent  Materials.  A  ccnposite  material  is 
composed  of  two  basic  conponents  -  a  reinforcing  fiber  vAiich  is  surrounded  by 
a  soft  light-weight  matrix.  The  majority  of  advanced  composite  materials 
utilize  continuous  fiber  reinforced  materials  fabricated  from  ceramics  and 
metals.  The  fiber  must  have  a  high  specific  strength  to  carry  the  load  in  the 
composite  structure  and  a  hi^  specific  stiffness  to  sujply  the  required 
structural  rigidity.  Designs  require  the  fibers  to  be  miformly  aligned  to 
maximize  the  contribution  of  fiber  strength  properties  to  the  composite 
structure.  Ideally,  the  fiber  axis  shotiLd  be  aligned  parallel  to  the  maximum 
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loading  axis  of  the  structtire.  Since  a  single  row  of  fibers  is  called  a 
lamina  (Figure  12.6)  and  one  must  stack  up  lamina  in  a  sequence  to  form  a 
laminated  structxire  (Figure  12.7)  vhere  the  fiber  orientation  can  be  varied, 
the  materials  engineer  can  design  in  a  wide  variety  of  material  properties. 
Another  benefit  the  fiber  offers  is  that  it  will  prevent  crack  propagation 
throu^  the  matrix. 

The  strength  and  stiffness  requirements  for  fibers  limits  the  material 
selection  primarily  to  ceramics  and  metals.  Table  12.2  shows  the  primary 
materials  being  investigated;  for  metals,  boron,  tungsten  and  steel  are  the 
most  promising.  For  ceramics,  graphite,  glass  and  silicon  carbide  hold  great 
prcmxse.  Boron  fibers  are  the  most  versatile  metal  fibers.  Research  and 
developnent  have  concentrated  on  these  fibers.  Tungsten  and  steel  are  of 
interest  for  hi^  taooperature  applications.  Graphite  is  the  most  attractive 
ceramic  fiber  and  it  can  be  made  into  a  variety  of  shapes  and  lengths. 
Its  primary  disadvantage  is  cost.  Glass  is  cheaper  by  comparison  and  does 
have  seme  very  attractive  properties;  however,  it  requires  special  surface 
coating  protections  called  sizes  to  protect  the  fiber  frem  degradation  by  the 
environment  and  from  handling  during  manufacture.  Ceramic  fibers  exhibit  the 
greatest  strengths  and  are  very  stable,  but  they  can  be  difficult  to  bond  and 
usually  require  a  coating  to  promote  better  interfacial  bonding.  The  metals 
tend  to  exhibit  lower  strengths  but  are  easier  to  bond.  However,  inter  facial 
reactions  can  degrade  the  fiber  properties. 


+ 

LAMINA 

THICKNESS 


o  o  o  o  o  o  o 
oooo  oooo 
o  o  o  o  o  o  o 


4 


LAMINA  (ONE  LAYER) 


FIGURE  12.7. 
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TABLE  12.2 


PROPERTIES  OF  FIBER  AND  MATRIX  MATERIALS 


Strength 

Modulus 

Melt  Point 

Density 

Tensile  St. 

Density 

E 

Density 

(°F) 

.  (Ib/in^) 

(Ksi) 

(xlO^)  (xlO^  Psi) (xlO^  in) 

AMORPHOUS 

FIBERS 

S-Glass 

1540 

.090 

650 

7.2 

12.6 

14.0 

Si02 

3020 

.079 

850 

10.8 

10.5 

13.3 

DISCCNTINUOUS  FIBERS 

(SINGLE  CRYSTAL) 

3700 

.143 

3000 

21.4 

62 

43.4 

Sic 

4870 

.116 

3000 

26.1 

70 

60.8 

Graphite 

6600 

.060 

2845 

47.4 

102 

170 

POLYCRYSTALLINE  FIBERS 

M2O3 

3700 

.144 

300 

2.6 

25 

21.9 

Graphite 

6600 

.054 

200 

3.7 

30 

56.5 

Tungsten 

6150 

.697 

580 

.8 

59 

8.5 

Steel 

2550 

.280 

600 

2.1 

29 

10.3 

MULTIPHASE  FIBERS 

Boron/W 

4170 

.095 

400 

4.2 

55 

57.8 

SiC/W 

4870 

.148 

300 

2.0 

70 

43.3 

Cont.  Tensile  , 

Use  Tern.  (^C)  Elongation  (%)  Strength  (Ksi)  Modvilus  (xlO”  Psi) 


MATRIX  MATERIALS 


Polycarbonate 

120 

80 

9 

.34 

Polysulfone 

150 

50 

10 

.36 

Epcxy-BPA 

145 

4.8 

13 

.36 

Epoxy-CA 

150 

3.5 

19 

.78 

Phenolic 

160 

2 

9 

.45 
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The  fibers  can  esdiibit  a  variei^  of  crystalline  structures  viiich  affect 
their  properties.  Single  crystal  discontinuous  fibers  are  of  an  order  of 
magnitude  stronger  than  all  the  other  fibers;  however,  they  are 
microscopically  small  and  difficult  to  fabricate.  Polycrystalline  fibers  have 
lower  strengths  due  to  the  inperfections  present  in  the  larger  fibers. 
Amorphous  fibers  are  strong  but  hi^ly  susc^tible  to  damage  and  degradation. 

Ttiere  are  a  variety  of  manufacturing  processes  used  to  fabricate 
reinforcement  fibers.  Glass  fibers  are  typically  prepared  by  making  the 
tnat-p>rial  molten  and  then  extruding  it  through  small  orifices  (s  .005  inch) 
vMch  removes  vater.  The  fibers  are  fired  at  a  tesiperature  to  densify.  Metal 
fibers  are  madp  using  the  wire  drawing  process;  however,  to  reduce  the  fiber 
down  to  .005  inch  requires  diamond  dies  which  are  very  expensive.  There  are 
other  processes  such  as  electronechanical ,  liquid  metal,  vapor  deposition,  and 
numerous  proprietary  processes. 

The  fabrication  of  boron  fibers  is  very  interesting  and  unique.  A 
tungsten  wire  .005  inch  in  diameter  is  passed  through  a  long  glass  tube  that 
contains  boron  gas  such  as  boron  trichloride  (BCL2) .  Pressure  is  maintained 
in  the  tubes  to  force  the  boron  to  deposit  on  the  tungsten  filament  vhich  is 
heated  electrically.  Ihe  entire  process  takes  less  than  two  minutes  and 
produces  a  filament  between  .003"  -  .005"  in  diameter.  The  time  can  be 
adjusted  to  yield  other  diameters  but  .008  inch  is  presently  the  upper  limit. 
The  fibers  are  usually  produced  in  continuous  lengths  of  10,000  feet.  The 
only  supplier  of  boron  fibers  in  quantity  is  the  AVCO  Corporation  in  lowell, 
Massachusetts.  One  advantage  of  composite  fibers  is  that  they  can  be 
processed  by  winding  around  mandrels.  They  are  then  woven  into 
two-dimensional  cloths  and  three-dimensional  mats  that  are  then  impregnated 
with  the  matrix.  Ihis  process  can  be  a  definite  advantage  in  producing 
coiplex  shapjes. 

12.3.5.3  Matrix  Materials.  Ihe  second  basic  corponent  material  is  the 
binding  or  matrix  material.  The  fibers  carry  load  and  give  structural 
rigidity.  The  matrix  material  bonds  the  fibers  and  holds  them  in  a  specific 
alignment.  It  transfers  loads  through  shear  mechanisms  to  the  fibers.  The 
matrix  material  must  be  li^t  wei^t  so  that  the'  benefits  of  the  specific 
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strength  and  specific  stif&iess  are  not  depleted.  Similarly,  a  soft  ductile 
material  vdll  assure  full  fiber  strength  is  utilized  and  provide  good  plastic 
behavior  for  transferring  load  by  shear  mechanisms.  The  matrix  also  serves  to 
protect  the  fibers  fron  the  environment  which  may  tend  to  attack  the  fibers 
and  degrade  properties.  This  is  especially  true  of  glass  fibers.  Uie  matrix 
will  arrest  crack  propagation  occinnring  in  fibers  since  it  acts  as  a  foreign 
material  interface  and  will  deflect  the  cracks.  Other  properties  the  matrix 
can  contribute  inclvide  toughness,  fatigue  strength,  oxidation  resistance, 
corrosion  resistance,  etc.  The  matrix  material  requirements  have  limited 
selection  to  polymers  and  metals.  Polymers  such  as  epoxies,  phenolics,  and 
polysulfones,  and  metals  such  as  titanium  and  aliiminum  possess  the  lower 
strength/hi^  ductility  behavior  that  is  necessary  in  a  matrix.  E^xies  and 
aluminum  have  been  extensively  researched  and  most  ocirponents  that  are  flying 
on  aircraft  or  are  in  the  development  stages  are  made  from  these  matrices. 
The  phenolics  and  titanim  are  still  in  their  early  stages  but  hold  great 
promise  for  higher  terrperature  applications.  Unlike  fibers,  these  materials 
can  iisually  be  fabricated  xising  conventional  techniques  alreacty  established 
vdiich  reduce  costs.  The  pol^ners  are  temperature  limited,  can  be  sensitive  to 
humidity,  and  possess  a  wide  variation  in  properties  such  as  chemical  and 
corrosion  resistance.  Important  considerations  in  design  are  the  stresses  and 
life  expectancies  required  of  these  matrix  materials.  These  polymers  also 
exhibit  shrinkage  vdiich  must  be  accounted  for  in  the  processing  since  it  can 
be  as  great  as  15%.  lypical  matrix  properties  are  shown  in  Table  12.2. 
12.3.5.4  Manufacturing  of  Conposite  Materials.  Careful  consideration  must  be 
given  to  the  decision  concerning  vdiich  fiber  to  carbine  with  vhich  matrix  to 
yield  the  optimum  design.  All  design  parameters  must  be  considered  in  order 
for  the  materials  engineer  to  design  in  the  required  properties.  This  further 
requires  that  material  selection  and  assessment  be  made  early  in  the 
preliminary  design  phase. 

Regardless  of  the  selection  of  fiber  and  matrix,  they  both  must  be 
properly  bonded  together.  Ihe  bonding  surface  between  fiber  and  matrix  is 
called  the  interfacial  bond,  and  it  has  an  enormous  effect  upon  conposite 
material  properties  and  fabrication.  The  metallxirgical  reactions  that  occur 
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are  cxJitplicated  and  difficult  to  control.  Poor  bonding  at  the  interface 
results  in  premature  failtires.  The  fiber  shape  is  also  irrportant.  The 
greater  the  surface  area  available  for  bonding,  the  better  the  bonding  will 
be.  Thvis,  non-circular  cross  sections  pronote  si^jerior  bonding. 

Qliere  are  a  number  of  techniques  ettployed  to  control  interfacial  reaction 
zones  to  limit  property  degradation  vAiile  pronoting  good  bonding.  The 
ceramics  tend  to  be  thermodynamically  stable  and  show  smaller  reaction  zones. 
They  do  not,  however,  bond  as  well  to  the  matrix  and  therefoie  require  a 
coating  such  as  polytieric  or  metallic  coating.  These  coatings  are  especially 
required  in  metal  matrix  cotposites.  When  boron  was  initially  used,  it  often 
had  to  be  coated  with  silicon  carbide  to  bond  it  with  aluminum  and  titanium. 
Advances  in  cciiposite  technology  and  fiber  treatment  have  eliminated  these 
requirements. 

Interface  reactions  result  in  the  destruction  of  the  fiber.  The  three 
most  prevalent  reactions  that  occur  include  diffusion  between  fiber  and  matrix 
vhich  can  leave  a  penetration  zone  with  recrystallization  that  causes 
premature  failure  of  the  fibers.  In  addition,  precipitation  reactions  and/or 
solid  solution  reactions  can  produce  fracture  initiation  sites  leading  to 
failure. 

A  variety  of  manufacturing  processes  are  utilized  that  allow  the  designer 
latitudes  which  conventional  materials  do  not  possess.  The  fibers  can  be 
woimd  around  a  mandrel  and  then  impregnated  with  the  resin  matrix  material  or 
a  "prepreg"  (partially  cured)  composite  can  be  wound  on  a  mandrel.  The 
"prepreg"  composite  is  often  used  to  allow  forming  of  lamina  into  desired 
shapes  and  then  using  a  mold  under  heat  and  pressure  to  produce  final  curing. 
The  metal  matrix  composites  can  be  handled  in  this  same  manner.  Often,  large 
sheets  of  composite  material  are  precut  into  ply  patterns  which  are 
hand-assembled  and  then  finish  processed.  Composites  can  be  woven  into 
two-dimensional  and  three-dimensional  shapes  vhich  are  then  impregnated  with 
the  matrix  and  cured.  This  type  of  fabrication  gives  the  designer  tremendous 
capabilities  that  otherwise  would  not  be  available. 

12.3.5.5  Design  Applications.  As  composites  become  more  common,  it  is 
necessary  for  engineers  accustomed  to  dealing  with  metals  to  become  more 
versatile  in  their  thinking.  The  various  properties  of  metals  are  generally 
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considered  to  be  isotropic  -  the  sane  in  all  directions.  The  design 
trade-offs  mentioned  earlier  are  thus  fairly  strai^tforward.  It  is  ccninon 
vAien  dealing  with  metals  to  consider  such  properties  as  modulus,  E,  as 
constants.  This,  in  fact,  is  not  the  case  with  conposite  materials. 
Cjonposites  are  anisotropic;  that  is,  having  different  strength,  stiffness  and 
thermal  properties  in  the  various  directions.  This  produces  an  added 
ccnplexity  for  the  designer  but  also  provides  him  with  more  flexibility  and  a 
chance  to  truly  optimize  a  structure.  If  the  designer  fully  xmderstands  the 
anisotropic  natiare  of  a  material,  he  is  able  to  design  the  material  to  meet 
the  specific  design  reguiranents  (strength,  stiffness,  etc.)  of  his  structure 
in  its  specific  directions.  This,  of  course,  forces  the  designer  to  analyze 
the  strengths  and  stiffnesses  required  in  his  structure  based  upon  such 
factors  as  the  es^jected  loads  and  allowable  deflections.  Since  structures  are 
often  loaded  in  only  one  or  two  directions,  a  ccrtposite  can  be  itade  strong  in 
those  directions.  This  can  be  compared  with  fabricating  steel  or  aluminum 
vhere,  in  order  to  meet  a  maximum  strength  requirement  in  one  direction,  the 
material,  because  of  its  basic  nature,  is  equally  strong  in  other  directions. 
It  can  thus  be  seen  that  a  composite  can  be  used  more  efficiently  and  produce 
a  second  weight  savings  in  addition  to  the  one  inherent  in  its  hi^ 
strength-to-weight  or  stiffness-to-^wei^t  ratios. 

12.3.5.6  Econcmic  Factors.  Composite  materials  have  gained  a  reputation  for 
being  ejqjensive.  To  seme  extent,  this  reputation  is  deserved.  As  graphite 
and  boron  filaments  were  first  being  developed  in  the  1960 's,  costs  of  $600 
per  pound  and  $400  per  poeind,  respectively,  were  common.  With  increased 
production,  these  costs  are  now  on  the  order  of  $100  per  pound  for  boron  and 
$50  per  pound  for  graphite,  vduLch  still  does  not  favorably  conpare  with  a  few 
dollars  per  pound  for  steel.  However,  the  cost  to  be  considered  is  not  the 
material  cost  but  the  cost  of  the  finished  product  or  the  dollars  per  payload 
of  the  final  system. 

The  cost  of  cotpositing  includes  the  tape  fabrication,  laying  the  tape 
into  layers  and  the  final  assembly  of  the  part  including  its  joining  to  other 
parts  of  the  structxnre.  Mental  parts,  of  course,  have  forming  and  machine 
costs  associated  with  them.  The  Navy  built  a  cenposite  wing  for  a  Target 
Drone  using  graphite  epoxy.  While  the  material  costs  were  hi^er  than  the 
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aluminum  it  replaced,  the  manxifacturing  steps  in  fabricating  the  wing  were 
reduced  fron  approximately  50  for  the  aluminum  to  3  for  the  graphite.  The 
savings  resulting  from  this  are  quite  obvious. 

C.E.  Cataldo,  in  an  article  on  the  use  of  ccnposites  for  the  space 
shuttle,  estimates  that  ccsiposites  can  save  5-9%  of  the  structural  weight  of 
approximately  200,000  pounds.  This  savings  of  10,000-18,000  pounds,  vdien 
compared  to  the  original  payload  of  25,000-65,000  pounds,  produces  a  p^load 
increase  of  up  to  40%.  The  cost  savings  resulting  from  this  are  so  huge  that 
they  overshadow  ary  material  costs. 

12.3.5.7  Analysis  of  Conposite  Materials.  A  technique  used  to  determine  the 
properties  of  a  composite  is  based  \:pon  the  volume  fraction  of  its 
constituents.  This  rule  of  mixtures  is  illustrated  in  the  following 
paragraphs. 

If  we  stipulate  that  the  fibers  are  uniformly  aligned,  as  shown  in  Figure 
12.7,  and  that  a  pjerfect  bond  exists  between  the  fiber  and  matrix,  then  the 
strain  is  the  same  for  the  conposite,  the  matrix,  and  the  fiber 

=  ^4:  - 

c  t  m 

The  total  load  in  the  conposite  will  be  the  sum  of  the  loads  carried  by  the 
fiber  and  the  matrix. 


Utilizing  the  axial  stress  relationship  and  solving  for  P,  the  equation  can  be 
rearranged  as  follows 

0  =  J  P  =  oA 

0  A  =  OxA-  +  0  A 
c  c  r  f  m  m 

Dividing  by  the  area  of  the  conposite  and  noting  it  is  equivalent  to  the  sum 
of  the  areas  of  the  fiber  and  matrix  yields 
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(12.7) 


=  V. 


f'^f 


+  V;nO„ 

m  m 


v^re  and  V  are  voliime  fractions.  Note  V _  +  v  =1.  Dividina  bv  the 
strain  yields  ™  f  m 


.  +  V  E 

c  f^f  m  m 


EXAMPLE:  Design  a  ccirposite  naterial  with  a  modul\as  of  elasticity  of  11 
GPa  and  a  yield  strength  of  160  MPa.  The  prop^ies  of  the  fiber  and  matrix 
are  shewn  belcw  and  have  been  measured  for  the  limiting  strain. 


=  14.  GPa 

E  = 
m 

7.  GPA 

=  210.  GPa 

0  = 
m 

105. MPa 

£  =  .015  m/m 

@  e 

m 

=  .015  m/m 

Consider  both  equations 


a  =  V,c-  +  Vo 
c  IX  mm 


160  =  210 (V^)  +  150  (V^) 


b.  E^  =  V,E.  +  V  E 
c  f  f  mm 


11  =  14  (V^)  +  7  (V^) 


I-tote;  =  1  - 

m  r 


160  =  210  (V^)  +  150  (1-V^) 


=  0.524 

a 


11  =  14  (Vj)  +  7  (1  -  Vj) 
Vj  =  0.571 
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NDtice  there  are  two  answers  and  in  order  to  meet  both  design  requirements  you 
must  choose  the  largest  value  where  the  volume  fraction  of  fibers  would  be 
57.1%. 

Ihis  type  of  analysis  is  often  applied  to  caiposites  and  is  adequate  for 
properties  such  as  density  and  certain  mechanical  properties  such  as  the 
longitudinal  stiffness.  Note,  however,  that  this  relation  says  that  the 
properties  of  a  coirposite  are  bounded  by  those  of  its  constituents. 

As  mentioned  earlier,  a  simple  analysis  such  as  the  rule  of  mixtures  or 
•idp.as  that  one  might  be  comfortable  with  in  dealing  with  metals  is  not 
adequate  in  the  analysis  and  design  of  ccnposite  structures.  Since  the 
mathematics  of  anisotropic  elasticity  and  of  ccnposite  laminate  analysis  are 
beyond  the  level  of  this  course,  this  chapter  will  simply  present  an  overview 
of  seme  of  the  techniques  used. 

There  are  more  steps  with  more  ccnplexity  in  the  design  of  a  structural 
corrposite  than  exist  in  the  design  of  metal  structures.  The  design  cycle 
starting  with  the  fiber  and  matrix  constituents  are  illustrated  in  Figure 
12.8.  Note  that  this  cycle  is  continuous.  This  implies  that  the  material's 
design  must  be  simultaneous  with  the  structure's  design.  The  final  objective 
of  the  design  process  is,  of  course,  to  achieve  an  optimum  structure.  This 
can  often  be  translated  into  ininiinum  wei^t. 

The  design  cycle  shown  in  Figure  12.8  combines  two  basic  approaches  vdiich 
have  been  labeled  microscopic  and  macroscopic.  The  microscopic  approach 
considers  the  mechanical  and  physical  characteristics  of  possible  fiber  and 
matrix  materials.  These  properties  are  normally  considered  to  be  isotropic; 
however,  the  structure  as  a  vhole  is  heterogeneous  (non-hemogeneous) .  Using 
various  assumptions  concerning  the  state  of  stress  and  strain  between  the 
fiber  and  matrix,  the  designer  must  select  the  proper  constituent  materials 
and  decide  upon  the  packing  arrangement  of  the  fibers  and  their  proper  volume 
fraction.  As  seen  in  Figure  12.9,  the  paroperties  of  the  unidirectional 
laminate  depend  upxDn  this  microscopic  analysis  of  the  constituents.  The 
reason  the  laminate  or  single  layer  of  composite  is  made  of  fibers  all 
oriented  in  one  direction  is  that  this  allows  ease  of  manufacture  and 
analysis. 

The  next  step  is  to  determine  the  mechanical  and  physical  characteristics 
of  the  building  block  laminate.  This  involves  vhat  is  called  a  macroscopic 


12.26 


analysis.  Here,  the  matrix  and  fiber  are  considered  smeared  together  to  form 
a  homogeno\3S  but  anisotropic  material  for  analysis  purposes.  Such  a  prooedxnre 
is  advantageous  since  existing  plate  and  shell  theories  can  be  applied  to  the 
analysis  of  these  thin  homogeneous  anisotropic  laminatae. 


FIGUEE  12.8.  CDMPOSITE  STRUCTURAL  DESIGN  AND  ANALYSIS  CYCLE 


Laminate  analysis  involves  determining  the  optimum  nirriber  and  orientation 
of  the  lamina  necessary  to  provide  the  desired  properties.  Presently  there  is 
not  a  unique  mathematical  solution  to  this  problem  and  the  ccntoination  of 
ccirputer  iteration  and  ejperience  guides  the  designer  to  a  final  laminate 
arrangement.  An  arrargement  of  0,  45,  90,  0  for  a  four-ply  structure  is 
illustrated  in  Figure  12.7.  Usually  the  optimization  procedure  produces 
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angles  vMch  might  be  difficult  to  manufacture;  thus,  it  is  comon  to  change  a 
calculated  0,  39,  87,  +5  arrangement  to  the  0,  45,  90,  0  mentioned  earlier  for 
ease  and  econorr^  of  manufecture. 

The  final  result  of  this  analysis  is  a  structure  of  lower  cost,  omproved 
maintainability  and  -increased  performance.  Table  12.3  shows  seme  of  the 
current  applications  of  cenposites.  During  development,  the  aerospace 
industry  tested  conposites  in  non-critical  structures  such  as  doors  in  the 
F-111  and  in  the  wing  of  a  Navy  Target  Drone.  The  use  of  materials  and 
fabrication  techniques  previoiKly  untried  was  often  done  to  the  reluctance  of 
both  designers  and  users.  Today,  however,  such  materials  are  seeing  use  in 
primary  aircraft  structures  such  as  the  longitiriinal  spine  of  the  B-1  and 
horizontal  and  vertical  stabilizer  sections  of  the  F-15. 

FIBER  VOLUME 


FIBER  STIFFNESS  (E,)  NORMALIZED  TO 
MATRIX  STIFFNESS  (E  J 

FIGURE  12.9.  TRRNSVERSE  STIFTNESS  IN  UNIDIRECTIC»IAL  LAMENME 

12.3.5.8  Discontinuous  Fiber  Reinforcement  Materials.  There  are  various 
types  of  discontinuous  fiber  reinforcement  materials.  In  this  paragraph, 
glass,  whisker  formation,  and  binary  eutectic  fiber  reinforced  coiposites  will 
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be  discussed.  The  discontinuous  glass  fiber  reinforced  catposites  are  the 
most  widely  used  ccxrposites  in  industry  tod^.  The  primary  advantages  of 
these  ccnposite  systems  are  in  their  manufacture.  They  utilize  the  technology 
applied  to  the  polymeric  and  textile  industries.  This  eliminates  the  need  to 
make  large  coital  investments  in  exotic  manufacturing  equipment.  In 
addition,  the  fabrication  of  the  glass  fibers  and  polymeric  matrices  are  well 
established  and  inexpensive.  The  glass  can  be  chopped  into  fibers  and 
randcraly  distributed  throughout  the  matrix  giving  three-dimensional  stiffening 
and  strengthening.  The  procedure  consists  of  mixing  a  polymeric  with  glass 
fibers  and  using  conventional  injection  molding  techniques  to  produce 
fiberglass  types  of  structures.  The  fibers  can  be  chopped  up  and  blown  into  a 
meshed  structure  with  a  matrix  added  to  form  a  "prepreg"  that  can  be  final 
cured  into  a  finished  product.  The  glass  fibers  can  be  woven  into  mats  and 
inpregnated  with  resin  to  form  a  continuovis  fiber  reinforced  ccnposite.  There 
are  mai^  epcxy  and  phenolic  resins  that  can  be  used  with  the  chopped  fibers 
but  these  composites  do  not  exhibit  hi^  strength  levels.  They  are  vised  for 
corrosion  resistance,  li^t  wei^t,  low  cost  manufacturing,  op>tical 
properties,  etc. 

It  is  obvious  fron  Table  12.2  that  the  single  crystal  vhiskers  p»ssess  an 
order  of  magnitude  greater  psropjerties  over  the  polycrystalline,  amorphous  and 
multiphase  fibers.  These  whiskers  are  free  frcm  defects  and  dislocations,  and 
their  properties  approach  theoretical  values.  The  obvious  question  is,  "Why 
aren't  these  coitposite  fibers  being  used  throui^out  the  industry?"  The  answer 
is  technically  complicated  but  economically  sirnple.  The  vhiskers  are  only 
100-200  microns  long  vhich  is  the  reason  they  can  be  fabricated  defect- free. 
Bit  they  are  also  difficult  to  bond  to  the  matrix  and  infinitely  more 
difficult  to  align.  The  handling  and  fabrication  of  these  vhiskers  requires 
ejpensive  and  sophisticated  equipment.  Economically  they  are  impractical.  A 
typical  exaitple  is  silicon  carbide  vhiskers.  They  are  grown  using  a 
sublimation  technique  that  involves  the  prolysis  of  organosilanes,  silicon 
conpounds,  and  hydrocarbons  in  a  lydrogen  atmosphere  at  1500-2000°C.  Graphite 
can  be  formed  into  scroll-type  fibers  by  ccxtbining  sublimation  techniques  with 
DC  arcing  under  hi^  pressure  in  an  inert  atmosphere.  Other  techniques  use 
evaporation  and  vapor-liquid-solid  deposition  processes,  but  all  are 
prohibitively  expensive. 
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TABLE  12.3 


COMPOSITE  APPLICATIONS  AND  DEVELOPMENT 


SPACE  SHUTTLE  APELICATiaJS 

llirust  Hot 


Structvire 

Fuselage 

Bulkheads 

Tanks 

Structure 

Shrouds 

SKINS 

Boron-Epoxy 

Ring  Frames 

Glass-Ip)oxy 

Glass-  Carbon-Carbon 

Filament 

Wound  with 

Liners 

Boron- 

Epoxy 

Boron-Alian. 

Longerons 

Graphite 

Epoxy 

Graphite 

Coated 

Boron- 

Alum. 

Graphite- 

Epoxy 

Stringers 

Graphite 

Epoxy 

Skins 

Beams 

Glass- 

Honeyccnb 

]^AMS 

Fittings 

Graphite- 
Epoxy  for 
small 
tanks 

Boron- 

Polvitnide 

650* 

Same 

Boron-Epoxy 

Graphite- 

Polyimide 

650° 

Boron-Alum 

Boron-Alum. 

600° 

STRUTS  Graphite-Epoxy  Metal-Ifetal 

2000-2500° 

Same  Glass-Epcc^ 

Borsic-Ti  800° 


12.30 


TJ®LE  12.3  (continued) 
COMPOSITE  APPLICATICNS  AND  DEVELOPMENT 
DEVELOPMENT  STATUS 


Current  Developments  Advanced  Developnent  Futvire  Development 

(Near  Term  Payoff)  ■  (Long-Term  Payoff)  (New  Concepts-New  Materials) 


Boron-Epoxy 


Boron-Polyimide  In-Situ  Reinforcements 

Graphite-Polyiraide 


Graphite-Epojy 

Boron-Aluminum 

Carbon-Carixjn 


Wire  Reinforced  Ifetals 

Graphite-Aluminum  High  Tertperature  Polmers 

Beryllium  Conposites  Low  Cost  High  Speed 

Production 


lybrids 


Advanced  I^brids 


COMPOSITE  TEMPERATURE  CAPABILITIES 
Maximum  Tenperature  Capjabilities 

Polymer  Matrices  Tenperature  (*^F)  Metal  Matrices 

S-Glass/Epoxy 


Boron/Epcxy 

350* 

Graphite/Epoxy 

Graphite/Polyimide 

600 

Boron/Alinninum 

650 

Graphite  /  Aluminum 

700 

800 

Borsic/Titanium 

Carbon /Carbon 

1000 

Beryllium/Titanium 

Carbon/Carbon 

(Uncoated) 

1800 

Borsic/Nickel 

Wire /Super alloys 

Carbon /Carbon 
(Anti-Oxidant) 

2400 

Molybdenum/Colunibium 

Tantalum 

*  Assuming  Stable  Care 
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Ihe  binaiy  eutec±ic  system  offers  the  only  real  viable  alternative  to 
this  problem  vMch  is  available  with  today's  technology.  A  binary  eutectic 
occurs  v^en  a  liquid  transforms  into  a  two-phase  solid  at  constant 
temperature.  Ihe  reaction  is  invariant  and  occurs  only  at  a  specific 
conposition  and  a  specific  tenperature.  It  is  a  specific  point  on  the  phase 
diagram.  It  should  be  noted  that  not  all  binary  systems  will  exhibit  this 
jiiencraenon  and  of  those  that  do,  not  all  can  be  controlled  to  produce  an 
aligned  microstructure. 

Ihe  directional  solidification  (DS)  of  the  eutectic  is  the  key  processing 
parameter.  It  is  the  directional  cooling  vdiich  permits  the  alignment  of  the 
vdiisker  fibers  within  the  matrix  material.  Controlling  the  cooling  gradient 
by  slowly  withdrawing  molten  material  from  the  furnace  orients  the  fiber 
growtii  direction  parallel  to  the  uniaxial  heat  flow.  Hie  fiber  and  matrix 
materials  freeze  into  a  solid  solution  forming  a  discontinuous  fiber 
reinforced  carposite  material  in  one  step.  The  growth  rates,  however,  are 
only  about  a  quarter  of  an  inch  per  hovur  vhich  means  the  manufacturing  costs 
are  quite  hi^.  The  primary  application  of  interest  is  in  the  turbine  of  an 
aircraft  engine.  The  enormous  property  increases  gained  and  the  potential 
engine  performance  irtprcvements  justify  the  high  manufacturing  costs.  The 
eutectic  system  receiving  the  most  attention  is  the  tantalum  carbide  fiber  in 
a  nickel  base  super  alloy  or  cobalt  based  alloy.  There  are  many  shapes  the 
fiber  can  take;  however,  it  has  been  shown  that  rods  or  platelets  give  the 
best  properties.  These  alloys  cure  very  stable.  They  compete  with 
conventional  alloys  at  room  taiperature  but  are  vastly  superior  at  high 
torperatures.  Their  properties  do  not  degrade  even  at  levels  close  to  their 
melting  point.  Since  the  vMskers  grow  naturally  within  the  matrix,  there  are 
no  bonding  problems.  Experiments  conducted  on  the  interfacial  bond  strength 
show  it  is  stronger  than  the  fiber  itself.  This  is  probably  due  to  the  unique 
matcliing  of  preferred  crystallographic  planes  that  occur  during  processing. 
Figure  12.10  center  shows  a  DS  turbine  blade.  The  DS  blade's  carposite 
structure  is  coaxpared  to  other  forming  techniques  —  a  conventional  casting 
left  and  a  single  crystal  turbine  blade  right. 

A  low  volume  fraction  of  fiber  will  give  a  rod-like  reinforcement 
structure,  vbereas  near  equal  fractions  of  matrix  and  fibers  produce  a 
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platelet  structure.  Hie  rod  fibers  are  superior  and  more  predictable; 
however,  they  only  supply  viniaxial  strengthening.  Ihe  platelets  have  a 
two-dimensional  strengthening  effect.  In  addition,  studies  on  solute  atoms, 
dislocation  mobility,  tenperature,  strength  and  plastic  properties  shew  that 
intermetcillic  corepounds  make  better  fibers  than  do  ceramics  or  metals.  These 
materials  offer  great  advancements  to  aircraft  and  engine  performance. 


FIGURE  12.10.  MICKBTRUCTURE  OF  CONVENTICHALLy  CAST  (LEFT)  , 
DIRECTIONALLY  SOLIDIFIED  (CENTER)  AND  SINGLE 
CRYSTAL  (RIGHT)  TURBINE  BLADES 
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GLOSSARY 


12.3.5.9 

Fiber  -  Iteinforcing  ccnponent  of  the  ccnposite  that  supplies  the  strergth 
and  stiffness  prc^rties. 

Matrix  -  Soft  li^t-wei^t  bonding  catponent  of  the  ccnposite  that  holds  the 
fibers. 

TaTTiina  -  A  single  layer  of  fibers  surroimded  by  matrix  material. 

TflTtiinate  -A  structure  build  up  using  lamina  at  various  orientations  to  yield  a 
structural  element. 

Prepreg  -  A  partially  cured  ccnposite  laminate  that  is  pliable  and  can  be 
shaped  and  then  finish  cured. 

Interfacial  Bond  -  Surface  bond  betsreen  the  fiber  and  the  matrix. 

Eutectic  composite  -  A  two-phase  solid  vAiere  one  phase  acts  as  fiber  rein¬ 
forcement  and  the  second  phase  is  the  matrix. 

Whisker  -  A  smel  1  single  crystal  fiber  that  is  defect  free  (length  is  100-200 
microns) . 
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12.4  FUNDAMENTS  OF  STRUCTORES 


This  section  on  the  fundamentals  of  structures  has  been  prepared  to 
provide  a  practical  conbination  of  theory  and  application  to  the  operational 
problems  of  aircraft/  missiles /  and  space  craft.  For  this  reason  only  the 
minimum  mathematical  relationships  are  provided. 

Ihere  are  three  fundamental  objectives  that  will  be  addressed; 

(1)  Acquaint  the  reader  with  the  basic  properties  of  structural 
materials  and  the  particular  qualities  of  these  materials  that  make 
them  suitable  in  particular  structures. 

(2)  Furnish  the  reader  the  fundamental  reasons  for  operating  strength 
limitations  and  good  maintenance  practices. 

(3)  Equip  the  reader  with  the  ability  to  recognize  and  diagnose  the 
causes  of  structural  and  mechanical  failures. 

These  objectives  are  first  served  by  describing  the  principal 
requirements  of  any  aircraft,  missile,  or  space  craft  structure.  Tie  most 
important  basic  requirement  is  that  the  primary  structure  should  be  of  the 
lowest  possible  veight.  All  of  the  basic  items  of  performance  and  efficiency 
of  a  configuration  are  seriously  affected  by  the  structural  weight.  This  is 
especially  true  when  the  extremes  of  performance  are  demanded  of  a 
configuration.  For  example,  during  preliminary  design  of  a  long  range  jet 
aircraft,  a  configuration  v^i^t  growth  factor  of  twenty  may  be  typical.  In 
other  words,  if  the  wei^t  of  ary  single  item,  e.g. ,  landing  gear  structure, 
vere  to  increase  one  pound,  the  gross  weight  of  the  aircraft  must  increase 
tweniy  pounds  to  maintain  the  same  performance.  Any  additional  weight  would 
require  more  fuel,  more  thrust,  larger  engines,  greater  wing  area,  larger 
landing  gear,  heavier  structure,  more  fuel,  etc.,  until  the  aircraft  gross 
weight  had  increased  tv^ty  times  the  original  weight  change. 

long  range  missiles  and  space  craft  usually  encounter  a  design  growth 
factor  viMch  is  considerably  in  excess  of  any  typical  aircraft.  Some  typical 
long  range  ballistic  missiles  have  demonstrated  preliminary  design  growth 
factors  on  the  order  of  80  to  200.  Of  course,  such  configurations 
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represent  an  extreme  of  performance  but  serve  notice  of  the  great  significance 
of  structural  weic^t.  A  limiting  situation  can  exist  vdien  the  demands  for 
performance  exceed  the  "state-of-the-art".  If  performance  demands  are  extreme 
and  basic  powerplant  capabilities  are  relatively  lew,  the  growth  factors 
approach  infinite  values  and  iirpractical  gross  weights  result  for  the 
configuration. 

While  idle  structure  must  be  of  the  lowest  possible  wei^t,  it  must  also 
be  easily  accessible  for  repair,  inspection,  and  maintenance.  There  must  be 
qdfypiate  protection  from  the  environment  to  prevent  corrosion,  ionizing 
radiation,  etc.  In  many  instances,  the  accommodation  for  easy  access  for 
siirple  maintenance  must  be  forsaken  siitply  to  obtain  reduced  structural 
wsight. 

Uie  primary  structure  must  be  the  minimum  weight  structure  vdiich  can 
safely  sustain  the  loads  typical  of  operation.  Bie  actual  nature  of  the  most 
critical  loads  will  depend  to  a  great  extent  upon  the  design  mission  of  the 
vehicle.  During  design  and  development,  the  mission  must  be  thoroughly 
analyzed  to  define  the  most  critical  loads  that  will  determine  the  minimum 
necessary  (size  and  weight)  of  the  struct^3ral  elements.  From  an  apparent 
infinite  number  of  possible  situations,  the  most  critical  conditions  must  be 
defined,  (fenerally,  there  are  three  important  areas  of  structural  design,  any 
one  of  which  (or  conbination)  could  provide  the  most  critical  requirements  of 
the  structure. 

12.4.1  Static  Strength  Considerations 

Static  loads  refer  to  those  loads  vdiich  are  gradually  applied  to  the 
structure.  The  effects  of  the  onset  of  loading  or  the  repsetition  of  loading 
deserve  sep)arate  consideration.  Throu^out  the  op)eration  of  its  mission,  a 
vehicle  structure  encounters  loads  of  all  sorts  and  magnitudes.  Various  loads 
may  originate  during  manufacture,  transport,  erection,  launch,  flight  gusts 
and  maneuvers,  landing,  etc.  These  various  conditions  may  be  encountered  at 
various  gross  wei^ts,  e.g. ,  positions,  altitudes,  pressurization,  etc.  If 
p)articular  elCTients  of  the  structure  are  separated  for  stuefy,  it  is 
appreciated  that  these  elements  are  si±)ject  to  a  great  spectrum  of  varying 
loads. 

For  the  considerations  of  static  strength,  it  is  important  that  this 
^ctrum  be  analyzed  to  select  the  maximum  of  all  loads  encountered  during 
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normal  (or  intended)  operation.  This  maxiicium  of  all  normal  service  loads  is 
given  ^cial  significance  by  assigning  the  nonraenclature  of  "limit"  load.  The 
specific  requirement  of  the  structure  is  that  it  must  be  able  to  vdthstand 
"limit"  load  vathDut  ill  effect.  Most  certainly  the  structure  must  withstand 
limit  load  without  c±)jectionable  permanent  deformation. 

Specific  requirements  are  different  for  various  structureil  applications, 
and,  in  sane  cases,  a  "yield"  factor  of  safety  of  1.15  must  be  incorporated. 
This  requirement  would  demand  that  the  primary  structure  be  capable  of 
withstanding  a  load  15%  greater  than  limit  without  "yielding"  or  deforming 
sane  objectionable  amount.  If  such  requironents  vrere  specified  for  a  fighter 
aircraft,  the  aircraft  could  be  safety  maneuvered  to  limit  "G"  without  causing 
the  aircraft  to  be  permanently  deformed.  If  such  requirements  were  specified 
for  a  typical  missile,  the  missile  could  be  fueled  and  static  tested  without 
causing  the  structure  to  be  permanently  deformed.  Of  coxjrse,  the  nxrnber  of 
times  this  action  could  be  repeated  without  ill  effect  would  not  be  part  of 
the  static  strength  consideration. 

A  separate  provision  must  be  made  to  account  for  the  possibility  of  a 
one-time  application  of  some  severe  load  greater  than  limit.  For  exaiiple,  the 
previously  mentioned  filter  aircraft  may  require  sane  fli^t  maneuver  load 
greater  than  limit  in  order  to  avert  a  disaster  of  collision.  The  same  idea 
applies  to  the  missile  vAiere  malfunction  of  equipment  may  cause  higher  than 
normal  tank  paressurization.  In  either  of  these  exanples,  seme  load  greater 
than  limit  is  always  a  (remote)  pxDssibility  and,  within  reasonable  limits, 
should  not  cause  a  catastrophic  failure  of  the  primary  structiure.  There  must 
be  some  provision  for  the  possibility  of  a  single  critical  load  greater  than 
limit. 

Experience  with  piloted  aircraft  has  shown  that  an  "ultimate"  factor  of 
safety  of  1.5  is  satisfactory.  Thus,  a  primary  structural  element  should  be 
capable  of  withstanding  one  load  50%  greater  than  limit  without 
failvure.  Of  course,  loads  vdiich  generate  stresses  greater  than  the  "yield" 
point  will  cause  objectionable  psermanent  deformation  of  the  structure  and 
render  it  unsuitable  for  continued  operation.  The  principal  concern  is  that 
the  primary  structure  withstand  the  "ultimate"  load  without  failure.  To  be 
sure,  the  ultimate  load  can  be  resisted  only  by  a  sound  structure,  i.e. ,  no 
cracks,  corrosion,  eccentricity,  etc. 


12.37 


In  the  previoijs  discussion  the  yield  factor  of  safety  of  1.15  and 
ultimate  factor  of  safety  of  1.5  have  been  selected  for  exartple  since  these 
values  are  representative  of  piloted  aircraft.  On  the  other  hand,  certain 
missile  configurations  may  have  factors  of  safety  well  below  that  of  piloted 
aircraft,  e.g.,  yield  factor  of  safety  of  1.0,  ultimate  factor  of  safety  of 
1.2.  In  order  to  ccraplete  the  picture,  seme  ground  support  equipment  may  have 
an  ultimate  factor  of  safety'  of  five  or  six,  and  a  typical  bridge  structure 
may  have  an  ultimate  factor  of  safety  of  twenty.  The  factors  of  safety  for 
airborne  vehicles  must  be  as  low  as  is  consistent  with  the  safety  and 
integrity  of  the  structure. 

When  specific  limit  loads,  yield  and  ultimate  factors  of  safety  are 
defined,  there  will  be  no  deliberate  addition  of  strength  above  these 
^jecified  minimums.  Ihe  reason  is  sinple:  undesirable  structural  weight 
would  be  added.  However,  the  normal  variation  of  material  strength  properties 
must  be  accounted  for  by  designing  to  minimum  guaranteed  strength  or  specific 
levels  of  probability.  In  either  case,  it  is  possible  that  a  considerable 
percentage  of  the  structures  will  exceed  the  strength  requirements  by  slight 
margins.  This  is  an  e:q)ected  result  vben  the  structure  is  required  to  meet  or 
exceed  the  minimum  specified  values. 

As  a  result  of  these  static  strength  considerations,  the  primary  structure 
must  withstand  limit  load  without  objectionable  permanent  defontation  and 
ultimate  load  without  failure.  Because  of  the  yield  factor  of  safety  and 
certain  material  characteristics,  objectionable  permanent  deformation  does  not 
necessarily  take  place  inmediately  above  limit  load.  This  could  lead  to 
difficulty  in  appreciating  over-stress  conditions  since  objectionable 
permanent  deformation  does  not  necessarily  occur  just  beyond  limit  load. 
However,  at  ultimate  load,  failure  is  imminent. 

12.4.2  Rigidity  and  Stiffness  Considerations. 

"Strength"  coiiLd  be  defined  as  the  resistance  to  applied  loads.  On  the 
other  hand,  "stiffeess"  coiiLd  be  defined  as  the  resistance  to  applied 
deflections.  This  particular  distinction  between  strength  and  stiffness  is  a 
necessary  consideration  since  the  development  of  adequate  strength  does  not 
enstire  the  attainment  of  adequate  stiffness  and  rigidity.  In  fact,  the 
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particular  requirements  of  stiffness  must  be  given  consideration  vAiich  is 
separate  (but  not  ccrpletely  independent)  fron  the  basic  strength 
considerations . 

The  stiffness  characteristics  of  a  structure  are  very  important  in 
defining  the  re^nse  of  the  structure  to  dynamic  loads.  In  order  to 
distinguish  a  c^amic  load  fron  the  ordinary  static  load,  inspect  Figure  12.11 
where  a  weight,  W,  is  suspended  by  a  spring  which  has  a  stiffness,  K. 

GRADUALLY  APPLIED  LOADING  SUDDEN  LOADING 


EQUILIBRIUM 

FIGUKE  12.11.  TIME  DEPENDENCE  OF  LOAD  APPLICATION 

If  the  weight  ;^re  lowered  slcwly  and  the  force  gradually  applied  to  the 
spring,  the  spring  would  deflect  slowly  until  the  spring  supports  the  entire 
veight  with  an  equilibrium  deflection  of  6  =  W/K.  Alternatively,  if  the 
weight  is  suddenly  dropped  onto  the  spring,  the  input  energy  of  the  sudden 
loading  will  cause  the  spring  to  deflect  twice  as  greatly,  26.  Then  the 
weight  will  oscillate  back  and  forth,  finally  caning  to  rest  at  the  same 
equilite’ium  deflection  as  for  the  gradually  applied  load.  During  the  first 
plunge  vten  the  spring  is  deflected  26,  the  spring  is  subject  to  an 
instantaneous  load  vrfiich  is  twice  the  weight.  Under  the  dynamic  loading 
illustrated,  the  <^amic  load  is  twice  as  great  at  the  static  load.  Of 
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coiirse,  if  the  wei^t  had  been  projected  da-mward  onto  the  spring  with  an 
initial  velocity,  the  input  energy  would  be  greater  and  the  dynamic 
airplification  of  load  woiild  be  increased  considerably. 

The  previous  exaitple  serves  only  to  point  out  the  serious  nature  of 
dynamic  loads.  In  a  more  typical  (and  cotplex)  aircraft  or  missile  structure, 
itany  degrees  of  freedcm  exist  and  the  response  may  shew  the  coupling  between 
various  possible  modes  of  oscillation.  In  any  case,  the  energy  of  load  input, 
the  rate  of  onset,  and  the  characteristic  response  of  the  structure  must  be 
examined  to  determine  the  critical  airplification  of  loads.  The  stiffness 
characteristics  and  the  existence  of  danping,  either  natural  or  synthetic, 
must  be  tailored  -  if  possible  -  to  minimize  critical  airplification  of  loads. 

Vibration  of  structures  may  provide  critical  situations  and  create  a 
source  of  damaging  loads.  Ihe  fundamental  nature  of  a  vibrating  system  is 
best  illustrated  by  the  simple  spring-mass  system  of  Figure  12.12.  If  a 
wei^t,  W,  is  suspended  on  a  spring  of  stiffness,  K,  a  small  disturbance  of 
the  system  will  cause  the  weight  to  oscillate  at  seme  frequency  vhich  is 
dependent  on  spring  stiffness  and  weight. 


PLATFORM  EXCITATION 
IS 

asin  27r  f(t) 


SPRING-MASS  NATURAL  FREQUENCY 
IS 

•n  271  V  W 


CYCLIC  AMPLITUDE  OF 
OSCILLATING  WEIGHT  IS  A 


FIGUFE  12.12.  FUNDAMENTAL  NATURE  OF  A  VIBRATING  SYSTEM 
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This  "natviral"  frequency  is  related  by  Equation  12.9. 


(12.9) 


where  f^=  natural  frequency,  cps 
K  =  spring  stiffness,  Ibs/w. 

g  =  acceleration  due  to  gravity,  386 

W  =  v^ight,  lbs. 

In  order  to  consider  the’  possibili-ty  of  a  forced  vibration  of  this 
system,  suppose  the  platform  is  moved  back  and  forth  with  a  sinusoidal  motion, 
a  sin  2  irf  (t) ,  v^ere  f  is  the  frequency  of  excitation  and  a  is  the  airplit'ude 
of  excitation.  The  cyclic  motion  of  the  platform  will  induce  a  cyclic  motion 
of  suspended  veight.  The  airplitude  of  motion  of  the  v^ight.  A,  is  related  to 
the  platform  excitation  and  the  natural  frequency  of  the  spring-mass  system. 
Equation  12.10  defines  this  relationship; 


A 


1  - 


(f/f 
'  /  n 


(12.10) 


where  A  =  airplitude  of  vreight  motion 

a  =  airplitude  of  platform  motion 
f  =  excitation  frequency  of  platform 
f^  =  natural  frequency  of  spring-mass  system 

A  careful  inspection  of  Equation  12.10  points  out  one  of  -tbe  imdesirable 
possibilities  of  a  forced  vibraticxi  of  a  structure.  As  the  excitation 
frequency  of  the  platform  equals  the  natural  frequency  of  the  system  (f/f^=l) , 
a  resonant  oondition  develops  and  the  amplitude  of  vreight  motion  approaches  an 
infiojite  value.  Of  course,  this  resonant  condition  could  cause  sudden  failure 
of  the  spring. 

If  the  platform  is  subjected  to  an  excitation  frequency  which  is  well 
below  the  natural  frequency,  the  weight  anplitude  is  very  nearly  the  same  as 
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platform  anplitude.  TJio  vJBight  would  move  along  with  the  platform  with 
only  slightly  greater  than  static  deflection  of  the  spring.  VJhen  the 
excitation  firequency  is  considerably  greater  than  the  natural  frequency ,  the 
weight  is  essentially  isolated  v^ile  the  platform  oscillates.  Uie  cyclic 
deflection  of  the  spring  would  approach  the  cyclic  displacement  of  the 
oscillating  platform.  If  dairying,  or  resistance  to  motion,  is  introduced  into 
the  system,  the  resonant  condition  will  sirtply  produce  less  than  infinite 
motion  of  the  oscillating  wei^t.  Above  and  below  the  resonant  condition, 
daitping  will  alter  the  motion  depending  on  the  amount  of  danping  present  in 
the  system. 

labile  the  simple  ^tem  illustrated  does  not  portray  the  bdiavior  of 
cotplex  structures,  the  fundamental  relationships  are  the  same.  During  design 
of  a  structure  the  stiffness  and  mass  distribution  must  be  tailored  to  ensure 
that  the  ordinary  environment  of  vibration  does  not  allow  any  approach  to 
resonant  conditions  and  produce  damaging  loads. 

Aeroelastic  problems  are  encoxmtered  due  to  the  interaction  of 
aerodynamic  forces  and  elastic  deflection  of  the  structures.  Since  elastic 
deflections  are  involved,  the  inherent  stiffness  and  rigidity  of  the  structure 
is  a  principal  qualil^  determining  the  extent  of  aeroelastic  problems. 

Static  aeroelastic  problens  involve  only  the  relationship  of  the 
aerodynamic  forces  and  elastic  deflections  without  the  generation  of  inertia 
forces.  A  typical  static  aeroelastic  problem  encountered  in  aircraft  is  the 
phencmenon  of  "aileron  reversal".  Deflection  of  an  aileron  produces  a  section 
pitching  mcment  tending  to  twist  the  wing  in  torsion.  Thus,  if  an  aileron  is 
deflected  down  at  hi^  speed,  the  wing  may  develop  such  significant  twisting 
deflection  that  the  aircraft  may  roll  opposite  to  the  direction  desired.  Of 
course,  sufficient  stiffness  must  be  provided  in  the  structure  to  prevent 
aileron  reversal  or  any  significant  loss  of  control  effectiveness  within  the 
intended  range  of  flight  speeds. 

A  more  disastrous  sort  of  aeroelastic  problem  is  referred  to  as 
"divergence."  Suppose  that  a  surface  is  subjected  to  a  slight  gust  vdien 
at  very  hi^  speed.  If  the  change  in  lift  occurs  fon^d  of  the  elastic  center, 
the  surface  will  tend  to  twist  leading  edge  up  as  well  as  bend  up.  The  twist 
represents  additional  angle  of  attack,  more  lift,  more  bend,  more  twist,  etc. , 
until  a  sudden  failure  results.  Such  a  failure  is  sudden  and  catastrophic 
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without  warning.  It  is  obvious  that  divergence  could  not  be  tolerated  and 
sufficient  stiffness  must  be  present  to  prevent  divergence  within  the 
anticipated  flight  range.  In  addition,  below  the  divergence  speed  there  must 
be  sufficient  stiffness  to  ensure  no  serious  change  in  load  distribution  due 
to  this  sort  of  interaction  between  aerodiTiamic  forces  and  elastic 
deflections. 

All  of  the  static  aeroelastic  problems  are  specific  to  the  stiffness 
qualities  of  the  structure  and  the  <^niamic  pressure  of  flight.  Thus,  any 
specific  operating  limitation  iirposed  will  be  relative  to  a  certain  dynamic 
pressure  hence,  indicated,  calibrated,  or  equivalent  airspeed. 

Ihe  dynamic  or  oscillatory  aeroelastic  problons  introduce  an  additional 
variable;  inertia  forces.  Thus,  ^mamic  aeroelastic  problems  involve  some 
combination  of  aeroc^^arnic  forces,  elastic  deflections,  and  inertia  forces. 
"Flutter"  is  one  such  problem.  If  a  surface  with  particular  mass  and 
stiffness  distribution  were  exposed  to  an  airstream,  the  oscillatory 
aeroc^amic  forces  may  combine  with  the  various  natural  oscillatory  modes  of 
the  surface  to  produce  an  unstable  motion.  Flutter  is  essentially  an 
aero<^amically  excited  oscillation  in  vMch  airstream  energy  is  extracted  to 
amplify  the  eneigy  of  the  stixictural  oscillation. 

Flutter  is  not  necessarily  limited  to  control  surfaces  or  wing  surfaces. 
Structural  panels  m^  encounter  flutter  conditions  vdiich  are  just  as  critical 
and  damaging. 

During  design,  the  review  and  analysis  of  possible  flutter  behavior 
constitutes  one  of  the  roost  highly  ccmplex  studies.  The  mass  and  stiffness 
distribution  must  be  arranged  to  prevent  flutter  from  occurring  during  normal 
operation.  There  is  the  irtplication  that  any  alteration  of  stiffness  or  mass 
distribution  due  to  service  operation  could  cause  a  possible  dangerous 
reduction  of  the  speed  at  vAiich  flutter  would  occur. 

Hiis  fact  is  irrportant  with  respect  to  all  of  the  conditions  requiring 
adequate  stiffness  of  the  primary  structure.  Any  alteration  of  stiffiiess  nay 
produce  a  dangero\as  change  in  dynamic  reponse,  vibration,  or  aeroelastic 
behavior. 

12.4.3  Service  Life  Considerations 

When  considering  the  service  life  of  a  structure,  the  entire  gamut  of 
loads  must  be  taken  into  account.  To  achieve  satisfactory  performance  during 
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service  operation,  a  structure  must  withstand  the  cumulative  effect  of  all 
varieties  of  load  that  are  typical  of  normal  use.  Normal  periods  of  overhaul, 
inspection,  and  maintenance  must  be  anticipated. 

Creep  is  the  continued  plastic  straining  of  a  part  subjected  to  stress. 
Of  course,  creep  is  of  a  particularly  serious  nature  vdien  the  part  is  subject 
to  stress  at  high  temperature  since  elevated  temperatures  reduce  the 
resistance  to  plastic  fiav.  Gas  turbine  cotponents,  reentry  configurations, 
rocket  conbustion  chambers  and  nozzles  represent  sate  of  the  typical 
structures  in  which  creep  is  important. 

If  a  part  is  exposed  to  stress  at  high  terrperature,  the  part  will 
continue  to  strain  at  a  constant  stress.  If  the  exposure  time  is  increased  to 
sane  critical  point,  the  creep  rate  will  suddenly  increase  and  failure  will 
occur.  Such  failures  due  to  creep  can  take  place  at  well  below  the  static 
iiltimate  strength  of  the  material.  Of  course,  the  creep  stress  must  be 
supplied  for  sufficient  time  to  generate  the  condition  of  failure.  It  is 
topical  of  all  metals  that  any  increase  in  applied  stress  or  tenperature  will 
increase  the  creep  rate  and  reduce  the  time  required  to  cause  failure. 

Ihe  creep  damage  is  accumulated  throu^out  the  life  of  a  structure,  with 
the  times  at  high  stress  and  tenperature  causing  the  most  rapid  rate  of 
accumulation.  In  order  for  a  structure  to  perform  satisfactorily  in  service, 
the  spectrum  of  varying  loads  and  temperatures  must  not  cause  a  critical 
accumulation  of  creep  damage.  In  other  words,  service  use  should  not  cause 
either  creep  deformation  sufficient  to  prevent  operation  or  creep  failure  by 
fracture  or  buckling.  In  sane  applications  of  turbines,  machinery  and 
mechanism,  the  limit  of  creep  deformation  may  be  the  appropriate  design 
consideration.  The  primary  airframe  structure  may  not  be  adversely  affected 
by  such  creep  deformations  and  the  final  failure  by  fractiare  or  buckling  may 
be  the  critical  consideration. 

In  special  high  tenperature  structures  the  anticipated  service  life  has 
distinct  limitations.  For  example,  gas  turbine  powerplants  may  have  turbine 
structtaral  elements  vdiich  must  be  replaced  at  regular  intervals  vhile  the 
shaft,  case,  and  compressor  withstand  only  ordinary  inspection.  If  certain 
turbine  elanents  were  required  to  demonstrate  the  same  time  life  as  less 
highly  stressed  or  heated  parts,  such  life  may  not  be  at  all  possible.  As  a 
result  the  design  seir^ice  life  of  such  high  tenperature,  hi^  stress  parts  may 
be  set  by  design  limitations  rather  than  arbitrary  desirable  values. 
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fatigue  is  the  result  of  repeated  or  cyclic  loads.  If  a  metal  part  is 
subjected  to  cyclic  stress  of  sufficient  magnitude,  a  crack  will  eventually 
form  and  propagate  into  the  cross-section.  When  the  remaining  cross-section 
cannot  withstand  the  existing  loads,  a  final  rupture  occurs  as  if  by  static 
load.  The  most  inportant  aspect  of  fatigue  is  that  the  failure  is  progressive 
by  the  accumulation  of  fatigue  damage.  When  a  critical  level  of  damage  is 
accumulated,  a  crack  foms  and  propagates  until  final  failure  takes  place. 
VJhen  a  part  is  exposed  to  the  variety  of  repeated  loads  during  service 
operation,  the  cumulative  fatigue  must  be  limited  so  that  failure  does  not 
occur  within  the  anticipated  service  life. 

In  order  to  prevent  fatigue  failures,  the  structural  design  must  bring 
into  consideration  many  inportant  factors.  First,  a  reasonable  estimate  of 
the  service  life  must  be  made  and  the  typical  spectrum  of  service  loads  must 
be  defined.  Then  the  fatigue  characteristics  of  the  materials  must  be 
determined  by  laboratory  test  of  specimens.  The  effect  of  stress 
concentrations,  corrosion  environment,  residual  stresses,  and  manufacturing 
quality  control  must  be  analyzed.  With  these  factors  known,  the  concept  of 
cumulative  damage  can  be  applied  to  determine  the  dimensions  of  the  part 
necessary  to  prevent  fatigue  failiare  during  the  anticipated  service  life. 

The  normal  scatter  and  variation  in  material  properties  encountered  in 
fatigue  tests  will  not  allow  prediction  of  the  specific  life  of  individual 
parts .  A  more  appropriate  consideration  would  be  to  account  for  the 
variabili'ly^  of  material  characteristics  and  load  spectrum  by  the  definition  of 
failure  prcbability.  In  other  words,  as  parts  in  service  approach  the  design 
service  life,  the  probability  of  fatigue  failure  increases.  If  parts  are 
exposed  to  service  well  beyond  the  design  service  life,  failure  probability 
will  be  quite  high  and  the  incidence  of  failures  and  malfunction  of  equipment 
will  increase. 

The  use  of  pieriodic  in^)ection  and  maintenance  is  very  unnecessary  to 
ensure  failure  free  operation.  Regular  inspection  must  guarantee  that  parts 
do  not  incur  excessive  deformation  or  cracks  during  exposure  to  fatigue  and 
creep  conditions.  There  is  always  the  possibility  of  short  periods  of  high 
stress  or  teirperature  vhich  could  cause  acceleration  of  creep  or  fatigue 
damage  and  precipitate  a  premature  failure.  This  is  a  very  inportant 
obligation  of  the  maintenance  facility. 
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Ihe  various  considerations  of  static  strength,  stiffness,  and  service 
life  vd.ll  all  contribute  specific  demands  on  a  structvire.  Just  vMch  of  these 
considerations  predominate  v^ill  depend  upon  the  exact  nature  of  the  structure. 
An  aircraft  shew  that  any  one  of  the  static  strength,  service  life,  or 
stiffness  requironents  may  predominate.  In  the  design  of  most  pcvi^rplant 
systems,  the  service  life  considerations  of  creep  and  fatigue  visually 
predominate.  Very  short  life  missile  structures  may  shew  that  the  static 
strength  considerations  prevail  during  design.  However,  if  the  missile  must 
vd-thstand  considerable  transportation,  handling,  and  continuous  functional 
checks,  service  life  considerations  may  be  inportant. 

12.4.4  load  and  Stress  Distribution 

Fundamentally,  there  are  two  types  of  loads  which  can  be  applied  to  an 
elenent  of  structure.  Ihere  are  "axial  loads"  vhich  are  applied  along  the 
axis  of  the  part  (Figure  12.13)  and  there  are  "transverse"  loads  -  often 
referred  to  as  "shear  loads"  -  vhich  are  applied  normal  or  perpendicular  to 
the  axis  of  the  part  (Figure  12.14).  Any  conplicated  load  condition  to  vhich 
a  structure  is  subjected  can  be  resolved  to  the  various  axial  or  transverse 
loads  acting  on  a  particular  part.  Since  the  inherent  strength  properties  of 
a  structural  material  are  based  xipon  the  element  strengths  of  its  crystals  and 
grains,  a  rcre  appropriate  definition  of  load  condition  on  a  part  is  the 
airo'unt  of  load  per  unit  of  cross  section  area.  Thxos  load  per  unit  area  is 
referred  to  as  "stress"  and  all  basic  properties  of  structural  materials  are 
based  upon  stress  or  load  per  unit  of  cross  section  area.  Of  course,  as  there 
are  two  basic  types  of  loads  —  axial  and  transverse  —  there  are  two  basic 
types  of  stress  which  result  frem  these  loads  —  axial  or  "normal"  stresses  of 
tension  and  coipression  and  the  transverse  or  "shear"  stresses. 


TENSION  COMPRESSION 


FIGURE  12.13.  AXIAL  LOADS 
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Axial  tension  applied  to  a  material  will  produce  certain  types  of  effects 
and  certain  types  of  failures,  \^le  axial  cotpression  in  the  same  type  of 
material  will  prod\Jce  ccxnpletely  different  effects  in  different  types  of 
failvires. 
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FIGUEE  12.14.  TEANSVERSE  LOADING 

A  shear  stress  applied  to  a  particular  material  results  in  a  scmev^iat 
unxisual  pattern.  If  a  shear  is  ajplied  to  an  element  of  material  in  a 
vertical  direction,  that  elonent  will  experience  balancing  shears  in  a 
horizontal  direction  of  an  equal  magnitude.  This  must  be  so  in  order  that  the 
element  be  in  equilibri-um  in  rotation. 

In  any  case  in  vdiich  the  basic  loads  —  either  axial  or  transverse  —  are 
resolved  on  a  cross-section  of  a  part,  the  stresses  are  then  ccrrputed  for  each 
element  of  the  part  as  the  amount  of  load  per  unit  of  area. 


Itor  nonrial  stress.  For  shear  stress, 

0  =  P/A  Og  =  P/A 

a  =  normal  stress,  psi  o  =  shear  stress  (12.11) 

P  =  loads  in  lbs  P  =  loads  in  lbs 

.  2  2 
A  =  area,  in  A  =  area,  in 

As  an  exaitple  of  the  typical  stress  distribution  in  a  loaded  structure, 
the  follcwing  exairple  problem  is  provided.  Ihis  exairple  will  best  furnish  an 
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interpretation  of  the  idea  of  stress  and  the  function  of  certain  structural 
coiponents. 

Figure  12.15  illustrates  a  typical  beam  structure.  In  this  case  there  is 
a  siitplified  spar  or  beam  subjected  to  a  concentrated  load  at  the  right  end. 
The  left  end  is  mounted  or  fixed  to  a  rigid  wall  surface.  For  sittplicity,  it 
is  assumed  that  the  spar  (or  beam)  is  the  entire  effective  structure  and  that 
all  loads  will  be  resisted  by  this  part.  Ihe  problem  will  be  to  investigate 
the  stresses  at  Points  A  and  B  in  this  spar  structure  vdiich  result  from  the 
application  of  the  shear  load  at  Point  C.  Ihe  beam  has  a  constant 
cross-section  throughout  the  span  as  shown  in  Figure  12.15.  The  spar  flanges 
of  this  type  structure  furnish  the  primary  bending  resistance,  vAiile  the  spar 
vreb  connecting  the  spar  flanges  provides  the  primary  resistance  to  shear 
loads.  The  rivets  attach  the  spar  web  to  the  spar  flange,  and  there  are 
vertical  stiffeners  attached  to  the  web  to  maintain  the  form,  shape,  and 
stabilily  of  the  structure.  The  distribution  of  stress  at  Point  A  in  the  beam 
is  best  visualized  by  taking  a  section  through  the  beam  at  Point  A  and 
supplying  the  internal  loads  at  the  Section  A  vhich  are  necessary  to  resist 
the  applied  external  shear  load  at  C,  thus  maintaining  ecpiilibrium  of  the 
structure.  This  is  shown  in  Figure  12.16. 
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WEB  TO  FLANGE  RIVETS  ARE  SPACED  %  IN. 


FORCE 
OF  10,000  LB 


10  IN.  BETWEEN  CENTROIDS 
OF  FLANGE  AREAS 


WEB  TO  FLANGE  RIVETS 
ARE  %  IN.  DIAMETER 


UPPER  SPAR  FLANGE 
X-SECTION  AREA  =  4.0  IN.^ 


SPAR  WEB 
IS  0.1  IN.  THICK 


LOWER  SPAR  FLANGE 
X-SECTION  AREA  =  2.0  IN.2 


FIGURE  12.15.  TYPICAL  BEAM  STRUCTURE 
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First,  to  maintain  an  equilibrium  in  a  vertical  direction,  there  must  be 
a  resisting  shear  load  on  Section  A  of  10,000  lbs  in  a  vertical  direction 
down,  vdiich  resists  the  applied  shear  load  up  at  Point  C.  Due  to  the  lever 
arm  of  force  at  Point  C,  there  is  a  bending  moment  produced  in  the  structure 
at  Point  A.  The  magnitude  of  this  bending  moment  is  the  product  of  the  force 
and  the  lever  arm,  i.e.,  10,000  lbs  times  100  in  =  1,000,000  in  lbs  of  moment 
at  Section  A.  As  the  spar  flanges  provide  the  primary  resistance  to  bending, 
there  vrLll  be  a  compression  axial  load  developed  in  the  upper  spar  flange  and 
a  tension  axial  load  produced  in  a  lower  spar  flange.  Ihese  axial  loads  in 
the  spar  flanges  for  this  untapered  beam  must  be  equal  to  maintain  equilibrium 
of  a  structure  in  a  horizontal  direction.  These  axial  forces  in  the  ppa-r 
flanges  will  be  referred  to  as  "P"  pounds  of  load.  These  two  forces  of  P, 
acting  as  a  couple  at  a  distance  of  10  in  apart  (the  distance  between  the 
centers  of  gravity  of  areas  of  the  ujper  and  lower  spar  flanges)  must  provide 
internal  balance  in  the  structure  to  the  ^plied  external  bending  moment  of 
1,000,000  in-lbs  In  other  words,  P  x  10  in  must  equal  1,000,000  in-lbs  The  P 
pounds  of  load  in  a  flange  is  then  ccmputed  to  be  100,000  lbs. 

Plb  x  10"  =  1,000,000  in-lbs 
P  =  100,000  lbs 

To  determine  the  stress  in  the  upper  flange,  the  compression  load  is 
distributed  over  the  compression  area 

=  P/A  =  100,000  lbs/4  sq  in  =  25,000  Ibs/in^ 

The  tension  stress  in  a  lonrer  flange  is  the  load  divided  by  the  area 

=  P/A  =  100,000  lb/2  sq  in  =  50,000  Ibs/in^ 

Since  the  spar  w^  furnished  primary  resistance  to  shear  loads,  there  is 
a  shear  load  of  10,000  lbs  acting  on  the  effective  area  of  the  web.  The 
effective  area  of  this  vreb  is  the  depth  times  the  thickness.  In  this  case, 
depth  is  taken  as  10  inches  and  the  thickness  is  1/10  of  an  inch  which  produce 


-  100  IN. 

UPPER  FLANGE 


FIGURE  12.16.  STATIC  STRUCTURAL  EQUILIBRIUM 

one  square  inch  of  cross  sectional  area.  iMs  shear  stress  in  the  web  is  then 
load  divided  by  area 

Og  =  P/A  =  10,000  Ib/1  sq.  in.  =  10,000  Ibs/in^ 

To  investigate  the  stress  distribution  at  Section  B,  the  same  fundaitental 
procedure  is  aTtployed  (Figure  12.17).  That  is,  the  section  at  Point  B  is 
furnished  with  the  loads  on  the  cross-section  necessary  to  place  the  structure 
in  equilibrium.  The  bending  momant  to  be  resisted  by  the  spar  flanges  is  the 
10,000  lbs  of  force  acting  at  the  50  inch  lever  arm.  This  produces  a  bending 
of  500,000  in  lbs  and  results  in  axial  loads  in  the  spar  flanges  of 
50,000  lbs  each. 


FIGCIRE  12.17.  SHEAR  STRESS  DISTRIBUTION 

The  ccnpression  stress  in  the  upper  flange  is  then 

o  =  P/A  =  50,000  Ih/  sq  in  =  12,500  Ihs/in^ 
c 

The  tension  stress  in  the  Icv^r  flange  is 

0^  =  P/A  =  50,000  lb/2  sq  in  =  25,000  Ihs/in^ 

Since  the  same  amount  of  shear  load  must  be  supplied  on  Section  B  to 
provide  equilibriim  in  a  vertical  direction  and  resistance  to  the  applied 
10,000  lbs  shear  load,  the  shear  stress  in  the  spar  v?eb  remains  the  same 
10,000  psi  throu^out  the  span  of  the  beam  and  it  does  so  as  long  as  there  is 
no  change  in  the  shear  load  across  the  beam. 

To  determine  the  stress  in  the  rivets  attaching  the  spar  web  to  the 
flanges,  the  function  of  the  rivets  must  be  made  clear.  The  point  of  load 
application  at  Section  C  on  the  beam  has  a  vertical  load  of  10,000  lbs 
applied.  This  concentrated  load  of  10,000  lbs  must  be  appropriately 
distributed  to  the  vieib  by  a  fitting.  The  desired  result,  as  in  Figure  12.18, 
is  to  distribute  the  concentrated  shear  load  of  10,000  lbs  to  the  edge  of  the 
spar  such  that  for  the  10  inch  effective  depth  of  the  web  there  is  1,000 
lbs  of  load  for  each  1  inch  along  the  edge. 
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PURPOSE  OF  THE  FITTING  IS  TO 
DISTRIBUTE  THE  CONCENTRATED  LOAD 
OF  10,000  LB  TO  A  UNIFORM  SHEAR 
IN  THE  WEB  OR  1 000  LBS/IN. 


FIGURE  12.18.  DISTRIBUTIO^T  OF  CCNCENTRATED  SHEAR  LOAD 

Ihe  ElCTient  1  in  the  upper  right  hand  edge  of  the  web  at  Section  C  has 
applied  on  its  edge  1,000  lbs  shear  load.  Figure  12.19  illustrates  the  manner 
in  vhich  this  1,000  lbs  of  load  applied  to  this  1  inch  Element  1  is  resisted. 

At  tile  lefr  hand  side,  there  is  a  shear  load  dov-n  of  1,000  Lbs  balancing 
the  applied  1,000  lbs.  This  shear  load  on  the  left  hand  edge  is  furnished  by 
the  adjacent  elenent  of  the  web  to  the  left  of  Element  1.  Since  an  element  of 
structure  with  an  applied  shear  load  must  also  have  balancing  shears  at  90°, 
there  v/ill  exist  (or  must  exist)  on  Element  1  shear  loads  on  the  upper  and 
lower  edges  as  shown .  The  shear  load  on  the  upper  edge  of  Element  A  is 
supplied  by  the  next  piece  of  structure  in  contact  with  the  edge  of  the  web. 
This  load  must  cone  fran  the  spar  flange  and  is  transmitted  by  the  "web-to- 
flange"  rivets.  The  primary  function  of  the  "web-to-flange"  rivets  is  then  to 
provide  a  continuity  of  shear  and  to  balance  the  applied  vertical  shear  load 
in  a  horizontal  direction.  The  shear  load  of  the  lower  side  of  Element  1  is 
supplied  by  the  next  element  imnediately  underneath  Elonent  1  in  the  web.  If 
the  spacing  of  rivets  attaching  the  web  to  the  flange  is  three-quarters  of  an 
inch,  the  load  for  each  rivet  in  shear  will  be  3/4  of  1,000  lbs  (rivets  spaced 
3/4  inch)  or  750  lbs  per  rivet.  With  the  diameter  of  the  rivet  given  as  3/8 
of  an  inch,  the  rivet  stress  in  shear  could  then  be  ccnputed  as  the 
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load  divided  by  area 


Rivet  stress 


load 

area 


Og  =  750  lbs 

{Tr/4)  (3/8)^in.2 

750  lbs 

s  - r“ 

.1104  in^ 
o  =  6800  psi 

W 


1000  LB  OF 
DISTRIBUTED 
SHEAR  LOAD 
(1000  LBS/IN.) 


FIGURE  12.19.  SHEAR  ELEMENT 
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FIGURE  12.20.  SPAR  FLANGE  AND  WEB  LOADS 

Figure  12.20  gix’es  axi  ilii;istration  of  the  upper  and  lower  spar  flange 
removed  frcm  the  structure  and  the  loads  applied  to  spar  flanges  and  the  web. 

The  flange  loads  of  100,000  lbs  at  Section  A  are  the  result  of  the 
accumulated  axial  force  from  the  distributed  shear  loads  of  1,000  lbs  per  inch 
for  the  entire  length  of  100  inches.  This  distributed  load  at  Section  A  (100 
inch  length)  and  50,000  lbs  at  Section  B  (50  inch  length). 

The  vertical  stiffeners  attached  to  the  web  have  no  particular  stress  — 
either  catpression  tension  or  shear  -  until  the  web  of  the  spar  begins  to 
buckle .  (Figure  12.21) 

The  primary  function  of  these  vertical  stiffeners  is  to  maintain  the  form 
and  stability  and  to  provide  support  for  the  web  panel,  thereby  preventing  or 
delaying  the  shear  buckling  of  the  web.  When  buckling  occurs  in  the  web,  the 
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FLANGE 


FIGURE  12.21.  VERTICAL  STIFFENERS 

vertical  stiffeners  must  withstand  conpression  loads  to  prevent  collapse  of  the 
structure. 

The  function  and  iitrooirtance  of  the  spar  web  can  be  best  oiphasized  by  two 
exaitples  shown  in  Figure  12.22.  If  there  were  no  web  between  the  pinned 
flanges  and  a  shear  load  applied,  no  resisting  loads  would  be  developed  in  the 
flanges.  The  structure  woiild  collapse  to  the  shape  shown  with  no  resistance 
to  the  applied  load. 


12.56 


STRUCTURE  OFFERS 
NO  RESISTANCE  TO  AN 


t 


BENDING  OF  FLANGES 
PRODUCES  LITTLE 
RESISTANCE  TO  AN 
APPLIED  LOAD 


THE  WEB  RESISTS  SHEAR 
AND  ALLOWS  BENDING  TO  BE 
RESISTED  BY  AXIAL  LOADS 
IN  THE  FLANGES 


FIGURE  12.22.  SPAR  VEB  FUNCTION 


jf  V7sir0  no  wsb  bstwson  ths  fixsd  flsngss,  3  shssir  load,  applisd  would 
produce  "secondary"  bending  of  the  flanges.  Since  the  bending  resistance  of 
the  flanges  is  quite  smalls  prohibitive  stresses  and  deflections  would  result. 
Obviously,  the  more  efficient  structure  is  the  web-flange  combination  vhich 
resists  shear  in  the  vreb  and  bending  by  axial  loads  in  the  flanges. 

Actually,  there  is  very  rarely  such  a  thing  as  a  "lightening  hole"  in  a 
shear  web.  More  usually  these  holes  are  for  access  in  maintenance  and 
production  and  are  a  structural  penalty  for  anything  but  minimum  gage 
thickness  structures  (li^t  planes,  gliders,  airships,  etc.). 

The  previous  example  problem  of  the  stress  distribution  in  a  simplified 
structure  has  considered  that  there  were  no  particular  corplications  to 
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praiuce  anything  other  than  pure  axial  and  shear  stresses.  In  a  more  detailed 
analysis,  consideration  vroiiLd  be  made  of  the  contribution  of  the  web  to 
bending  resistance,  the  ccnplication  and  magnification  of  stress  due  to  rivet 
holes,  etc. 

Ihere  are,  of  course,  exaitples  in  which  the  distribution  of  stress  in  a 
structural  element  is  complicated  by  the  particular  manner  of  loading.  There 
follow  certain  examples  of  the  particular  stress  as  distributed  in  typical 
structures  if  ttere  is  bending,  torsion,  etc. ,  and  the  resulting  stresses 
remain  in  the  elastic  range  of  the  material  characteristics.  "elastic 
range",  it  is  implied  tiiat  stresses  may  be  applied,  then  released,  and  no 
pennanent  deformation  of  the  structure  would  be  incurred. 
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Figure  12.23  illustrates  the  use  of  a  solid,  rectangular  bar  with  pure 
bending  nicments  applied.  The  stresses  at  Section  A  will  be  distributed  as 
shown  with  the  upper  portion  of  Section  A  subjected  to  a  ccmpression  stress 
v^ch  will  be  a  maximum  at  the  outer  surface  and  the  lower  s\irface  of  Section 


A  sxibjected  to  a  tensile  stress  v^ch  will  be  a  maximum  at  the  outer  surface. 

The  point  at  vdiich  the  stress  is  zero  -  neither  tension  nor  ccmpression  - 
is  referred  to  as  the  "neutral  axis".  For  a  symmetrical,  rectangular  section 
this  point  would  be  midway  between  the  upper  and  Icwer  surfaces.  The  neutral 
axis  for  a  hcanogenous  inaterial  subjected  to  elastic  bending  is  always  located 
at  the  center  of  gravity  (or  "centroid")  of  cross  sectional  area. 


NEUTRAL 

AXIS 


FIGURE  12.23.  PURE  BENDING  OF  A  SOLID,  RECTANGULAR  BAR 

In  Figure  12.24  there  is  an  vinsymitietrical  cross-section  subjected  to  pure 
bending.  For  the  cross-section  shown,  the  neutral  axis  will  be  closer  to  the 
upper  surface  of  the  part  than  the  lc^^7er  surface.  The  stress  distribution 
illustrated  in  Figure  12.24  will  continue  to  be  a  linear  variation  of  stress 
between  the  two  maximum  stresses  at  the  upper  and  lower  surfaces.  However,  as 
the  neutral  axis  is  closer  to  the  upper  surface,  the  magnitude  of  carpression 
stress  will  be  smaller  than  the  tensile  stress  on  the  lower  surface.  This  must 
be  so,  since  the  ccxtpression  load  produced  by  the  smaller  cotpression  stress 
distributed  over  the  larger  area  above  the  neutral  axis  will  be  equal  to  the 
tension  load  produced  by  the  higher  tensile  stress  acting  over  the  smaller 
tensile  area  below  the  neutral  axis.  Thus,  equcil  and  opposite  compression  and 
tensile  loads  exist  which  furnish  equilibrium  to  the  cross-section  in  the 
horizontal  direction. 


12.59 


(b 


SECTION  A 


T* 


NEUTRAL 

AXIS 


TENSION 


FIGURE  12.24.  PURE  BENDING  CF  AN  UNSYMMBTRICAL  CROSS-SECTION 
Par  beams  or  other  structures  subjected  to  bending  ironents, 


a 


b 


max 


to 

I 


(12.12) 
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=  maxiirtum  bending  stress,  psi 


Where  a, 
b 

max 

M  =  bending  moment  on  cross-secrtion,  lb  in 
c  =  distance  from  neutral  axis  to  outermost  fiber,  in 
I  =  mcment  of  inertia  of  cross-section,  in^ 


If  the  stress  on  some  other  fiber  is  desired,  the  equation  is  simply  = 

f^/I,  vtere  y  is  the  distance  from  the  neutral  axis  to  the  fiber  under 
consideration. 

The  term  "mcment  of  inertia"  appears  in  the  previous  equations.  As  has 
been  explained,  the  bending  stress  in  a  fiber  depends  on  its  distance  from  the 
neutral  axis.  If  the  material  can  tolerate  a  given  amount  of  stress,  the 
largest  mcment  can  be  resisted  vAien  the  area  resisting  it  is  as  far  as 
possible  frcm  the  neutral  axis,  itie  moment  of  inertia  is  a  quantity  which 
takes  the  shape  of  the  cross-section  into  account  in  determining  the  amount  of 
bending  mcment  which  can  be  resisted  by  a  given  cross-section  without  exceed¬ 
ing  a  specified  stress.  For  several  cross-section  shapes  the  following  apply; 
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mcment  of  inertia  about  the  neutral  axis,  in'^ 
rectangle  or  triangle  width,  in 
rectangle  or  triangle  hei^t,  in 
outside  diameter,  in 
inside  diameter,  in 
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12.4.6  Pure  Torsion 

Figure  12.25  illustrates  the  condition  of  pure  torsion  applied  to  a 
solid,  circular  shaft.  With  torsion  applied  to  the  part,  the  stress  produced 
on  Section  A  is  primarily  that  of  shear  -  a  shear  stress  which  is  a  maximum  at 
the  outer  surface  and  varies  linearly  to  zero  at  the  axis  of  the  part. 

For  circular  shafts  subjected  to  torsion  matients. 


a 

s 

max 


J 


(12.13) 


Where  a  =  maximum  shearing  stress,  psi 
s 

max 

T  =  torsional  monent  or  twisting  monent,  lb  in 

c  =  radius  frcm  axis  to  outermost  surface,  in 

.  4 

J  =  "polar  moment  of  inertia"  of  cross-section,  in 


Use  f  =  Mr/J  if  the  stress  at  a  distance  r  from  the  axis  is  desired. 

The  polar  moment  of  inertia,  J,  is  used  for  determining  strength  of 
sections  subjected  to  torsion. 

Only  the  circular  sections  have  simple  equations  for  this  sort  of  stress 
distribution.  These  are  as  follovre: 


J  =  ^2^0^  (circular) 

~  §2  ^0^  “  ^i^  (doughnut) 

.  4 

Inhere  J  =  polar  moment  of  inertia,  in 
Dq,  =  outside  and  inside  diameters,  in 

For  sections  other  than  circular,  reference  may  be  made  to  any  of  the 
more  standard  texts  on  strength  of  materials. 
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FIGURE  12.25.  PURE  TORSIOSf  OF  A  SOLID,  CIRCULAR  SHAFT 

Figure  12.26  illustrates  the  conditions  in  vMch  a  continuous,  hollow 
cross  section  is  subjected  to  torsion.  In  this  instance,  the  shear  stress  is 
a  constant  value  around  the  periphery  of  the  cross  section.  Since  the  shear 
load  distributed  around  the  periphery  must  provide  an  internal  resisting 
moment  equal  to  the  external  applied  moment,  the  shear  stress  for  this  contin¬ 
uous  shell  structxure  may  be  computed  by  use  of  the  following  relationship: 
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o_  = 


■2aE 


v^iere  c  = 
s 

T  = 
A  = 
t  = 


shear  stress,  psi 

applied  torque,  in-lbs 

enclosed  area  of  cross-section,  sq  in 

shell  thickness,  in 
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FIGURE  12.26.  PURE  TORSION  OF  A  HOLLOW  TUBE 
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If  this  section  were  slotted  in  a  longitudinal  direction,  a  very  large 
amount  of  the  torsional  rigidity  would  be  lost.  There  would  be  no  continuous 
shear  stress  distributed  around  the  periphery  of  the  shell  and  high  local 
shear  stresses  with  great  deflections  would  be  encountered.  The  maximum  shear 
stress  in  this  case  could  be  calculated  by  application  of  the  following 
equation: 


a 


max 


where  a 

s 

max 


maximum  shear  stress,  psi 


(12.15) 


T  =  applied  torque,  in-lbs 

b  =  unrolled  width,  in 

c  =  shell  thickness,  in 

Any  shell  structxure  subjected  to  torsional  loading  vhich  has  a  cutout  or 
slot  will  have  a  tendency  to  develop  much  higher  stresses  and  m^  be 
excessively  flexible. 


12.4.7  Bolted  or  Riveted  Joints 

Figure  12.27  illustrates  a  typical  type  of  bolted  or  riveted  joint  v^iich 
is  encountered  in  more  conventional  structures.  In  such  a  joint,  the  load 
applied  to  element  A  is  transferred  by  the  bolt  and  distributed  to  elements  B 
and  C.  This  transfer  of  load  by  the  bolt  produces  a  shear  stress  on  the  bolt 
cross-section.  The  bolt  shear  stress  is 


a 

s 


load 

area 


P/2  _  _ ^ 

ir/ 4  d^  ird^ 


(12.16) 


The  tensile  load  developed  in  the  plate  creates  a  critical  tensile  stress 
along  section  x.  The  average  tension  stress  at  section  x  is 


load 

area 
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P 

(w-d)t 


(12.17) 


12.65 


FIGUPE  12.27.  DOUBLE  LAP  JOINT,  BOLT  IN  DOUBLE  SHEAR 

In  addition,  the  bolt  bearing  on  the  s\rrface  of  the  hole  (surface  Y)  may 
create  a  critical  oarpressive  stress.  The  average  bearing  stress  in  the  plate 
at  the  bolt  hole  is: 

“hr  =  ^  <12.18) 

“hr  '  W 


The  application  of  the  bearing  load  creates  a  shear  stress  along  Sections 
z  vdiich  tends  to  tear  out  the  edge  of  the  plate. 


Lilij 


The  "edge  tear-out"  stress  is 

a  = 

s  area 

°s  =  2^  (12.19) 

These  sirrple  equations  represent  only  the  siirple  average  stresses  in 
order  to  appreciate  seme  of  the  fundamental  requirements  of  a  joint.  Ihere  is 
the  obvious  possibility  that  friction  between  the  plates  may  acccirplish  part 
of  the  shear  transfer  and  the  hole  may  create  considerable  stress  concentra- 
cause  peak  stresses  vrell  above  the  cenputed  average. 

12.4.8  Pressure  Vessels 

Hie  use  of  hi^ly  pressurized  containers  in  various  aircraft  and  missiles 
creates  significant  structural  problems.  Seme  of  the  more  siirplified 
situations  are  represented  by  the  stresses  created  in  pressurized  spherical 
and  cylindrical  shells.  These  vessels  are  illustrated  in  Figure  12.28. 


t 


FIGURE  12.28.  STRESS  IN  PRESSURIZED  VESSELS 
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The  stress  in  a  pressurized  spherical  shell  is  vmifom  and  constant  in 
all  directions  along  the  surface.  The  resulting  tensile  stress  is  due  to  the 
pressure  load  being  distributed  over  the  effective  shell  area 


load 
-  area 


(12.20) 


load  =  (pressure)  (area) 
=  (P)  (irr^) 

area  =  2trrt 
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2 

pirr 

2irrt 
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In  order  to  consider  the  stresses  in  a  pressurized  cylindrical  shell,  the 
existence  of  two  separate  stresses  must  be  noted.  The  longitudinal  stress,  o 
is  related  as  follows: 


_  load 
X  area 
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(12.21) 


Ibis  is  identical  to  the  relationship  developed  for  the  pressurized 
spherical  shell.  However,  an  additional  stress,  o^,  is  developed  vMch  is 
referred  to  as  the  "hoop"  stress. 
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load  per  unit  length  =  2pr 


area  per  unit  length 
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Oy  =  ^  (12.22) 

(NOTE:  Gy  =  2f^) 

For  this  situation  the  hoop  stress  incurred  is  twice  the  longitudinal 
stress  developed.  inportant  fact  is  concluded  frcm  this  relationship:  If 
a  failure  due  to  pressure  occurs  in  a  \iniform  cylindrical  shell,  the  hoop 
stresses  will  predominate  in  the  mode  of  failure. 

12.4.9  Catppnent  and  Principal  Stresses 

Twd  factors  determine  the  strength  and  manner  of  failure  of  any 
structural  mentoer.  Cne,  of  course,  is  the  property  and  character  of  the 
structural  material.  Hie  other  is  the  maximum  normal  stresses  ("principal 
stresses")  and  maximum  shear  stresses  vMch  exist  in  various  areas  of  a  part. 
T^enever  there  is  a  normal  stress  applied  to  a  part  there  will  exist  various 
magnitudes  of  normal  and  shear  stresses  at  planes  different  to  the  direction 
of  loading.  Also,  any  time  a  shear  stress  is  applied  to  a  part,  there  will 
exist  at  certain  planes  various  magnitudes  of  normal  and  shear  stresses.  Hie 
various  ccmponents  of  the  applied  primary  stress  must  be  investigated  to 
determine  the  influence  upon  strength  and  failure  type. 

A  typical  example  of  catponent  stress  is  illustrated  in  Figure  12.29. 
Figure  12.29A  shows  a  specimen  of  material  with  a  pure  axial  tension  load 
applied.  If  a  section  or  cut  is  made  in  this  specimen  at  Section  X,  it  is 
seen  that  a  constant  tension  stress  exists  on  this  plane  perpendicular  to  the 
direction  of  loading.  Figure  12.29B  shows  the  same  specimen  subjected  to  the 
same  pure  axial  tension  load  but  with  a  cut  made  along  Section  Y.  An 
investigation  of  the  stresses  acting  on  Section  Y  shows  that  there  still 
exists  a  uniform  tensile  stress  across  the  section.  However,  the  tensile 
stress  vAiich  exists  on  plane  Y  will  have  caiponents  vMch  are  perpendicular 
and  parallel  to  the  surface  of  the  cut.  One  of  these  conpanents  which  is 
perpendicular  to  the  face  of  the  section  will  be  a  normal  tensile  stress  of  a 
smaller  magnitude  than  the  tensile  stress  applied  to  the  cross-section.  Along 
the  surface  of  Section  Y  will  be  a  corponent  of  stress  parallel  to  the 
surface.  Hus  is  a  shear  stress  —  it  is  a  conponent  of  the  primary  ^plied 
tensile  stress. 
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since  the  tensile  stress  conponent  perpendicular  to  the  face  of  the  cut 
will  be  of  a  sinaller  magnitude  than  the  applied  tensile  stress,  it  will  hcive 
no  particular  or  imnediate  influence  on  strength  or  irode  of  failure.  Ifawever, 
the  shear  stress  conponent  parallel  to  the  face  of  the  section  is  a  stress  of 
an  altogether  different  nature.  It  may  have  a  decided  effect  on  strength  and 
failure  type  as  seme  materials  are  much  more  critical  in  shear  than  in  normal 
stress  (particularly  ductile  itetals) . 


FIGURE  12.29.  COMPONENT  STRESSES 

Figure  12.30A  shews  a  part  loaded  in  tension  with  a  section  cut  along  the 
direction  of  applied  stress.  It  is  obvious  there  would  be  no  shear  stress 
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(due  to  force  conponents)  along  this  Section  Z.  Since  shear  stresses  do  not 
exist  either  at  a  section  perpendicular  to  the  direction  of  applied  stress,  or 
at  a  section  parallel  to  the  direction  of  applied  stress,  it  is  reasonable  to 
assvjme  that  between  these  liinits  the  existing  shear  stress  will  be  a  maxirnum 
on  a  section  at  45°  -to  the  direction  of  primary  load  application. 


PARALLEL  TO  APPLIED  NORMAL  STRESS 


T,  AXIAL  OR 
NORMAL  FORCE  ON 
ELEMENT  DUE  TO 
APPLIED  STRESS 


(b) 


S,  SHEAR 
COMPONENT 
OF  FORCE  T 


FIGURE  12.30.  CXDMPCXgENT  SHEAR  STRESSES 


To  determine  the  magnitude  of  the  maximum  shear  stress,  it  is  best  to 
take  a  one  inch  cube  of  material  frcxn  the  basic  specirten  and  examine  the 
forces  existing  in  rhis  element.  There  is  applied  to  the  sides  of  this  one" 


inch  element  (as  in  Figure  12.29B)  a  load  vMch  is  the  tensile  stress  on  the 
speciroen  cross-section.  Since  the  section  or  cut  is  to  be  located  at  45°  to 
the  applied  stress,  the  ccqponent  of  force  distributed  along  the  diagonal  vdll 
be  as  follows: 

S  =  T  sin  45° 

S  =  .707  T 

This  shear  force  will  be  distributed  over  the  diagonal  surface  which,  for 
the  one  inch  element,  will  be  1.414  inches.  The  shear  stress  is  then 

=  .707  T 

1.414 

=  .50  T 

Thus,  an  applied  axial  stress  will  produce  a  maximum  coiponent  shear  stress 
which  is  one-half  the  magnitude  of  the  applied  stress.  Knowledge  of  the 
presence  of  the  cotponent  shear  stress  and  its  existence  as  a  maximum  at  45° 
is  basic  to  a  discrimination  between  ductile  and  brittle  failure  types. 

An  exanple  of  the  existence  and  orientation  of  the  maximum  cotponent 
shear  stress  is  provided  in  Figure  12.31  vdiere  two  specimens  of  metal  -  one 
very  ductile  and  one  very  brittle  -  are  subjected  to  failing  tensile  loads. 
The  brittle  type  material  will  fracture  with  a  clean  break  at  90°  to  the 
direction  of  loading.  The  ductile  specimen  will  exhibit  fracture  planes  at 
45°  to  the  direction  of  the  applied  tensile  stress  thus  verifying  the 
existence  of  conponent  shear  stresses.  This  45°  of  fracture  fron  pure 

tensile  loads  is  referred  to  as  the  "ductile"  or  "shear"  type  of  failure. 
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FIGURE  12.31.  EFFECT  OF  DUCTILnY  W  FRACTURE 

The  condition  of  an  applied  pure  shear  presents  a  different  and  sli^tly 
nore  ccxrplex  problem  concerning  ccmponent  and  principal  stresses.  Consider  a 
one  inch  cubic  element  subjected  to  pure  shear  (Figure  12.32) . 

if  s  Section  A  is  talcen  at  a  diagonal  of  45^  in  one  direction,  it  is 
apparent  that  there  must  be  a  catpression  force  on  the  diagonal  to  statically 
balance  the  action  of  the  two  shear  forces  applied  along  the  edges  of  the 
element.  The  two  shear  forces  have  ccnponents  at  45*^  v^^ich  are  additive  and 
must  be  balanced. 


C  =  S  sin  45*^  +  S  cos  45° 

C  =  (2)  (.707)S 

C  =  1.414  S 

Since  the  ccttpression  force  is  distributed  along  the  diagonal  the  ccrpression 
stress  may  be  found  as 


carpression  stress 


1.414  S 
1.414 


=  1.000  S 
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FIGURE  12.32.  APPLIED  PURE  SHEAR  -  SECTICN  AA  ANALYSIS 

Thus,  for  the  case  of  pure  shear  applied  to  an  element,  there  ^ists  at  one 
45®  section  a  cctrpression  stress  equal  in  magnitude  to  the  applied  shear 
stress.  The  existence  of  this  conponent  ccnpression  stress  rs  evident  in  the 
failure  mode  of  a  thin  walled  tube  subjected  to  a  failing  torsion  load 
(torsion,  of  course,  produces  a  uniform  piure  shear) .  If  the  walls  of  the  tube 
are  sufficiently  thin,  the  tube  will  fail  primarily  in  buckling  due  to  the 
princip>al  corpression  stress  existing  in  the  45®  direction. 


FOLDS  OR  BUCKLES 
OCCUR  PERPENDICULAR 
TO  PRINCIPAL  COMPRESSION 


FIGURE  12.33.  BUCKLING  FAILURE  DUE  TO  TORSION 
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If  the  same  one-inch  element  of  Figure  12.32  (shown  in  Figure  12.34)  is 
subjected  to  the  same  pure  shear  condition  -  but  sectioned  along  Plane  B,  a 
different  stress  situation  results.  If  the  element  is  cut  along  Section  B,  it 
is  apparent  that  a  tension  force,  T,  must  be  sufficient  to  balance  the 
cotponents  of  the  -two  shear  forces  along  the  edge.  Ihis  situation  will 
produce  a  tension  stress  distributed  along  the  diagonal  v^ch  is  pqi-ia)  in 
magnitude  to  the  applied  shear  stress.  The  presence  of  this  tension  stress  is 
by  "ths  mode  of  failure  of  a  brittle  shaft  in  torsion  (twisting  a 
piece  of  chalk  produces  the  same  result)  (Figure  12.35).  In  this  case,  a 
fracture  will  begin  as  a  45°  spiral  surface  \^ich  is  perpendicular  to  the 
principal  tension  stress. 

It  must  be  remenfcered  that  catponent  stresses  may  have  a  definite  bearing 
on  the  strength  and  mode  of  failure  of  any  structure.  Any  normal  or  shear 

stress  applied  to  a  part  will  produce  ccnponent  and  principal  stress  that 
cannot  be  neglected. 


FIGURE  12.34.  APPLIED  PURE  SHKAR  -  SECTION  B-B  ANALYSIS 
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TORSION  APPLIED 
CLOCKWISE 
FROM  THIS  END 


FRACTURE  SURFACE  BEGINS  PERPENDICULAR 
TO  THE  PRINCIPAL  TENSION  STRESS 


FIGURE  12.35.  BKETITiE  FAILURE  DUE  TO  TORSION 


12.4.10  Strain  Resulting  fron  Stress 

Any  structure  v^iich  is  subjected  to  stress  must  deform  under  load,  even 
though  the  deformation  may  not  be  visible  to  the  naked  eye.  Recall  frcjn  the 
previous  sections  that  stress  is  the  true  measure  of  state  for  a  part 
subjected  to  load.  Strain  is  a  similar  means  of  measure.  In  order  to  fully 
evaluate  the  state  of  being  of  a  stressed  material,  all  deformations  must  be 
considered  on  a  unit  basis.  Strain  is  thus  defined  as  deformation  per  \mit  of 
length  and  in  most  engineering  terminology  is  referred  to  by  "e"  or  "e" 
(epsilon) . 

Suppose  a  steel  bar  100  inches  in  length  is  subjected  to  a  tensile  stress 
of  30,000  psi  -  as  in  Figure  12.36.  The  total  change  in  length  throughout  the 
100-inch  length  would  be  about  0.1  inches.  If  subjected  to  uniform  stressing, 
the  part  woxild  be  subjected  to  a  uniform  strain  which  is  as  follows: 


strain 


total  deformation 
original  length 


e  = 


AL 

L 


0.1  in 
100  in 


£  =  .001  inches  per  in 


(12.23) 
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To  produce  a  total  deformation  of  0.1  inches  in  100  inches  of  length,  the 
strain  must  be  .001  in/in  or  0.1%. 
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FIGURE  12.36.  AXIAL  STRAIN 

In  addition  to  the  longitudinal  strain  of  the  part  shewn  in  Figure  12.36, 
there  will  be  a  lateral  contraction  of  the  metal.  For  most  metals  there  is  a 
definite  relationship  between  this  transverse  strain  and  the  longitudinal 
strain  and  it  is  iirportant  vdien  considering  cotibined  stresses  that  are 
dependent  upon  deflections.  The  proportion  between  the  lateral  strain  and  the 
longitudinal  strain  has  been  given  the  name  "Poisson's  ratio,  y  (mu.)"  and  for 
most  hcmogenoiis  metals  this  ratio  has  the  approximate  magnitude  of  0.3  (y  = 
0.3).  In  the  case  of  the  previous  exanple  the  longitudinal  extension  strain 
of  .001  would  be  aocaipanied  by  a  lateral  contraction  strain  of  .0003. 

Thus,  it  is  uiportant  to  remeirber  that  any  iretal  subjected  to  stress  must 
strain.  The  amount  of  strain,  vAiile  not  necessarily  visible,  must  be  present 


12.77 


and  is  very  iitportant.  Just  as  large  stresses  may  produce  numerically  small 
strain,  any  small  strain  forced  on  a  structure  may  produce  large  stresses. 

Shear  stress  also  produces  shear  strains.  Shear  strain,  vAiile  not 
necessarily  denoting  a  change  in  length,  does  describe  a  change  in  relative 
position  of  the  part.  Figure  12.37  should  illustrate  this  fact. 


FIGURE  12.37.  SHEAR  SI'RAIN 

Shear  strain  is  most  properly  described  by  the  angle  of  strain,  y  (gamma)  in 
radians.  When  a  shear  stress  is  applied  to  an  element,  the  shear  strain,  y, 
can  be  ccnputed  as  the  proportion  of  the  change  in  position  cf  cne  side.  A*,, 
to  the  original  length,  £. 

T  =  4^  (12.24) 

One  of  the  most  inportant  properties  of  an  aircraft  material  is  its 
stiffness.  If  a  part  vrere  subjected  to  a  particular  level  of  stress,  small 
strains  would  indicate  a  stiff  or  rigid  material;  large  strains  would  indicate 
a  flexible  material.  The  accepted  method  of  measuring  the  stiffness  or 
rigidity  of  a  material  is  to  cotpute  a  proportion  between  the  applied  axial 
stress  and  the  resulting  axial  strain.  This  proportion  is  kncwn  as  the 
"Moduliis  of  Elasticity"  or  "Young's  Modulus"  and  is  denoted  by  the  letter,  E. 

E  =  -  (12.25) 

e 
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vtere  E  =  modulus  of  elasticity,  psi 
0  =  stress,  psi 
e  =  strain,  in/in 

Typical  values  for  tMs  Modulus  of  Elasticity  are 

Steel  E  =  30,000,000  psi 

Aluminvim  Alloy  E  =  10,000,000  psi 
Magnesium  Alloy  E  =  7,000,000  psi 

By  a  conparison  of  these  values  it  is  seen  that  an  aluminum  alloy  part 
subjected  to  a  given  stress  vTould  strain  three  times  as  greatly  as  a  steel 
part  subjected  to  the  same  stress  level. 

There  is  no  true  modulus  of  elasticity  in  shear.  However,  for  ccnputing 
shear  deflections  there  exists  an  equivalent  quantity  known  as  the  "modxilus  of 
rigidity"  or  "shear  modulus"  and  is  denoted  by  the  letter  G. 


G 


a 


s 


Y 


vdiere  G 
a 

s 

Y 


modulus  of  rigidity,  psi 
shear  stress,  psi 
shear  strain,  in/in 


(12.26) 


For  most  metals  the  shear  modultis  is  approximately  40%  of  the  elastic  modulus 
-  e.g. 


Steel  G  =  12,000,000  psi 

Aluminum  Alloy  G  =  4,000,000  psi 


For  homogenous  materials  the  Modulus  of  Rigidity  may  be  determined  by  the 
following  equation: 


2(l+y) 


(12.27) 
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vtere  G  =  Modulus  of  Rigidity 

E  =  Modulus  of  Elasticity 

y  =  Poisson's  ratio 

Bending  stresses  will  caiase  bending  deflections.  In  this  case,  no  general 
strain  relationship  can  be-  defined  that  is  similar  to  siitple  axial  strain. 
Since  bending  stresses  do  vary  throughout  the  cross-section  there  will  be  a 
variation  of  axial  strains  proportional  to  the  axial  stress. 

For  the  initially  straight  beam  that  is  shown  in  Figure  12.38,  an  applied 
pure  bending  incment  will  produce  a  deflection  of  pure  curvature. 


FIGURE  12.38.  PURE  CURVATURE 


Pure  bending  imposed  will  caxise  catpression  strains  on  the  lower  surface,  and 
zero  strain  at  the  neutral  axis.  The  result  of  these  strains  is  a  bending  of 
the  beam  to  the  arc  of  a  circle  vrith  no  change  in  length  of  the  neutral  axis. 
Of  course,  if  one  end  of  the  beam  is  held  stationary  the  other  end  deflects 
upward  -  but  only  because  of  the  curvature  of  bending  (Figure  12.39) . 


12.80 


FIGURE  12.39.  BENDING  DEFliBCTION  OF  CANTILEVER  BEAM 


A  scmewhat  similar  condition  exists  for  a  length  of  shaft  subjected  to  a 
pure  torsion  loading.  The  shear  stress  distributed  on  the  cross-section  vd.ll 
produce  shear  strains,  which  are  angular  displacements.  The  net  effect  is  to 
produce  a  laniform  tv^ist  throu^out  the  length  of  the  shaft.  If  a  straight 
line  vere  to  be  drawn  on  the  shaft,  as  in  Figure  12.40,  this  line  would  be 
displaced  upon  load  application  and  vrould  finally  occupy  the  position 
indicated  by  the  dotted  line.  The  helix  angle  of  displacorient  wotild  depend 
upon  the  shear  strains  developed  at  the  surface  of  the  shaft. 


FIGURE  12.40.  ANGULAR  DISPLACEMENT  DUE  TO  TORSICW 

Actiaally,  only  tj^o  factors  determine  the  strain  for  a  particular  material 
subjected  to  stress.  One,  of  course,  is  the  magnitude  of  stress;  the  second 
is  the  type  of  material  and  characteristic  stiffness.  The  actual  amount  of 
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deflection  of  a  loaded  structure  will  depend  on  the  j^ysical  arrangement  of 
the  structure  and  the  cumulative  effect  of  the  local  strains  existing  in 
various  coiponents  of  the  structure. 

12.4.11  Stress-Strain  Diagrams  and  Material  Properties 

In  order  to  evaluate  the  properties  of  a  material  and  the  possible 
structural  application,  it  is  necessary  to  determine  strains  corresponding  to 
various  levels  of  applied  stress.  Laboratory  tests  are  then  conducted  v^ich 
subject  a  specimen  of  material  to  various  magnitudes  of  stress  while  strains 
are  recorded  at  each  stress.  If  the  stresses  and  corresponding  strains  are 
then  plotted  on  a  graph,  many  useful  and  inportant  properties  of  the  material 
may  be  observed. 


FIGURE  12.41.  STRESS-STRAIN  DIAGRA^5  OF  TYPICAL  MUD  STEEL 

Figure  12.41  shows  a  typical  "stress-strain"  diagram  for  a  mild  steel  at 
roan  temperature.  As  the  stress  is  first  applied  and  increased,  the  strain 
begins  and  increases  in  direct  proportion  to  the  stress  -  i.e. ,  the 
stress-strain  diagram  is  a  strai^t  line  fron  A  to  B;  therefore,  this  material 
woiiLd  snap  back  to  zero  strain  upon  the  release  of  stress.  Because  of  the 
elastic  nature  of  the  material  in  this  range  of  stress,  the  stress  range  fron 
A  to  B  is  referred  to  as  the  "elastic  or  proportional  range"  of  the 
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material.  This  definition  then  iirplies  that  stresses  up  to  Point  B  will  not 
caxise  permanent  deformation  of  the  material. 

If  the  applied  stress  is  gradually  increased  above  the  value  at  Point  B, 
the  plot  of  stress  vs.  strain  will  deviate  sli^tly  from  a  straight  line. 
There  will  then  be  sane  small  -  but  measurable  -  permanent  strains  thus 
incurred.  Hence,  Point  B  is  the  end  of  the  elastic  or  proportional  range  of 
the  material  and  the  value  of  stress  at  Point  B  is  termed  the  "elastic  or 
proportional  limit;  Cp." 

Should  the  stress  be  gradually  increased  up  to  Point  C,  a  noticeable 
yielding  of  the  material  will  be  apparent;  at  Point  C  the  strain  suddenly 
increases  without  further  increase  in  stress.  In  fact,  V7ith  most  ductile 
steels,  there  may  be  a  decrease  in  stress  resisted  by  the  material  as  large 
plastic  strain  takes  place.  It  is  obvious  that  any  stress  above  the  value  at 
Point  C  will  produce  very  large  and  objectionable  permanent  strains.  To 
verify  this  condition,  assume  that  a  stress  is  applied  up  to  Point  X  on  the 
diagram  in  Figure  12.42.  The  strain  at  this  point  is  quite  large.  If  the 
stress  were  to  be  released,  the  material  would  "relax"  along  the  dotted  line 
shovjn  which  is  parallel  to  the  original  straight  line  frcm  B  to  A.  At  Point  Y 
the  stress  is  again  zero  but  a  large  permanent  deformation  has  taken  place. 
The  value  of  stress  at  Point  C  is  logically  termed  the  "yield  point"  or  "Yield 
stress,  material  and  is  the  stress  beyond  which  large  and 
objectionable  permanent  strains  take  place. 

The  stress-strain  diagram  of  Figure  12.42  does  show  that  the  material  is 
capable  of  withstanding  stresses  greater  than  the  yield  stress  -  but  not 
without  large  permanent  strains.  If  the  stresses  above  Point  C  are  gradually 
applied,  the  material  will  continue  to  withstand  higher  and  higher  stresses 
imtil  the  very  ultimate  strength  capabilify  is  reached  at  Point  D.  If  any 
attempt  is  made  to  subject  the  specimen  to  a  stress  greater  than  the  value  at 
Point  D,  failure  will  begin  and  will  be  catplete  at  Point  E.  Since  the  value 
of  stress  at  Point  D  is  the  very  highest  stress  the  material  can  withstand 
without  failxore,  it  is  termed  the  "ultimate  strength"  or  "ultimate  stress, 
a^",  of  the  material. 

The  two  most  iitportant  strength  properties  vMch  are  derived  fran  the 
stress-strain  diagram  are  the  "yield  strength"  and  the  "ultimate  strength." 
There  is  a  direct  analogy  between  these  two  properties  and  the  operating 
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strength  liinitations  of  an  airframe  structure.  If  the  material  shown  on  the 
stress-strain  diagram  is  never  subjected  to  a  stress  above  the  yield  point,  no 
significant  or  objectionable  permanent  strains  will  take  place;  if  an  airframe 
structure  is  never  subjected  to  a  load  condition  greater  than  the  "limit" 
load,  no  significant  or  objectionable  permanent  deformations  will  be  incurred. 
If  the  material  shown  in  the  stress-strain  diagram  were  subjected  to  a  stress 
above  the  yield  point,  large  and  undesirable  permanent  strains  will  take 
place;  if  an  airframe  is  subjected  to  a  load  condition  greater  than  the 
"limit"  load,  undesirable  permanent  deformation  of  the  structure  may  be 
anticipated,  e.g. ,  permanently  distorted  fuselage,  bent  wings,  deformed  tanks, 
etc.  If  an  attempt  is  made  to  subject  a  material  to  a  stress  greater  than  the 
ultimate,  failure  will  then  occur;  if  any  flight  condition  is  attoipted  vdiich 
produces  loads  greater  than  "ultimate"  load,  actual  failure  of  the  airframe 
is  imminent. 


FIGURE  12.42.  PEIMMENT  SET 

The  basic  stress-strain  diagram  of  Figure  12.43  readily  defines  five  of 
the  roost  important  static  strength  properties  of  a  material: 

1.  The  elastic  or  proportional  limit  is  the  end  of  the  elastic  region  of  the 
matarlall  A  part  subjected  to  stresses  at  or  below  the  proportional 
limit  will  experience  no  permanent  deformation.  Upon  the  release  of 
stress  the  part  will  snap  back  to  the  original  vmstressed  shape. 
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2.  Ihe  yield  strength  is  the  hipest  practical  value  of  stress  to  v^ich  a 
material  should  oe  subjected.  Stresses  between  the  proportional  limit 
and  the  yield  point  will  cause  only  slight  and  hardly  measurable 
permanent  deformation.  Any  stress  above  the  yield  point  will  result  in 
large  and  c±>jectionable  permanent  deformation. 

3.  The  ultimate  strength  is  the  very  maximum  of  stress  vMch  a  material  can 
withstand  without  failure.  Extrenely  large  and  undesirable  permanent 
deformations  will  ordinarily  result  viien  this  point  is  approached. 

4.  The  fractiice  point,  or  the  effective  stress  at  time  of  failure,  is 
determined  primarily  to  evaluate  the  manner  of  fracture  and  the  ductile 
quality  of  the  material. 

5.  Ihe  total  strain  or  total  elongation  of  the  material  at  the  point  of 
fracture  is  an  indication  of  the  ductility  of  the  material.  Any  metal 
vhich  would  have  less  than  five  per  cent  elongation  in  a  two  inch  test 
length  is  considered  to  be  brittle  -  ordinarily  too  brittle  to  be 
applicable  to  an  ordinary  aircraft  or  missile  structure. 
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At  this  point,  it  is  appropriate  to  present  seme  of  the  variations  in  the 
stress-strain  diagrams  due  to  manner  of  loading  or  material  type.  One  point 
to  consider  is  that  stress-strain  diagrams  are  not  usually  used  in  connection 
with  shear  properties  of  a  material.  Ihe  proportional  limit,  ultimate 
strength,  etc.,  in  shear  are  not  true  properties  because  of  "section  or  form" 
considerations.  Cross-section  dimensions  and  area  distribution  will  cause 
significant  variations  in  the  shear  strength  capabilities. 

Many  materials  used  in  aircraft  construction  (aluminim  alloys,  magnesium 
alloys,  and  seme  steels)  do  not  exhibit  a  definite  yield  point  or  a  distinct 
proportional  limit.  Figure  12.44  illustrates  this  fact.  In  such  an  instance 
it  is  necessary  to  define  the  yield  point  and  proportional  limit  as  a  given 
departiire  from  the  original  straight  line  of  0.0001  in/in.  The  yield  stress 
is  determined  as  the  stress  v^ich  produces  a  departure  frem  the  original 
straight  line  of  0.002in/in  (Figinre  12.45).  In  this  case  an  applied  stress 
equal  to  the  yield  stress  woxild  result  in  a  permanent  set  of  0.2%.  This 
amount  is  considered  admissible  for  ordinary  purposes. 


Corrpression  stress-strain  diagrams  are  similar  to  tension  stress-strain 
diagrams  except  that  the  departure  frem  proportionality  generally  occurs 
sooner  and  more  gradually.  Catpression  stress-strain  diagrams  are  more 
difficult  to  obtain  correctly  because  of  buckling  of  the  specimen.  As 
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FIGURE  12.45.  DEFINITION  FOR  VARIANT  STRESS-STE^AIN  REILATICWSHIP 


cctipression  buckling  of  a  specimen  constitutes  a  failure,  it  is  obvious  that 
compression  iiltimate  strengtli  is  not  true  property  and  could  not  be  determined 
as  a  specific  quality  of  a  metal.  The  stiffness  of  a  metal  and  the  physical 
arrangement  of  the  structure  carbine  to  decide  the  stability  of  the  structure 
vhen  subjected  to  compression  loads  (Figure  12.46). 
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In  order  to  determine  other  material  properties  ej^ressed  by  the 
stress-strain  diagrams  a  closer  examination  must  be  made  of  certain  areas  of 
the  stress-strain  diagram. 

Ely  an  inspection  of  the  straight  line  portion  of  the  stress-strain 
diagram,  an  evaluation  may  be  made  of  the  inherent  stiffness  of  the  material. 
In  the  elastic  range  of  a  material  the  proportion  between  stress  and  strain  is 
the  Modulias  of  Elasticity  (sonetimes  called  "Young's  Modulus")  and  the 
magnitude  of  this  proportion  is  a  direct  measure  of  the  inherent  stif&iess.  A 
high  value  for  the  Modulus  of  Elasticity  will  indicate  a  very  stiff  or  rigid 
material  while  a  low  value  indicates  low  stiffness  or  greater  flexibility 
(Figure  12.47). 


FIGURE  12.47.  DEFINITION  OF  MODULUS  OF  ELASTICITY 

When  the  stress-strain  diagrams  for  three  different  materials  are 
ccnpared,  the  difference  in  the  slopes  and  the  proportions  of  inherent 
stiffness  are  readily  apparent  as  shown  in  Figure  12.48. 

Only  two  solutions  exist:  (1)  lower  the  value  of  the  operating  stress  or 
(2)  change  to  a  material  type  vdiich  has  a  hi^er  elastic  stiffness  (higher  E) . 
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FIGURE  12.48.  MIXOLES  OF  ELASTICITY  COMPARISCN 

One  iirportant  point  to  consider  is  that  the  Modulus  of  Elasticity,  E, 
cannot  be  altered  by  heat  treatment  or  changed  any  significant  amount  by  use 
of  a  different  alley.  Thus,  the  Modulus  of  Elasticity  is  an  intrinsic 
property  of  the  t^^  of  metal.  Figure  12.49  shows  the  typical  stress  strain 
diagrams  for  steel  in  various  conditions  of  heat  treatment.  Notice  that  the 
origin  of  each  has  the  same  slope.  While  the  strength  and  ductility  are 
changed  by  heat  treatment,  the  stiffness  of  the  elastic  material  remains 
unaltered.  Increasing  hardness  by  heat  treatment  will  sinply  increase  the 
stress  at  vdiich  plastic  flow  begins.  Thus,  if  elastic  deflections  of  a  part 
viiich  are  excessive  alreac^  exist,  the  problem  will  not  be  solved  by  heat 
treatment  of  changing  alley. 
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FIGURE  12.49.  EFFECTS  OF  HEAT  TREAIMENT 

Actiially,  the  elastic  modulus  of  a  given  metal  vd.ll  vary  only  v/ith 
tenperature,  elevated  tenperature  producing  a  lower  Modulus  of  Elasticity. 
For  exaitple,  a  hi^  strength  aluminum  alloy  at  600°  F  vdll  exhibit  a  Modulus 
of  Elasticity  V7hich  is  only  one-half  the  value  shewn  at  room  terrperature. 

Hie  elastic  range  of  a  material  is  the  area  of  most  general  interest  for 
ordinary  stractirral  investigation.  However,  vhen  imestigating  the  phenctnenon 
of  buckling  of  short:  columns  and  the  behavior  of  structures  at  loads  near 
ultimate,  the  stiffness  of  a  metal  in  the  plastic  range  is  of  particular 
interest.  As  was  previously  noted,  the  proportion  between  stress  and  strain 
is  defined  as  the  Modulus  of  Elasticity.  At  any  stress  below  the  proportional 
limit,  this  proportion  is  one  fixed  and  constant  value.  Above  the 
proportional  limit,  this  proportion  is  one  fixed  and  constant  value.  Above 
the  proportional  limit,  this  definition  results  in  a  proportion  vhich  is  knewn 
as  the  "secant  Modulus".  Once  beyond  the  proportional  limit,  the  proportion 
between  stress  and  strain  will  noticeably  decrease.  Hence,  the  Secant  Modulus 
vd.ll  have  a  value  lower  than  the  Modulus  of  Elasticity.  Another  method  of 
measuring  the  stiffness  in  the  plastic  range  is  the  slope  of  a  line  drawn 
tangent  to  the  stress-strain  diagram  at  sane  stress  above  the  proportional 
limit.  The  slope  of  this  tangent  line  defines  a  value  knavn  as  the  "Tangent 
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Modulus".  Cie  Tangent  Modulus  then  measures  an  "instantaneous"  stiffness 
while  the  Secant  Modulus  Ineas^Jres  a  "gross"  or  "cumulative"  stiffness.  Figure 
12.50  gives  the  procedure  of  calculation  of  these  properties  and  Figure  12.51 
illustrates  the  typical  variations  of  the  Secant  and  Tangent  Modulus  for  an 
aluminum  alloy. 

The  loss  of  stiffness  in  the  plastic  range  is  appreciated  when  it  is 
realized  that  small  changes  in  stress  will  produce  larger  dianges  in  strain 
than  in  the  elastic  region. 

The  energy  storing  and  energy  absorbing  characteristics  are  of  particular 
iirportance  in  defining  the  properties  of  a  material.  Any  stress  produced  in  a 
material  requires  the  application  of  a  force  through  a  distance  and,  as  a 
consequence,  a  certain  amount  of  work  is  done.  Consider  a  cubic  inch  elenent 
of  steel  with  a  gradually  applied  stress  of  30,000  psi  (as  in  Figure  12.52). 
As  the  stress  is  gradually  increased  from  zero  to  30,000  psi,  the  strain 
gradually  increases  fran  zero  to 

o  _  o  _  30,000 
^  E  30,000,000  “ 


FIGURE  12.50.  DEFINITION  OF  SECANT  AND  TANGENT  MODULI 
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FIGURE  12.51.  VARIATION  OF  SECANT  AND  TANGENT  MODULI 


FIGURE  12.52.  APPLICATION  OF  STE^S 


The  work  done  during  this  gradual  stressing  is  the  product  of  the  average 
force  and  the  distance. 
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Thus,  15  inch-pounds  of  work  would  be  required  to  produce  the  final  stress  of 
30,000  psi  in  the  cubic  inch  element.  This  amount  of  work  is  actually 
represented  by  the  area  enclosed  on  the  stress-strain  diagram  at  this 
particular  stress  level  (Figure  12.53) .  This  area  principle  is  then  used  with 
the  stress-strain  diagram  to  describe  the  energy  properties  of  a  material. 

If  a  material  is  stressed  to  the  proportional  limit,  all  work  done  in 
producing  this  stress  is  stored  elastically  in  the  material.  Therefore,  the 
area  under  the  straight  line  portion  of  the  stress-strain  diagram  is  a  direct 
measure  of  the  energy  storing  capability  and  is  referred  to  as  the  "Modulus  of 
Resilience".  Of  course,  the  manner  of  loading-tension,  coipression,  shear, 
etc. ,  must  be  specified  since  the  work  done  will  be  different  depending  on  the 
manner  of  loading. 


FIGURE  12.53.  WORK  DC^E  BY  STRESS  APPLICATION 


If  a  material  vrere  stressed  all  the  way  to  failure,  the  entire  area  under 
the  stress-strain  diagram  is  a  measure  of  the  work  required  to  fail  the 
material.  The  tou^iness  or  energy  absorbing  quality  of  the  material  is 
evaluated  in  this  manner. 


FIGURE  12.54.  ENERGY  STORING  CAPABILITIES  (RESILIENCE) 
OF  DUCTILE  AND  BRITTLE  STEELS 


A  catparison  of  the  stress-strain  diagrams  for  a  ductile  and  brittle 
steel  should  define  the  properties  of  "resilience"  and  "toughness".  The 
shaded  areas  shewn  in  Figure  12.54  denote  the  energy  storing  capabilities  for 
the  ductile  and  brittle  steels.  The  brittle  material  has  the  higher  elastic 
limit  and  consequently  a  higher  "Modulus  of  Resilience".  The  shaded  areas  of 
Figure  12.55  denote  the  energy  absorbing  capabilities  of  the  same  materials. 
The  ductile  material,  while  having  lower  strengths,  devlops  much  greater 
strains.  As  a  result,  the  ductile  material  will  require  a  greater  amount  of 
vrork  to  produce  failure  and  is  then  the  tougher  material.  The  additional  work 
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required  to  break  the  softer,  nore  ductile,  material  goes  into  forcing  the 
metal  particles  to  slip  relative  to  each  other.  When  a  specimen  of  ductile 
material  is  fractured  during  a  laboratory  test,  this  effect  may  be  appreciated 
by  handling  the  broken  specimen  of  ductile  material  inmediately  after  failure. 
A  ductile  specimen-  will  be  warn  to  the  touch  as  the  work  absorbed  by  the 
material  is  converted  into  heat. 


FIGURE  12.55.  ENER3Y  ABSORBING  CAPABILITIES  (TOUGHNESS) 
OF  DUCTILE  AND  BPJTTLE  STEELS 


It  seem  strange  that  it  is  possible  to  fail  the  more  ductile  metal  at 
a  losi^er  stress  but  still  require  that  more  work  be  done.  If  so,  remember  that 
vrork  is  the  product  of  average  force  and  distance. 

It  wDuld  be  desirable  for  a  structural  material  to  have  a  high  yield 
point,  high  resilience,  hi^  ultimate  strength,  and  also  high  tou^ness. 
However,  materials  vdiich  have  very  high  strength  usually  have  low  ductility 
and  lew  tou^iness.  The  balance  between  the  strength  requirements  and 
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toui^ess  (or  energy  absorbing)  requirements  will  depend  on  the  particular 
application  in  a  structure.  As  examples  of  the  two  extremes  consider;  (1)  a 
reciprocating  engine  valve  spring  and  (2)  protective  crash  helmet  or 
protective  headgear.  In  the  fabrication  of  a  valve  spring,  a  material  must  be 
selected  vAiich  has  .very  hi^  strength  and  great  resilience.  Such  a  material 
would  necessarily  have  low  tou^iness  but  in  such  an  ^plication  toughness  is 
unimportant  and  resilience  is  given  primary  consideration.  In  the 
construction  of  a  true  protective  crash  helmet  or  hard  hat,  sharp  impact  blows 
first  must  be  distributed  then  the  energy  of  impact  absorbed  as  far  as  is 
possible.  This  requires  a  thickness  of  a  crushable  material  for  the  inner 
lining  vdiich  has  high  tou^iness  per  unit  wei^t.  Actual  strength  is  not 
necessarily  a  factor  since  energy  absorption  requirements  predoninate. 

The  stress-strain  diagram  of  a  material  will  then  furnish  all  necessary 
information  to  determine  the  static  strength  properties.  Actually  there  are 
10  important  points  of  information  that  may  be  gained  from  an  inspection  of 
the  stress-strain  diagram.  As  shown  in  Figure  12.56  these  points  are  as 
follows: 


FIGURE  12.56.  STATIC  STRENGTH  PROPERTIES 
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(1)  The  proportional  limit  is  the  stress  v^ch  denotes  the  end  of 
proportionality  oetween  stress  and  strain.  If  not  clearly  defined, 
0.0001  in/ in  offset  is  used. 

(2)  The  elastic  range  of  stress  is  the  range  of  stress  and  strain  up  to  the 
proportional  limit. 

(3)  The  Modulus  of  ^asticity  is  the  slope  of  the  straight  line  portion  of 
the  stress-strain  diagram.  This  slope  measures  inherent  stiffness. 
Beyond  the  proportional  limit,  the  Secant  or  Tangent  Modulus  will  be 
appropriate. 

(4)  The  resilience  is  measured  by  the  area  xmder  the  straipt  line  portion  of 
the  stress-strain  diagram.  This  indicates  ability  to  store  energy 
elastically. 

(5)  The  yield  stress  is  the  value  of  stress  above  vdiich  objectionable  amounts 
of  permanent  strain  are  incurred.  If  not  clearly  defined,  0.002  in/in 
offset  is  used. 

(6)  Ihe  ultimate  strength  is  the  very  hipest  stress  that  the  material  can 
withstand  without  failure.  This  represents  the  maximum  load-carrying 
c^ability  for  static  loads. 

(7)  The  plastic  range  of  stress  is  between  the  proportional  limit  and  the 
ultimate  stress.  Permanent  strains  occur  in  this  area.  Below  the  yield 
stress  these  permanent  strains  are  relatively  small  and  insignificant; 
above  the  yield  stress  these  permanent  strains  are  large  and 
objectionable. 

(8)  The  fracture  point  is  the  effective  stress  at  time  of  failure.  It  is 

noted  primarily  to  evaluate  the  manner  of  failure  and  the  quality 

of  the  material. 

(9)  The  total  elongation  is  a  measure  of  ductility.  If  less  than  five 
percent  in  a  two  inch  ^cimen  length,  the  material  is  considered 
brittle. 

(10)  Ihe  tou^iness  of  a  material  is  represented  by  the  total  area  \mder  the 
stress— strain  diagram.  This  indicates  the  amount  of  work  required  to 
fail  the  material  and  denotes  energy  absorption  capability. 


An  interesting  phenomenon  in  connection  with  the  stress-strain  diagram  is 
"work  hardening"  or  "strain  hardening".  Suppose  that  a  material  is  subjected 
to  a  stress  beyond  the  yield  point  and  then  released.  If  stress  is 
subsequently  reapplied,  the  new  stress-strain  diagram  will  be  a  straight  line 
up  to  the  point  where  stress  was  released  and  then  it  continues  along  the 
original  curve.  Figure  12.57  illustrates  this  process.  Ihe  material  vhich 
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FIGURE  12.57.  WORK  HARDENING 


has  been  permanently  stretched  will  then  have  higher  proportional  and  yield 
strengths.  I'iany  of  the  topical  aircraft  materials  -  aluminum  alloy's, 
stainless  steels /  etc.  —  are  work  hardened  to  produce  these  beneficial  gains 
in  strength  properties.  Of  course,  the  v/crk  hardening  must  be  limited  to 
prevent  loss  of  ductility  and  the  formation  of  flaws  and  fissures. 

Ihe  effects  of  ductility  and  the  plastic  range  of  stress  are  quite 

significant  in  predicting  the  failing  load  of  a  structure.  The  ultimate 

stress  has  been  defined  previously  at  the  maximum  stress  a  material  will 

withstand  without  failing.  This  stress  is  designated  as  the  "tension 

ultimate"  )  or  "shear  ultimate"  (a  )  depending  on  the  manner  of  loading 
u  u 

and  is  determined  by  tests  of  small  specimens.  If  these  values  of  strength 
are  used  to  predict  the  ultimate  strength  capability  of  large  sections  in 
bending  and  torsion,  noticeable  errors  may  result.  An  elastic  stress 
distribution  in  bending  may  be  predicted  by  the  following  relationships; 


a 


b 


I 
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vtere  Oj^  =  bending  stress,  Ib/in^ 

M  =  applied  bending  itranent,  in-lb 

y  =  element  distance  frcm  neutral  axis,  in 

I  .  =  mcment  of  inertia  of  the  cross-section,  in^ 

Such  stress  distribution  is  linear  -  varying  directly  with  the  element 
distance  frcm  the  neutral  axis,  y.  A  typical  elastic  stress  distribution  is 
shown  in  Figure  12.58.  ^-Jhenever  stresses  are  produced  vhich  are  beyond  the 
proportional  limit  of  the  material,  the  bending  stress  distribution  tends  to 
remain  linear  but  due  to  the  loss  of  proportionality  between  stress  and  strain 
in  the  plastic  range,  the  stress  distribution  is  non-linear.  Figure  12.59 
shows  a  typical  stress  distribution  resulting  from  bending  loads  which  create 
stresses  in  excess  of  the  proportional  limit. 

Notice  that  for  a  given  maximum  stress  at  the  outer  fiber,  the  inelastic 
bending  stress  distribution  requires  a  greater  applied  bending  moment  than  the 
elastic  stress  distribution.  Ohe  reason  for  this  is  that  the  outer  fibers 
begin  to  yield  allowing  the  underlying  fibers  to  develop  a  stress  higher  than 
predicted  by  elastic  theory.  Thus,  the  use  of  the  equation 


FIGURE  12.58.  ELASTIC  STRESS  DISTRIBUTIOT 
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FIGURE  12.59.  PLASTIC  (NON-LINEAR)  STRESS  DISTRI3UTICH 

Me 

%  =  T 

to  conpute  the  maximum  bending  stress  is  valid  only  for  maximuiri  stresses  vAiich 
do  not  exceed  the  proportional  limit.  In  order  to  predict  the  failing  load  of 
a  structiare  in  bending,  the  same  form  of  equation  may  be  used  with  a 
fictitious  ultimate  stress  referred  to  as  the  "bending  modulus  of  rupture". 
This  bending  modulus  of  rupture  is  defined  by  the  following  equation: 

MbC 

I 

2 

bending  modulus  of  rupture,  Ib/in 
bending  monent  to  cause  failure,  in-lb 
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c  =  distance  to  critical  outermost  element,  in. 

I  =  section  moment  of  inertia,  in^ 

If  the  material  is  ductile  and  no  buckling  of  the  section  occurs,  the 
bending  modulvis  of-  rupture  will  be  seme  value  greater  than  the  tensile 
ultomate  strength.  In  a  perfectly  ductile,  stable  cross~section  the  bending 
modulus  of  rupture  could  be  1.5  times  the  tensile  ultimate;  in  a  very  brittle, 
stable  cross-section  the  bending  modulus  of  rupture  would  be  equal  to  the 
tensile  ultimate.  If  the  cross-section  is  canposed  of  thin  walled  unstable 
elements  which  are  apt  to  buckle,  the  bending  modulias  of  rupture  may  be  much 
less  than  the  tensile  ultimate  strength.  The  bending  modulus  of  rupture  is 
obviously  dependent  upon  the  type  of  material  (especially  ductility)  and  the 
shape  or  form  of  the  cross-section. 

An  analogoias  situation  exists  for  sections  ST±»jected  to  torsion  loading 
and  is  due  to  the  form  of  the  cross-section  and  the  character  of  the  material. 
In  order  to  predict  the  failing  torsion  load  of  a  structure,  another 
fictitious  stress  is  used  vMch  is  referred  to  as  the  "torsion  inod\ilus  of 
rupture" ,  .  This  torsion  modulus  of  n5>ture  is  defined  by  the  following 
equation; 


=  torsion  modulus  of  nature,  psi 
M 

tj^  =  torque  or  moment  to  cause  failure,  in-lb 
c  =  distance  to  critical  outermost  element,  in 
J  =  polar  moment  of  inertia  of  cross-section  in^ 

If  a  metal  is  subjected  to  a  stress  in  the  plastic  range,  the  strain  will 
continue  to  change  in  the  direction  of  stress  althou^  the  stress  is  held 
constant.  This  is  a  condition  vdiich  is  most  cannon  to  metals  at  high  stresses 
and  hi(^  tatperatures  and  the  phenomenon  is  known  as  "creep".  Figure  12.60  is 
a  "tg^itial  creep  curve  for  a  metal  at  a  constant  stress  and  teitperature. 
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FIGURE  12.60.  CREEP  CURVE  FOR  A  METAL  AT  CONSTANT  STRESS  AND  TEMPERATORE 

Upon  the  application  of  stress,  a  hi^  initial  creep  rate  will  develop, 
then  decrease  to  the  irimimum  value  of  steady  creep.  Daring  the  second  stage 
the  creep  rate  is  essentially  constant.  After  a  period  of  time  the  third  and 
final  stage  of  creep  will  take  place  with  an  increase  in  creep  rate  and  final 
rupture.  Since  the  stress  to  cause  rupture  varies  inversely  with  time  and 
terrperatxxre,  creep  problems  are  of  the  greatest  iitportance  in  light  weight, 
hi^  tenperature  structtires. 

In  the  design  of  hi^  tenperature  structures  the  amount  of  deformation 
allowable  may  be  a  more  severe  criteria  than  actual  rupture  strength.  This 
woxild  be  the  case  for  coiponents  v^ch  -  if  excessively  strained  -  would  not 
function  properly  or  would  fail  at  loads  lower  than  normally  anticipated. 

Stress-strain  diagrams  for  a  material  may  be  greatly  altered  when  the 
rate  of  stressing  is  very  high.  When  stress  is  applied  very  suddenly  (almost 
instantaneously)  the  result  is  "iiipact"  stresses.  Ihe  amount  of  energy 
absorbed  under  inpact  conditions  may  be  significantly  different  from  the 
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energy  absorbed  vjhen  the  load  is  applied  steadily  and  gradually.  Ihe  actual 
speed  or  rate  at  \diich  a  material  is  stressed  will  determine  what  changes  in 
energy  absorption  tate  place.  In  the  case  of  an  ordinary  structural  material, 
an  increase  in  the  rate  of  stressing  (above  that  of  very  gradual  load 
application)  will  initially  produce  a  slight  increase  in  the  energy  absorbed. 
With  continued  increase  in  rate  of  stressing,  a  "critical  speed"  will  be 
reached  and  the  energy  absorbed  will  be  at  maximum.  Above  this  speed  the 
toughness  will  be  greatly  reduced.  If  the  stress-strain  diagrams  for  impact 
stressing  v^re  recorded  and  ccrpared  with  the  stress— strain  diagram  obtained 
frcm  gradual  load  application,  the  result  would  be  similar  to  Figure  12.61. 


FIGURE  12.61.  EFFECTS  OF  APPLICATION  RATE  ON  STRESS 

When  stress  is  applied  above  the  critical  speed,  the  failing  stress  would  be 
higher  but  the  elongation  would  be  greatly  reduced.  In  fact,  the  effect  of 
impact  stresses  beyond  the  critical  speed  is  to  produce  brittle  type  failures 
in  tough  or  ductile  materials.  On  the  other  hand,  a  brittle  material  will  not 
show  any  great  effect  of  high  stress  rates  since  there  is  very  little  ductility 
or  energy  absorbing  capability  for  gradually  applied  loads. 
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While  ijipact  stresses  near  the  critical  speed  are  not  ordinarily 
encountered  in  airframe  structures,  due  consideration  most  be  given  to  the 
case  of  dynamic  mchinery  and  mechanisms.  The  effect  of  irrpact  stresses  is 
most  important  vdien  there  are  severe  discontinuities  in  the  shape  of  a  part  or 
vdien  a  part  is  operated  at  low  terrperatxures. 

Figures  12.62-12.69  are  illustrations  of  failure  of  different  types  of 
materials  under  tension,  torsion  or  shear. 
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FIGUEE  12.63.  MEDIUM  DUCTILITy  TENSION  FAILURE 
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ROUGH  GRANULATED 
TENSION  TYPE  ZONE 


FIGURE  12.64.  HIGHLY  DUCTILE  TENSICN  FAILURE  SHEET  OR  THIN  BAR  STOCK 
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DUCTILE 


SHEAR  TYPE  CIRCLES 


FIGURE  12.66.  TORSION  FAILURE 
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SMOOTH  SURFACE 


DUCTILE 


ROUGH  GRANULAR 
TENSION  TYPE  ZONE 


FIGUPvE  12.67, 


SHEAR  FAILURE 


45'  CRACK  DUE  TO 
RESOLVED  TENSION 


FIGURE  12.69.  SHEAR  IN  PANELS  TENSION  TYPE  FAILURE 
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12.5  AERDEIASTICITY 

12.5.1  INTRODUCTICN  AND  DEFINITIONS 


Except  for  a  few  isolated  exairples  in  the  TPS  curriculxm,  the  aircraft 
has  been  treated  as  a  rigid  bo<^.  That  is,  no  bending,  twisting,  or 
defoinnations  were  assumed  to  occur  on  the  structure  in  the  definition  and 
derivation  of  the  performance  and  stability  characteristics.  For  exairple, 
lateral  bending  of  the  fuselage  was  ignored  in  determining  the  steacfy  straight 
sideslip  equations.  Or  the  effects  of  wing  bending  and  subsequent  load 
redistribution  were  not  taken  into  account  in  the  lift,  drag,  and  mcment 
equations. 

In  reality,  thou^,  airplanes  deform  under  aerodynamic  loads.  Ihese 
effects  can  be  significant,  especially  in  high  speed  lightweight  aircraft  at 
hi^  maneuvering  accelerations,  or  in  the  large  category  aircraft  such  as  the 
C-5A,  Boeing  747  or  Lockheed  L-1011  during  nonral  flight  conditions.  The 
degree  to  w^ch  the  analytical  conputations  are  affected  is  dependent  upon 
flight  conditions.  However,  aeroelastic  effects  must  be  included  in  any 
precision  analysis.  We  therefore  must  treat  the  stability  derivatives  as 
functions  of  dynamic  pressure  (q) ,  as  well  as  Mach  (M)  and  angle  of  attack 
(a) .  Designing  and  building  a  totally  rigid  airplane  is  not  practical  because 
of  the  wei^t  penalty  (the  development  of  lightweight  caiposites  has  lessened 
the  problem  sonewhat  but  high  cost  becones  a  major  factor  in  that  area) .  The 
designer  is  thus  faced  with  the  dilertma  of  conflicting  requironents  for 
lightweight  (for  iitproved  aircraft  performance)  versus  structinral  rigidity  (to 
preclude  aeroelastic  effects) . 

The  study  of  aeroelasticity  is  important,  therefore,  simply  because  the 
aeroelastic  problem  is  a  reality.  A  decision  has  to  be  made  as  to  how  much 
rigidity  can  be  sacrificed  before  the  bending,  twisting  and  inertial  effects 
of  the  structure  restrict  the  performance  and  handling  qualities.  Is  there 
assurance  that  catastrophic  structural  failure  will  not  occur  somevdiere  in  the 
aircraft's  flight  regime? 

Aeroelasticity  is  often  defined  as  the  science  which  stvdies  the  inutual 
interaction  between  aerodynamic,  elastic  and  inertial  forces  of  an  airplane  in 
fli^t.  Again,  the  phenomenon  would  be  nonexistent  if  aircraft  structures 
v;ere  perfectly  rigid  but  the  weight/cost  penalties  for  that  privilege  would  be 
too  severe.  By  itself,  structural  bending  or  flexibility  is  not 
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objectionable.  The  problems  arise  viien  the  deformations  in  turn  cause  changes 
in  the  aercx^ainic  forces.  If  the  deformations  and  aerodynamic  forces  vary 
rapidly,  inertial  forces  become  iirportant. 

In  1946,  A.R.  Collar  presented  a  paper  to  the  Royal  Aeronauticcil  Society 
that  ingeniously  classified  problems  in  aeroelasticity  by  means  of  a  triangle 
of  forces.  Referring  to  Figure  12.70,  the  three  types  of  forces  (aerodynamic, 
elastic,  and  inertial)  are  placed  at  the  vertices  of  a  triangle.  Each 
aeroelastic  ftenomena  can  be  located  on  the  diagram  according  to  its  relation 
to  the  three  vertices.  For  exaitple,  dynamic  aeroelastic  phenomena  such  as 
flutter,  F,  lie  within  the  triangle,  since  they  involve  all  three  types  of 
forces  and  must  be  bonded  to  all  three  vertices.  Static  aeroelastic  phenomena 
such  as  wing  divergence,  D,  lie  outside  the  triangle  on  the  iipper  left  side, 
since  they  involve  only  aeixx^amic  and  elastic  forces.  Althou^  it  is 
difficult  to  define  parecise  limits  in  aeroelasticity,  the  classes  of  problems 
connected  by  solid  lines  to  the  vertices  are  usually  accepted  as  principsal 
ones.  Of  course,  other  borderline  fields  can  also  be  placed  on  the  diagram. 
For  exaiiple,  mechanical  vibrations,  V,  and  rigid  body  aerodynamic  stability, 
DS,  are  connected  by  dotted  lines.  It  is  very  likely  that  in  certain  cases 
the  dynamic  stability  problem  is  influenced  by  aircraft  flexibility  and  it 
would  therefore  be  moved  within  the  triangle  to  correspxond  with  DAS,  vhere  it 
would  be  regarded  as  a  dynamic  aeroelastic  problem. 

Collar's  Aeroelastic  Triangle  (or  more  ccitpletely  "Aeroinertia-Elastic" 
Triangle)  was  revised  to  a  tetrahedron  some  15  years  later  by  I.E.  Garrik  to 
include  Aerotherrooelasticity,  Figure  12.71.  Aeroc^amic  thermal  effects 
associated  with  hi^-sp)eed  fli^t  vehicles  introduce  deformations,  streses, 
and  changes  in  material  properties  that  can  greatly  extend  the  field  of 
aeroelasticity.  This  chapter  will  not  formally  address  Aerothermoelasticity, 
except  to  introduce  it  as  an  influencing  factor  of  aeroelasticity  in 
hi^-sp)eed  fli^t. 

Static  jteroelastic  Ihencmena; 

Ihese  phenomena  are  the  ones  vdiich  are  concerned  with  steacfy  state  (i.e. 
non-oscillatory)  aerodynamic  loads,  and  the  associated  steady  state 
distortion.  Static  implies  the  absence  of  inertial  forces,  thus  it  requires 
only  coupling  between  the  aerodynamic  and  elastic  forces  of  the  aircraft.  It 
includes  the  following  classical  cases: 
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A:  AERODYNAMIC  FORCE 
E:  ELASTIC  FORCE 
I:  INERTIAL  FORCE 

RELATED  FIELDS 

V:  MECHANICAL  VIBRATIONS 
DS;  DYNAMIC  STABILITY 


AEROELASTIC  PHENOMENA 

F:  FLUTTER 

B:  BUFFETING 

Z:  DYNAMIC  RESPONSE 

L:  LOAD  DISTRIBUTION 

D:  DIVERGENCE 

C:  CONTROL  EFFECTIVENESS  . 

R:  CONTROL  SYSTEM  REVERSAL 

DSA:  AEROELASTIC  EFFECTS  ON  DYNAMIC  STABILITY 

SSA:  AEROELASTIC  EFFECTS  ON  STATIC  STABILITY 


FIGURE  12.70.  CQIIAR'S  AEROELASTIC  TRIANGLE  OF  FORCES 
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A:  AERODYNAMIC  FORCE 
E:  ELASTICITY  FORCE 
I:  INERTIA  FORCE 
H:  HEATING 


FIGURE  12.71.  TETRAHEDRON  OF  AEROTHERMOELASTICITy 


^oelastic  Effects  on  Static  Stability,  SSA.  As  the  aerodynamic  and 
elastic  forces  ccnibine,  the  deformations  of  the  structure  can  be  of 
siafficient  magnitude  to  influence  the  static  stability  derivatives. 

load  Distribution,  L.  Influence  of  elastic  defonnation  of  the  structure 
on  the  distribution  of  aerodynamic  pressures  ov^er  the  structure. 

<^ntrol  Effectiveness,  C.  This  phenomsna  attempts  to  evaluate  the 
influence  of  the  elastic  deformations  on  the  controlability  of  the 
aircraft.  • 

Torsional  Divergence,  D.  A  static  instability  of  a  lifting  surface  at  a 
spe^  the  divergence  speed  vdiere  the  elasticity  of  the  lifting 

surface  plays  an  essential  role  in  the  instabilitv.  At  a  speed  slightly 
above  the  wing  divergence  speed,  the  elastic  restoring  mcfnent  about  a 
spanwise  with  elastic  axis  can  no  longer  balance  the  aerodynamic  irorent 
created  by  the  airloads. 

Control  Reversal,  R.  A  condition  occurring  in  flight  at  a  speed  called 
the  control  reversal  speed,  \diere  the  intended  effects  of  displacing  a 
control  surface  are  nullified  by  the  elastic  deformations  of  the 
structure.  For  example,  a  right  roll  aileron  deflection  may  result  in  a 
roll  to  the  left  because  of  a  wing  twist  caused  by  the  deflected 
®il®^ons.  Ihe  B-47  for  example  exhibited  an  aileron  control  reversal 
Edienonena  at  .80  Mach. 

Dynamic  Aeroelastic  Phenctnena; 

In  dynamic  aeroelastic  prcbleros,  we  are  concerned  with  the  oscillatory 
motion  of  various  parts  of  aircraft,  and  particularly  interested  in  the 
conditions  under  vAiich  these  oscillatory  modes  tend  to  diverge  (increase  in 
amplitude) ,  because  this  may  result  eventually  in  the  structural  failure.  As 
mentioned  previously,  ^amic  aeroelastic  phencnena  are  the  interaction 
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between  aeroc^amic,  inertial,  and  elastic  forces.  Sate  exaitples  of  this 
problan  are: 


Flutter,  F.  A  self-excited  t^amic  instability  of  the  structural 
coiponents  of  an  aircraft,  usually  involving  the  coupling  of  separate 
vibra-^on  modes,  vAiere  the  forcing  function  for  oscillation  is  drawn  fron 
the  airstream.  The  coupling  of  the  bending  and  twisting  mndes  of  a  wing 
results  in  a  b^ding-torsional  wing  flutter.  Other  examples  of  flutter 
modes  are:  wing-aileron,  tail- fuselage,  and  bending-torsion-aileron. 
Sinply  stated,  flutter  is  the  dynamic  instability  of  an  elastic  body  in 
an  airstream.  Flutter  speed  Up  and  corresponding  frequency  are  dpf j  np>d 

as  the  lowest  airspeed  and  frequency  vdiere  a  flying  structure  will 
exhibit  sustained,  simple  harmonic  oscillations. 

Buffeting,  B.  Transient  vibrations  of  aircraft  structural  coiponents  due 
to  aerodynamiic  impulses  produced  by  the  wake  behind  wings,  nacelles, 
fuselage  pods,  or  other  components  of  the  airplane.  The  problem  can  be 
serious  in  filter  aircraft  during  maneuvering  to  CL  at  hi^  speed, 

max 

often  resulting  in  rugged  transient  vibrations  on  the  tail  due  to  aero- 
<^amic  impulses  from  the  wing  wake.  Since  these  impulses  are  quite 
randan  there  is  no  analytic  theory  to  adequately  describe  the  phenomenon, 
Ciires  are  usually  made  as  necessary  by  prcper  positioning  of  the 
assembly. 

Aeroelastic  Effects  on  Dynamic  Stability,  PSA.  Changes  in  an  aircraft 
dynamic  stability  can  result  due  to  elastic  bonding,  twisting,  or 
deformation  of  its  structure  (e.g. ,  the  changes  in  the  short  period 
frequency  and  damping  due  to  fuselage  bending) . 

Mechanical  Vibration,  V.  A  mechanical  vibration  of  the  airplane 
structure  mainly  due  to  the  coincidence  of  a  powerful  engine  harmonic 
with  an  airframe  frequency.  Aeroctynamic  interaction  m^  not  be 
necessary. 

It  should  be  noted  from  the  above  discussion  that  flutter  and  divergence 
corresponded  to  conditions  of  aeroelastic  instability,  and  that  speeds 
b^ond  the  critical  flutter  and  divergence  speeds  will  result  in  an  eventual 
structural  failure.  However,  control  reversal  is  not  a  condition  of 
instability,  and  speeds  beyond  control  reversal  speed  will  result  only  in  a 
reversal  of  the  action  of  the  control  system  and  not  necessarily  in  a  failure 
of  the  structure.  Control  effectiveness  influences  maneuverability,  and  it  is 
important  that  the  control  system  designer  has  a  thorough  understanding  of 
this  phenanena. 

Thus,  the  three  critical  airspeeds  involving  aeroelasticity  are  flutter 
(U^) ,  divergence  (U^) ,  and  aileron  reversal  (Uj^)  speeds.  In  the  early  design 
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stages  of  an  aircraft,  the  ccrparative  values  of  these  speeds  must  be 
caipletely  analyzed.  Figure  12.72  shows  the  relation  between  the  critical 
speeds  for  a  typical  wing  with  varying  amounts  of  forward  and  backward  sweep. 


FIGURE  12.72.  EFTECT  CF  WING  SWEEP  ON  CRITICAL  SPEEDS 
12.5.2  Historical  Backgroxmd 

Aeroelastic  prcblems  were  relatively  unknown  until  World  War  II.  Prior 
to  that  time,  aircraft  speeds  vere  relatively  low  and  the  load  reguirements 
and  lack  of  design  refinements  produced  an  aircraft  structure  rigid  enou^  to 
preclude  most  aeroelastic  phenomena.  Except  for  the  time  spans  of  the  two 
World  Wars,  aircraft  top  speeds  had  increased  approxiraately  19  knots  per  year 
from  1910  to  1955.  As  the  speed  of  aircraft  began  to  increase  with  no 
increase  in  load  or  stiffeess  requirements,  the  designers  and  pilots  began  to 
encounter  problems  associated  with  aeroelasticity. 

When  airplanes  were  first  built  there  were  no  logical  stiffness  criteria 
for  design.  Hence,  increasing  speeds  led  to  the  wide  variety  of  aeroelastic 
problems.  Control  surface  flutter  first  occurred  at  about  110  knots;  wii^ 
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flutter  started  at  around  twice  that  speed.  Ihe  strength  questions  posed  by 
the  speed  increase  v^e  solved  partly  by  material  developments  and  partly  by 
constructional  techniques,  but  not  before  mar^  inr-f  light  failures  had 
occurred. 

Samuel  P.  Langley  was  probably  the  first  airplane  designer  affected.  His 
misfortune  occurred  just  prior  to  the  Wrights'  first  flight.  During  the 
launch  of  Langley's  monoplane  off  the  Potomac  River  houseboat,  catastrophic 
wing  divergence  occurred.  His  failure  with  the  nonowing  and  the  Wrights' 
success  with  the  biplane,  caibined  with  the  lack  of  a  torsional  stiffness 
criterion  for  the  moncwing,  resulted  in  the  favoring  of  the  dual  wing  design 
during  the  early  period.  Ihe  p)articilitY  was  understandable  since  no  designer 
vanted  to  over-stiffen  the  wing  at  the  escpense  of  added  weight. 

Althou^  a  few  externally  braced  monoplanes  were  built  prior  to  World  War 
I,  military  monoplane  design  virtually  stc^yped  between  1917  and  the 
mid- thirties.  During  the  biplane  era,  the  most  ccnmon  aeroelastic  problem  was 
tail  flutter.  One  of  the  first  documented  cases  occurred  on  the  tail  of  the 
British  twin-engined  Handley  Page  0/400  bcnber  at  the  start  of  the  first  war. 
The  noted  aerodynamicists  Lanchester  and  Bairstow  investigated  the  reasons  for 
the  flutter  incidents  and  the  following  is  a  quote  frcm  Lanchester 's  report; 

"The  difficulty  ej^jerienced  is  that  at  certain  critical  speeds  of  flight 
a  tail  wobble  is  set  up,  involving  heavy  torsional  stresses  on  the 
fuselage,  the  l^pe  of  vibration  being  an  angular  oscillation 
approximately  about  the  axis  of  the  fuselage;  I  am  informed  that  the 
angular  magni'tade  of  this  oscillation  amounts  at  times  to  something 
approaching  15  ,  and  is  undo;±>tedly  extremely  dangerous  to  the  structure 
of  the  machine.  I  gather  that  the  experience  of  the  pilots  vdien  this 
vibration  is  at  its  worst  is  terrifying." 

The  problem  ves  being  caused  by  the  coupling  between  the  fuselage  torsion 
mode  and  the  anti-symmetrical  elevator  excitation  mode.  In  the  latter,  the 
left  and  ri^rt  elevators  oscillated  about  their  hinge  lines  in  oppxDsite  phase. 
There  was  no  interconnecting  torque  tube  to  prevent  the  occurrence.  The  other 
mode  was  merely  the  low  frequency  torsional  oscillation  of  the  aircraft 
fuselage.  Since  the  vibrating  frequencies  of  these  two  separate  oscillatory 
modes  were  almost  the  same,  a  resonant  coupling  occurred,  with  the  resulting 
tail  wobble  phenomenon.  A  same  typ)e  of  tail  flutter  problCTi  was  e5^)erienced 
by  the  DH-9  in  1917,  resulting  in  several  fatalities.  The  fix  used  then  was 
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the  same  one  used  for  the  Handley  Page  botiber.  It  has  been  a  design  feature 
ever  since  for  reversible  flight  control  systems  -  the  left  and  right 
elevators  were  connected  with  a  stiff  torque  tube. 

Wing  problems  appeared  with  the  return  of  the  monowing.  Insufficient 
torsional  rigidity  -  led  to  divergence,  loss  of  aileron  effectiveness,  and 
flutter.  An  early  problem  occurred  on  the  Fokker  D-8.  Ihe  aircraft  was  built 
for  WWI  and  due  to  superior  performance  was  immediately  placed  in  coribat. 
Within  a  few  days  several  wing  failxnres  occurred  during  high  speed  dives.  The 
wing  torsional  stiffness  criterion  used  in  the  initial  design  was  the  same  as 
that  previously  used  for  biplanes. 

Ihe  Anty  ran  static  strength  tests  and  discovered  that  the  wings  were 
more  than  capable  of  withstanding  the  6g  design  limit.  Why  then,  the 
failures?  Confronted  with  the  dilenma,  Fokker  decided  to  conduct  his  own 
static  tests.  He  found  that  althou^  the  wing  did  indeed  have  the  design 
strength,  under  increasing  loads  the  angle  of  incidence  at  the  wing  tips 
increased,  relative  to  the  roots:  the  necessary  condition  for  divergence. 
In  the  high  speed  dives  the  air  loads  increased  faster  at  the  tips  and  the 
resulting  torsion  caused  the  wings  to  fail. 

After  the  war,  the  U.S.  Amy  encountered  a  violent  but  nondestructive 
case  of  wing-bending/aileron -rotation  flutter  on  the  same  Fokker  D-8.  The 
cure  for  the  problem,  as  with  virtually  any  control  surface  type  flutter,  was 
mass  balancing  of  the  control  surface. 

The  period  of  monoplane  development,  because  of  the  accotpanying 
resurgence  of  aeroelastic  problems,  initiated  the  first  serious  research  in 
the  field  of  aeroelasticity.  Early  day  techniques  were  mainly  cut-and-try. 
Many  analytic  theories  were  presented  on  wing  load  distribution,  wing 
divergence,  loss  of  lateral  control  and  aileron  reversal.  Potential  flow 
flutter  analysis  was  sufficiently  landerstood  by  1935  to  be  incorporated  into 
aircraft  design  but  the  majority  of  airplane  designers  were  reluctant  to  trust 
the  mathematicians  to  formulate  criteria  for  the  strength  and  rigidity  of 
aircraft  structural  oatponents. 

From  1934  to  1937,  the  perennial  arms  race  resulted  in  the  development  of 
many  new  types  of  aircraft.  Nimerous  cases  of  flutter  appeared,  mainly  of  the 
tail  type.  The  accidents  served  to  underscore  the  critical  nature  of 
mass  distribution  and  control  surface  mass  balancing.  In  1938,  a  panel  of 
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scientists  boarded  a  four-engined  Jvmker  aircraft  for  an  in-flight  observation 
of  a  planned  flutter  test.  All  perished  when  a  catastropMc  condition 
occurred.  The  accident  served  due  warning  to  the  coiinunity  of  the 
difficulties  and  hazards  of  flight  flutter  testing. 

With  the  development  of  irtproved  excitation  and  measuring  equipment,  and 
a  better  theoretical  understanding  of  the  nature  of  the  flutter  problem, 
flutter  testing  came  of  age  diaring  the  late  forties.  Analytic  predictions  are 
still  inprecise,  but  with  improved  model  and  wind  tunnel  technology,  and  with 
the  tranendous  surge  in  the  development  of  the  high  speed  ccitputers, 
investigation  of  aeroelastic  phenomena  is  beccming  more  and  more  a  controlled 
science. 


12.5.3  Mathematical  Analysis 

Hie  mathematical  treatment  of  aeroelastic  pjhencmena  is  sones^at  complex 
and  tedious.  It  requires  a  knowledge  of  the  structural  properties  of  the 
vehicle  and  of  the  nature  of  the  aerodynamic  forces  acting  on  the  surface  of 
the  airplane.  Many  textbooks  and  other  publications  exist  vdiich  treat  the 
subject  in  thorou^  detail  (References  12.1,  12.2,  12.10,  and  12.11).  This 
chapter  will  attenpt  to  only  touch  upon  the  essential  and  basic  ideas  behind 
the  mathematical  formulation. 

Specifically,  we  will  look  at  wing  torsional  divergence,  aileron 
reversal,  and  flutter.  The  first  two  are  static  phenomena  and  are  not 
significant  present-day  problems.  Vfe  treat  them  because  they  are  easy  to 
analyze  and  visualize.  Flutter  is  always  a  primary  design  consideration  and 
wLll  therefore  be  treated  in  detail. 

12.5.4  Wing  Torsional  Divergence 

If  a  wing  in  steady  flight  is  slightly  perturbed,  an  aerodynamic  mcment 
will  be  induced  which  tends  to  twist  the  wing  further.  Since  the  structure's 
stiffeess  is  independent  of  speed  of  flight,  and  the  aerodynamic  mcment  is 
proportional  to  the  flight  velocity  squared,  there  may  exist  a  critical  speed 
at  which  the  elastic  stiffness  of  the  wing  is  barely  able  to  sustain  the  v/ing 
in  its  deformed  state.  This  speed  is  called  the  divergence  speed  and  the  wing 
is  said  to  be  torsionally  divergent. 

More  simply  stated,  wing  divergence  is  the  condition  vAiere  the 
aerodynamic  mcment  exceeds  the  elastic  restoring  mcment  in  torsion.  Figure 
12.73  ill\istrates  a  typical  canSoered  wing  section  moving  with  a  flight  speed 

U. 
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a  - 

a,  =  INITIAL  ANGLE  OF  ATTACK  (RIGID) 
e  =  ANGLE  OP  TWIST 


FIGURE  12.73.  CAMBERED  WING  SECTION 

Acting  at  the  aerot^amic  center  is  the  lift  (L)  and  the  moment  about  the 
a.c.  elastic  restoring  mcment  on  the  wing  opposing  the  total 
aerod^Tiairdc  mcment  is  depicted  by  a  linear  coiled  spring  with  spring  constant 
Kg  attached  at  the  elastic  axis  (E.A. )  located  a  distance  e  behind  the  a.c. 
The  elastic  axis  is  defined  as  the  axis  about  which  the  wing  section  would 
twist,  subjected  to  pure  mcment.  A  vertical  force  applied  at  the  E.A.  would 
effect  only  a  vertical  deflection  (bending) .  No  twisting  would  result.  Rcr 
the  sake  of  definition,  the  bending  and  twisting  of  a  wing  are  illustrated 
separately  in  Figure  12.74. 
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FIGURE  12.74. 


TONG  BENDING  AND  IVJISTING 


“th©  wing  section  in  Figure  12.73,  the  angle  of  attack  on  the  airfoil 
is  ccfiposed  ox  tv70  parts:  c;^,  the  angle  of  attack  that  a  rigid  wing  would  see, 
and  6,  the  angle  resulting  fron  elastic  twist.  At  equilibrium,  the 
aerod^^namic  mctnent  about  the  elastic  axis  must  equal  the  elastic  restoring 
mcment  K-6. 


Aerodinnamic  moment  about  the  elastic  axis 


m 


AERO 


=  Le 


.m 


ac  =  ^6® 


Fran  subsonic  aeroc^amics,  we  know  that 


L  =  qSC^  =qSCL(a  +e)  and 
a  .  a 

"J  S  o 
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Equating  the  mcments, 


=  qS 


®  Cl  '“r  +  8>  -  =  ISlac’ 


Solving  for  the  twist  angle  6, 


6  = 


S  (®  \  “r  '  =  Stec) 


"e 


(12.30) 


At  the  divergence  condition  the  twist  angle  grows  without  bound  (i.e.  . 
This  condition  occurs  inathematically  vAien  the  denoninator  of  Equation  12.30 
equals  zero. 


Kq  '  ^SeC^^  =  0 

vdiere  q  =  q^^  (Dynamic  Pressure  at  Divergence) 

„  _  P  n  2  _  _ 1 

2  Se  Cj^ 


or,  solving  for  the  divergence  speed  Uj^, 


(12.31) 


The  design  parameters  affecting  the  divergence  of  straight  wings  are 
primarily  the  wing  torsional  stiffness  (Kg)  and  the  offset  distance  e. 
Increasing  Kg  is  a  costly  process  at  the  expense  of  considerable  weight.  An 
approach  more  frequently  enployed  is  to  proportion  the  wing  structurally  so  as 
to  move  the  elastic  axis  forward.  Aft  wing  sweep  lessens  the  divergence 
problem  since  swept  wing  tip  angles  of  attack  are  effectively  reduced  by  wing 


bending. 

If  a  real,  three  dimensional,  wing  were  considered,  the  actual  divergence 
speed  must  be  obtained  by  integration  across  the  span  of  the  wing  since  local 
spanwise  properties  will  vary. 
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12.5.5  Aileron  Reversal 


The  aeroelastic  prcblem  of  aileron  reversal  is  closely  related  to  wing 
torsional  divergence  in  that  both  depend  strongly  upon  the  torsional  stiffiiess 
of  the  wing.  The  history  of  this  subject  closely  parallels  that  of  wing 
torsional  divergence;  however,  it  was  during  World  War  II  that  the  problem 
came  into  iirportance.  The  increasing  speeds  and  the  requirements  for  rolling 
performance  precipitated  the  problem  greatly.  Following  WW  II,  and  with  the 
advent  of  thin-wings,  moderately  high  aspect  ratio,  and  sweepback,  the  aileron 
reversal  speed  became  of  prime  irtportance. 

Aircraft  with  conventional  planforms  may  suffer  serious  loss  of  aileron, 
elevator,  and  rudder  control  effectiveness  due  to  elastic  deformations  of  the 
structure.  The  aileron  controls  the  rolling  motion  of  an  aircraft,  and  when 
placed  downward,  the  lift  over  the  wing  is  increased  and  a  rolling  mcment  is 
produced.  The  down  aileron  also  produces  a  twisting  mcment  on  the  wing  vhich 
tends  to  twist  the  wing  nose  down  and  reduce  the  angle  of  attack  and  hence 
this  reduces  the  rolling  itonent.  A  similar  situation  exists  for  the  up 
aileron  wing.  The  up  aileron  produces  a  torque  vdiich  increases  the  angle  of 
attack  and  decreases  the  effective  rolling  inanent.  Since  the  wing  stiffness 
is  independent  of  the  fli^t  velocity  and  the  aerodynamic  force  ineffective  in 
producing  a  rolling  mcment  and  the  resulting  wing  segment  twist  and  aileron 
deflection  produce  no  effective  change  in  lift.  Beyond  this  critical  speed, 
the  effect  of  ailerons  is  actually  reversed.  This  analysis  also  applies  to 
the  other  control  surface  and  is  sometimes  known  as  "Control  Surface 
Effectiveness."  Figure  12.75  shews  how  aileron  effectiveness  as  measured  by 
the  ratio  of  rolling  velocity  to  aileron  angle  is  affected  by  forward  speed 
for  a  WW  II  filter-type  aircraft  at  sea  level. 

Avoidance  of  aileron  reversal  in  a  straight  wing  with  conventional 
ailerons  is  a  matter  of  providing  sufficient  wing  torsional  stiffness.  If  the 
wing  is  swept  back  aileron  reversal  is  a  serious  problem  and  wing  bending 
stiffness  must  also  be  increased.  This  semetimes  becomes  prohibitively  large 
in  a  weight  analysis,  hence,  the  other  means  are  eitployed.  Spoilers,  fully 
moving  wing  tips,  and  even  moving  the  aileron  inboard  all  produce  suitable 
means  to  corrbat  this  occurrence.  Figure  12.76  shows  that  the  aileron  reversal 
speed  can  be  increased  by  changing  the  configuration  of  the  aileron  controls. 
Mditionally,  this  figure  shews  that  control  effectiveness  can  be  increased 
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vd-th  decreasing  sweepback.  Inboard  ailerons  are  used  in  the  F-4,  vhile  the 
T-38  has  ailerons  at  the  send-span  position.  Often  aircraft  are  fitted  with 
two  sets  of  aileron  surfaces,  the  outboard  being  locked  out  at  high  speeds  and 
only  the  inboard  ailerons  being  used  for  roll  control.  The  cases  of  elevator 
and  rudder  control  effectiwness  and  reversal  are  usually  less  critical  than 
those  of  aileron,  they  are,  however,  more  ccrtplicated  due  to  the  large 
deformations  of  the  fuselage  and  attachments. 


AILERON 

REVERSAL 

SPEED 


FIGURE  12.75.  AILERON  EFFECTIVENESS 


12.125 


It  should  be  noted  that  flutter  and  divergence  correspond  to  conditions 
of  aeroelastic  instability,  and  that  speeds  beyond  critical  flutter  and 
divergence  speeds  will  result  in  an  eventual  structxnral  failure.  Control 
reversal  is  not  a  condition  of  instability,  and  speeds  beyond  the  control 
reversal  speed  will  result  only  in  a  reversal  of  the  action  of  the  control 
system  and  not  necessarily  a  failure  of  the  structure.  Hie  priire  concern  here 
is  the  loss  of  maneuverability.  With  the  advent  of  fully  boosted  systems, 
additional  cotplications  result  because  of  the  deformations  resulting  in  the 
controlling  mechanism  and  the  enormous  mechanical  advantage  available. 


DISTANCE  FROM  t  AIRPLANE  TO  AILERON 


FIGURE  12.76.  AILE30I  REVERSAL  SPEED  VS  AHJERCXJ  POSITION 
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Figure  12.77  illustrates  a  typical  cambered  vdng  section  with  ailerons. 
0  represents  a  downward  deflection  of  the  aileron.  If  the  wing  were  rigid, 
the  aileron  displacement  woiald  be  acconpanied  by  an  increase  in  lift. 
However,  for  an  elastic  wing  the  deflected  aileron  would  also  cause  a 
nose-down  twist  of  the  wing,  reducing  the  effective  angle  of  attack.  Ihe  nose 
down  twisting  nonents  increase  with  the  square  of  airspeed  whereas  the  elastic 
restoring  moment  stays  constant.  At  seme  critical  airspeed  Uj^,  called  the 
aileron  reversal  speed,  the  increase  in  lift  caused  by  the  deflected  aileron 
is  completely  negated  by  the  loss  in  lift  due  to  the  reduction  of  the 
effective  angle  of  attack  caused  by  the  twist.  The  mathematics  follow. 
Again,  the  elastic  restoring  moment  is  KqS. 


Aero  dynamic  moment  about  the  elastic  axis 
•  ^  ->^ac  - 


a 


where  m .  is  the  contribution  to  the  moment  about  the  a.c.  due  to  aileron 
ac^ 

deflection. 


L  = 


=  q  S  Cj^  (a^  -  6)  +  q  S 


or 


L  =  q  S 


'“r  -  *a] 


Substitution  into  the  moment  equation  yields 


a  J 


'^AEFO  =  g  S®  '“r  *‘=L,  ' 

-  q  Sc 

Again,  equating  elastic  and  aerodynamic  moments. 


[\ 

^lac  Slg  ^aj 


K  0  =  W.  =  Le  -771  -  771 

6  AEPO  ac  ac  6, 


=  q  Se  ■  ^3  SC 

a 
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Solving  for  the  twist  angle  9, 


e 


q  S 


(• 


Kq  +  q  Se 

a 


6 

a 


q  s 
+  — 


Kq  +  q  se 


(12.32) 


Ultimately,  we  will  be  interested  in  only  the  change  in  twist  angle  due  to  a 
change  in  aileron  deflection  i.e., 96/36^.  Consequently,  the  second  term  in 
the  6  equation  drops  out  since  it  remains  invariant  with  a  change  in  aileron 
deflection.  As  a  result. 


FIGURE  12.77.  CAMBERED  T-ffNG  SECTICN  PTiTH  AILERON 
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K„ 


+  q  Se  Cj^ 
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At  Uj^,  the  change  in  lift  due  to  a  change  in  aileron  deflection, 

AL^  ,  beccmes  zero, 
a 

AL^  =  0 

a 
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The  change  in  lift  AL^  is  r^resented  by 

a 

ALj  =  Lg  +  L^q,  where 
a  a 

Lgg  =  lift  increase  due  to  the  increase  in  caitiber  =  g  S  Cj^  6^ 

6 

a 

L^q  =  lift  decrease  due  to  the  leading  edge  twist  =  -  g  S  ||- 

a  a 

Therefore, 


.  =  q  S  ,  S  If 


a  a 


Substituting  Eguation  12.33  for  36/36  yields 

a 


\ _ a _ a 

Kq  +  qse  Cj. 


Rearranging  and  cancelling,  we  obtain 


qS^KeC^^  +qScC^^(^,^ 
Ke  +  q  se  ^ 


(12.34) 


At  the  reversal  speed  AL.  —  0.  Since  q,  S,  and  6  in  the  numerator  of 

a  ^ 

Equation  12.34  are  all  non-zero,  the  terms  within  the  parentheses  must  vanish 
to  satisfy  the  equation 


'e%  V  =0 

‘a  “ 


where  q  =  (Eynamic  Pressxure  at  Control  Reversal) 


®  Si/ 

a  0 
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notevrortl^  observation  is  the  fact  that  the  reversal  speed  is  independent  of 
twist  axis  location  (e) .  The  aerodynamic  mcment  at  the  reversal  condition  is 
a  pure  couple  and  therefore  independent  of  axis  position.  Also  notice  the  Up^ 
decreases  with  a  decrease  in  altitude. 

Preventing  aileron  reversal  in  a  straight  wing  is  a  matter^ of  increasing 
torsional  stiffness  (K„) ,  increasing  the  aileron  effectiveness 


or  of 


decreasing  the  magnitude  of  .  In  the  case  of  si'/ept  wings,  where  aileron 


a 

reversal  has  been  a  serious  problem,  bending  stiffness  must  also  be  increased. 
Since  veight  increases  acccnpany  stiffness  increases  other  methods  should  be 
€n5:loyed.  Alternative  methods  such  as  spoilers  ard  all  moving  wing  tips  have 
proved  beneficial.  Also,  the  effective  Kg  is  increased  as  the  aileron 
locations  are  moved  inboard. 


12.5.6  Flutter 

First  we  assume  a  cantilever  straight  wing  without  ailerons  mounted  in  a 
wind  tionnel  and  with  no  airflow.  I*7hen  the  model  is  disturbed,  oscillations 
are  induced  v^ich  gradually  danp  out  since  the  elastic  structure  provides  a 
daitping  ratio  of  its  own,  known  as  structural  damping,  usually  from  2%  to  8%. 
As  the  speed  of  the  wind  flow  is  increased,  the  rate  of  darrping  of  the 
oscillations  increases  due  to  aerodynamic  damping.  With  further  increases  of 
speed,  a  point  is  reached  vdiere  the  daitping  starts  to  decrease  rapidly.  At 
the  next  point,  the  Critical  Flutter  Speed,  the  oscillation  can  just  maintain 
itself  with  steady  amplitude.  Speeds  above  the  critical  flutter  speed  trigger 
a  violent  oscillation  and  subsequent  destruction  of  the  section  when  subjected 
to  a  small  disturbance.  This  airfoil  is  said  to  have  oscillatory  instability 
and  is  said  to  flutter.  Moreover,  once  the  oscillation  starts,  it  is  self 
sustained  and  no  further  external  forces  or  forcing  functions  are  required. 
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Mditional  types  of  flutter  can  involve  aileron  motion  v^ere  there  may  be 
one  or  more  ranges  of  speeds  for  vMch  flutter  occurs.  The  aileron  flutter 
induces  the  wing  to  flutter.  Usually  these  regimes  are  bounded  at  both  ends 
by  critical  speeds  where  one  has  an  oscillation  of  constant  amplitude. 

The  oscillatory  motion  of  a  fluttering  cantilever  wing  has  both  bending 
and  torsional  ccnponents.  If  the  airfoil  is  rigid  in  torsion  and  is 
constrained  to  have  only  a  flexural  degree  of  freedom,  it  will  not  flutter. 
With  only  torsional  degrees  of  freedcm  it  can  flutter  only  if  the  angle  of 
attack  is  near  the  stall  angle.  Thus,  coupling  of  several  degrees  of  freedcm 
is  a  necessary  portion  of  flutter.  Furthermore,  bending  movements  at  all 
points  across  the  span  are  approximately  in  phase  with  one  another.  The 
torsional  movements  are  all  approximately  in  phase;  however,  the  bending  mode 
is  lasually  considerably  out  of  jhase  with  the  torsional  movement.  This  jhase 
difference  is  responsible  for  flutter. 

An  airplane  wing  is  an  elastic  body  and  has  infinitely  many  degrees  of 
freedcm  in  vibration.  The  basic  construction  allows  any  elastic  deformation 
in  the  chordwise  section  to  be  described  by  (1)  deflection  of  a  reference 
point,  (2)  and  angle  of  rotation  about  that  point.  With  control  surfaces,  the 
fi^ssdcm  to  turn  about  the  hinge  line  is  so  much  more  iitportant  than  elastic 
deformation  that  this  deflection  is  best  described  by  an  angle  of  rotation 
about  its  hinge  line. 

One  must  consider  three  variables  in  wing  flutter  (this  does  not  include 
a  time  variable) : 

(1)  Bending 

(2)  Torsion 

(3)  Control  surface  rotation 

A  flutter  mode  vhich  contains  all  three  is  termed  ternary  flutter.  A  flutter 
mode  which  contains  only  two  (usually  the  first  two)  is  called  binary  flutter. 
Occasionally,  a  single  degree  of  freedom  oscillation  may  exist.  This  usually 
occurs  with  control  surfaces  and  is  referred  to  as  'buzz." 

In  general  the  design  criteria  requires  that  an  aircraft  must  be  able  to 
fly  near  without  the  appearance  of  undesirable  marginal  stability  in  its 
structural  vibrations.  Yet,  bending  and  torsional  flutter  can  emerge  suddenly 
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and  violently,  and  at  about  five  knots  above  the  wing  will  possibly  destroy 
itself  after  two  or  three  cycles  of  oscillations. 

Ihe  classical  type  of  flutter  nearly  always  involves  the  coupling  of  two 
or  more  degrees  of  freedom.  Ihe  analysis  can  be  very  cotplex,  depending  vpon 
the  desired  detail.  We  shall  attatrpt  only  an  elementary  treatment  here,  with 
etifiiasis  on  formulation  and  physical  interpretation.  Ihe  basic  second  order 
system  will  be  studied  in  depth  because  it  conveniently  demonstrates  the 
stability  problem,  the  basis  of  all  flutter  phenomena.  Vfe  will  be  mainly 
concerned  with  the  changes  in  systan  response  resulting  from  changes  in 
danping  and/or  excitation  frequency. 

Figure  12.78  shows  a  typical  wing  section  with  two  degrees  of  freedom; 
X,  a  vertical  displacement  downward  (bending),  and  a,  the  angular  deflection 
about  the  elastic  axis.  For  sinplification,  the  E.A.  and  the  c.g.  are 
coincident. 

Unlike  the  static  divergence  and  control  reversal  pr<±»lems,  the  flutter 
analysis  will  require  knowledge  of  the  mass  and  inertial  properties  of  the 
wing  since  an  oscillatory  motion  will  be  involved,  resulting  in  a  "generation" 
of  inertial  forces.  By  use  of  Newton's  second  law,  the  equations  of  motion 
for  the  wing  in  the  two  degrees  of  freedom  can  be  derived. 
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FIGURE  12.78.  TWO  DEGREE  OF  FREEDCFi  I'TING  SECTION 

NOTE:  f  (t)  =  time  varying  force  acting  at  c.g. 

M(t)  =  tome  varying  mcment  acting  on  wing  section 
m  =  mass  of  section 
I  =  mcment  of  inertia  about  c.g. 

k 

'x'  Ci  =  spring  constants  of  displacement  (bending) 
and  torsion,  respectively 

d,  b  =  bending  and  torsion  dairping  coefficients 
m  X  +  d  X  +  k^x  =  f(t) 

I  a  +  b  a  +  k^a  =  M(t) 

As  the  eguations  stand,  for  small  displacsnents  there  is  neither  inertial 
nor  elastic  coupling  between  the  two  degrees  of  freedcm.  That  is,  in  the 
absence  of  the  external  driving  forces  f{t)  and  M(t),  the  free  body  motion  is 
independent  in  both  x  and  a.  The  motions  do  not  interfere.  This  is  obvious 
fron  inflection  (no  a  in  the  x  equation,  and  vice  versa)  and  is  the  result  of 
our  collocating  the  c.g.,  the  axis  of  twist,  and  the  displacement  axis.  Ihe 
only  way  coupling  can  be  introduced  is  through  the  driving  forces  f(t)  and 
M(t) . 


12.133 


Before  shewing  hew  the  equations  can  be  coupled  aerodynamically,  let's 
first  look  at  the  solutions  to  the  equations.  As  they  are  presented,  the 
equations  are  linear  and  can  be  solved  separately.  The  solution  for  x,  for 
exairple,  is  the  sum  of  a  transient  (honogeneous)  and  a  steaefy-state  solution. 


x(t)  =  X^(t)  +  Xg{t) 
Placing  the  x  equation  in  standard  form 

X  +  2^(d^x  +  to^^x  =  f ' (t) 

2  ^x 

vdiere  u  ^  =  -ii 

n  m 


*1116  undairped  natural  frequency  is  the  frequency  the  system  would 

oscillate  in  the  absence  of  dairping.  Ihe  damping  ratio  c,  is  sinply  a 
convenient  way  of  representing  the  degree  of  damping  (due  to  the  dash-pot 
coefficient,  d)  relative  to  the  size  of  the  mass  and  spring  constants.  Ihe 
method  for  finding  the  transient  solution  can  be  found  in  numerous  texts  on 


differential  equations. 


x^  =  e 

L. 


(C^  cos  0)^  t  +  C2  sin  t) 


(12.38) 


0).  =  u 
d  n 


V  1  -  o  <  ^  < 


vtere  the  constant  is  the  damped  frequency,  and  for  obvious  reasons  is  less 
than  the  undamped  frequency  We  see  that  the  transient  response  has  the 
form  of  a  sinusoidal  oscillation  vhose  envelope  decreases  ejponentially  with 
time. 

To  solve  for  the  steaf^-state  solution  x  ,  knowledge  of  the  driving  force 
f(t)  is  required.  Flutter  has  been  observed  to  be  sinusoidal,  therefore  we 
c:an  assume  that  the  nature  of  the  driving  force  is  also  sini:isoidal. 


NOTE;  *In  this  example,  is  the  undamped  translational  structural 

frequency. 
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f(t)  =  F  cos  ut 


lie  can  then  assune  a  steady-state  solution  of  the  form 


Xg  =  A  cos  (OJ^  +  (j)) 


Substitution  into  the  differential  equation  yields  values  for  the  constants. 
Chapter  3  of  Reference  11  gives  the  details  to  the  solution. 


and 


The  steady-state  response  vdth  amplitude  A  and  phase  angle  <f,  is  of  the  same 
form  as  tlie  driving  force.  The  amplification  factor  A/  (F/kj^)  is  the  ratio  of 
the  anplitude  A  of  the  steady-state  solution  to  the  static  deflection  F/k  due 
to  a  constant  force  of  magnitude  F.  Figure  12.79  shows  a  plot  of  the 
anplification  factor  versus  frequency  for  different  values  of 
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FREQUENCY  RATIOj^. 

COit 


FIGUPE  12.79.  AT-IPLIFia^TICW  FACTOR  AND  PHASE  ANGLE  VERSUS  FREQUENCY  RATIO 

The  airplification  factor  is  ijnity  the  input  frequency  is  zero, 

regardless  of  the  daitping  ratio  ?.  It  is  important  to  note  that  when  the 
forcing  frequency  (co)  equals  the  vmdaitped  natural  frequency  (u  ) ,  the 
amplification  factor  is  1/(2  ?) .  Consequently,  any  amount  of  darrping  is  very 
critical  in  holding  down  the  response.  Lack  of  damping  results  in  an 
unbounded  output. 

Note  also  that  at  the  critical  frequency  (oj  =  u  )  the  phase  angle  st  is 
90^,  regardless  of  It  can  be  shown  that  during  this  condition,  the 

response  (bending)  velocity  is  exactly  in  phase  with  the  driving  force 
f  (t) ,  vhich  is  the  condition  for  maxirnom  pover  input  into  the  oscillating  wing 
system  fron  the  external  force. 

How  do  we  relate  the  above  analysis  to  a  meaningful  wing  flutter  problem? 

Looking  back  at  Figure  12.78  we  can  take  the  x  degree  of  freedom  as  the  basic 

wing  bending  mode  with  natural  frequency  of  oscillation  w  =  /k  /m,  a  value 

n  X 

obvio\3sly  dependent  upon  the  mass  and  elastic  properties  of  the  wing.  The 
external  force  f(t)  can  be  inferred  to  be  the  aerodynamic  lift  acting  on  the 
given  wing  section. 
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f{t)  =  -  L(t)  =  -q  S  a(t) 

a 

If  a(t)  were  to  vary  sinusoidally  at  frequency  u,  we'd  have 

f  (t)  =  -  q  S  Oq  cos  lot  (ttp  is  the  magnitude  of  the  angle  of 
a  attack  variations  with  time) 

or. 


f(t)  =  F  cos  tot 


This  is  exactly  the  same  type  of  external  driving  force  described  in  the 
exanple  above.  Sixrh  a  variation  in  a  can  also  result  from  the  rotational 
equation  of  motion. 


la  +  ba  +  k  a  =  M(t) 
a 

Solving  for  the  free  body  response  [M  (t)  =  0] ,  the  solution  is  similar  to  the 
displacement  transient  solution  (Equation  12.38) 


a(t)  =  e’^T^a"*^ 


Cj  cos 


In  most  wLng  structures,  bending  daitping  is  much  greater  than  torsional 
dairping.  Therefore,  without  loss  of  generality'  v;e  can  assume  torsional 
structural  danping  (s^)  to  be  small.  The  solution  then  has  the  appraximate 
form 


a(t)  =  C,  cos  u  t  +  C_  sin  oj  t, 

1  a  2  a  ' 

'^tere  0)^  =  -yj k^/I  represents  the  natural  undaitped  torsional  frequency  of  the 
wing.  If  the  initial  conditions  are  picked  appropriately,  a(t)  can  be 
represented  by 

a(t)  =  cos  03^t  (12.39) 
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The  critical  condition  now  occurs  when  o  =  o  =  u  i.e.,  vAien  the  natural 

an 

frequencies  of  wing  bending  and  wing  torsion  are  the  same.  Substituting 
0)  =  (0^  in  Equation  12.39  yields 


or. 


a(t)  =  Oq  cos  lot 


cos  tot  = 


Substituting  this  into  the  forcing  function  for  displacement  f(t)  gives 


f(t)  =|-a{t) 

0 

If  F/Oq  is  set  equal  to  some  constant  A,  then 

f  (t)  =  Aa(t) 

Equation  12.36  then  becomes 

mx  +  dx  +  kx  =  Aa(t) 

In  this  exaitple  the  result  is  aerodynamic  coupling  between  the  tvro  degrees  of 
freedom,  bending  and  torsion.  Under  certain  conditions  of  danping,  the  end 
product  from  the  coupling  will  be  classical  bending-torsional  flutter. 

A  side  view  of  the  flight  path  of  the  wing  in  bending-torsional  flutter 
is  shown  in  Figure  12.80. 
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WING  FLIGHT  PATH 


a  =  o 


FIGURE  12.80.  BENDING-TORSICmL  FLUTTER  (oj  =  to  =  to  ) 

a  Ti 

This  is  a  very  sirtplified  explanation  of  wing  bending-torsion  flutter. 
■Hie  bending  mode  is  coupled  aerodynamically  with  the  torsion  mode  to  effect 
the  dynamic  instability.  The  physical  interpretation  is  rather 
straightforward  but  precise  predictions  on  actual  wing  flutter  speeds  and 
frequencies  are  extremely  difficult  to  obtain  analytically,  in  our  exanple, 
several  specific  assurtptions  and  sinplifications  were  made:  (1)  the  elastic 
a>:is  was  placed  at  the  c.g.,  (2)  the  wing  section  was  treated  as  a  rigid  one, 
and  (3)  the  aerodynamic  dartping  and  air  mass  accelerations  were  ignored. 
Hcwever,  these  assunptions  did  not  detract  from  the  general  ideas  basic  to  the 
dynarriic  instability'. 

Other  types  of  flutter  vhich  may  occur  usually  involve  the  rotation  of 
seme  control  sinrface  about  its  hinge  line,  wing  bending/aileron  rotation 
flutter,  or  fuselage  torsion/rudder  rotation  flutter  are  exaitples.  It  is  also 
possible,  in  fact  more  ccainon  for  more  than  two  modes  to  be  present,  such  as 
in  wing  bending/wing  torsion/aileron  rotation  flutter.  Theodore  Theodorsen 
presented  a  mCk  paper  in  1935  which  treats  the  three^egree-of-freedom 
wing— aileron  airfoil  in  semewhat  ccnplete  and  tedious  detail. 

12.5.7  Structural  Modeling 

To  expand  the  analysis  to  treat  the  total  aircraft  structure  and  to  cover 
all  interactions  between  structural  corrponents  is  considerably  more  difficult. 
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The  problem  is  basically  twofold  in  the  derivation  of  the  equations  of  motion. 
First,  the  free  body  mass  and  structural  relationships  have  to  be  determined. 
Second,  the  aeroc^amic  forces  acting  on  the  entire  surface  have  to  be 
derived.  These  forces  will  in  turn  depend  upon  the  deflections,  velocities, 
and  accelerations  of  the  structure. 

General  equations  can  be  formulated  by  treating  the  aircraft  as  catposed 
of  a  large  number  of  discrete  masses.  Rewriting  Equation  12.36  in  matrix  form 
yields  the  general  equation 

[m]  Cx}  +  [d]  {x}  +  [k]  X  =  f(t) 

The  brackets  []  signiify  a  square  matrix  and  the  braces  indicate  a  column 
matrix.  Represented  are  n  simultaneous  equations  in  degrees  of  freedon  x^^* 
(i  =  1,  2,  ...,  n) .  The  equations  are  difficult  to  formulate  and  difficult  to 
solve.  The  main  problem  lies  in  determining  the  stiffness  coefficients  K.  .. 

To  obtain  the  harmonic  solution  to  the  above  equations  a  continuous 
forcing  function,  F^,  is  applied  equal  to  the  structural  danping  d^j.  As  a 
result,  the  forcing  function  and  structural  damping  cancel  and  the  equations 
are  simplified  to 


[m]  {x3  +  [k]  £x3  =  0 

Next,  a  matrix  of  flexibility  influence  coefficients  is  defined  as  the 
inverse  of  the  stiffness  matrix 

[c]  =  [k]"^ 

The  equations  now  beccme 

[c]  [m]  £x3  +  Cx3  =  0 

The  influence  coefficients  are  preferred  to  the  stiffness  coefficient  beca\ise 
they  are  more  conveniently  determined  by  actual  measurements  on  the  structure 
or  a  scaled  model.  Using  the  simple  spring  relationship  shown  in  Figure 
12.81,  the  force  F  required  to  deflect  a  spring  (or  wing  section  equivalent)  a 
distance  x  in  static  equilibrium  is  governed  by 
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F  =  kx 

Or ,  t^ing  a  different  view,  the  displacement  x  resulting  from  the  application 
of  a  force  F  is  given  by 


X  =  1/k  F  =  C  F. 


*lhe  coordinate  x^  can  represent  either  a  displacement  or  an  angular 
deflection. 


F*=kx 


I  I  K  u 

F  r  »n  — ^  ^ 


t 

F 


FIGURE  12.81.  SIMPLE  SPRING  REEATICNSHEP 

If  we  assime,  for  exaitple,  that  a  wing  in  Figure  12.82  is  approximated  by  five 
separate  sections,  the  displacement  at  section  3  due  to  forces  F2  and  F^ 
applied  at  stations  two  and  four,  respectively,  is  given  by 

^3  ^  ^32  ^2  ^4  ^4 
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FIGUBE  12.82.  CArJTIIjE'TER  T-JING 

In  general,  a  displacement  at  Station  (or  mass  point)  i  due  to  the 
existence  of  forces  at  Station  j  is  given  by  the  following  relationship: 


Xi  -  ‘^i?  ^2  ‘  ^in  ^n 

n 

X.  =  E  C.  .  F. 

^  j=l  => 

Conversely,  a  force  at  Station  i  will  cause  deflections  at  Stations  j, 
governed  by 

Fi  =  k. ^  X3_  +  k . 2  X2  +  .  .  .  +  k^  Xj^ 

n 


Ihe  C^j  can  be  determined 


X. 

■,  _  ^ 
'ij  F. 
—  ■  • 


F,  =  0,  k  #  j 


Therefore  will  not  beccme  undefined. 
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Rsr  exairple,  if  a  unit  force  were  applied  at  station  21,  (i.e.,  =  1)  and 

all  other  forces  were  zero,  all  the  a  ^7  can  be  determined  from 

C.  ^_  =  X. 

1,27  1 

On  the  other  hand,  the  k2^  ^  are  not  as  conveniently  determined  since  it  is 
difficult  to  isolate  a  force  vdien  a  unit  displacement  is  applied  at  any 
station.  Consequently,  it  is  much  easier  and  more  accurate  to  measure  a 
displacement  than  to  measure  a  force  and  the  influence  coefficient  is  more 
often  used. 

For  the  equations 

[c]  [m]  fx}  +  fx}  =  0 

we  can  assime  sinusoidal  solutions  of  the  form 

Xj^  =  sin  ujt,  or 
[x3  =  [A3  sin  ut 

Substitution  into  Equation  12.40  leads  to 

2  -  -  - 

-a  [cj  [m]  [A3  sin  ut  +  [A3  sin  wt  =  0 


or,  since  sin  cut  ^  0 

1  [A3  =  [c]  [m]  [A3 

2 

0) 

This  is  a  standard  eigenvalue  problem  \duch  can  be  solved  using  iterative 
techniques.  Values  of  the  model  vectors  or  mode  shapes  A  can  be  found, 
corresponding  to  a  specific  frequency  u.  Ihere  are  as  many  eigenvalue 
solutions  0)  and  corresponding  eigenvectors  A  as  there  are  degrees  of  freedom 
in  the  ^stem.  Ohat  is,  an  n  degree  of  freedom  system  has  n  mode  shapes  and 
frequencies.  The  concept  of  mode  shapes  is  so  basic  to  flutter  analysis  tliat 
further  discussion  is  justified  at  this  time. 
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12.5.8  StructTjral  Vibrations  -  Mode  Shape  Determination 

The  basic  discussion  of  mode  shapes  and  frequencies  will  use  the  uniform 


cantilever  beam  of  length  L  shewn  in  Figure  12.83  as  an  exairrole. 


X 


FIGURE  12.83.  UNIFORM  CANTILEVER  BEAM 

T-Jhat  do  the  natural  mode  shapes  represent  and  how  do  we  coirpute  them? 
Ihe  real  question  to  address  is:  If  the  beam  were  to  vibrate  freely  at  seme 
resonant  condition*  what  vibrating  shape  would  it  assume  and  what  would  the 
frequency  of  \dht:ation  be? 

Since  the  beam  is  continuous  it  has  an  infinite  number  of  mass  points 

and,  consequently,  an  infinite  nirriber  of  degrees  of  freedom.  We  state, 

without  proof,  that  there  are  also  an  infinite  number  of  distinct  vibration 

shapes,  each  of  vdiich  has  a  xmique  vibration  frequency. 

Ihe  Equation  of  Motion  for  this  s^^stem  leads  one  to  a  certplex  or  double 

eigenvalue  problem,  where  the  eigenvector  determines  the  modal  shape  of  the 

beam,  X  (y) ,  and  corresponding  eigenvalue,  w  ,  represents  the  square  of  the 

z 

vibration  frequency.  These  eigenvalues  are  derived  given  the  mass  and 
stiffness  distributions,  and  the  specific  boundary  conditions.  For  the  beam 

♦Definition:  Resonant  Condition  -  Vibrating  at  a  maximum  amplitude  in  phase 
with  an  oscillating  input  at  the  systems  undairped  natural  frequency, 


in  Figure  12.83  the  boundary  conditions  are:  zero  displacement  and  slope  at 
the  fixed  end  (y  =  0) ,  and  zero  shear  and  mcment  at  the  free  aid  (y  =  L) .  Ihe 
eigenvectors  are  solutions  to  harogeneous  equations  so  that  if  Xj^(y)  is  a  mode 
shape,  C  Xj^(y)  is  one  also.  As  a  result,  each  mode  shape  function  represents 
relative  displacements  along  the  beam.  Uie  absolute  values  are  determined 
frcm  the  initial  conditions. 

The  differential  equation  for  the  cantilever  beam  as  given  in  Reference  1 
is 


vhere 


El 


d^x 


+  m 


d  X 
dt.. 


=  0 


(12.41) 


E  =  modulus  of  elasticity 
I  =  area  mcment  of  inertia 
m  =  mass  distribution. 


Equation  12.41  is  a  separate  partial  differential  equation  with  a  solution 
x(y,  t)  a  product  of  a  function  of  y  only,  X(y) ,  times  a  function  of  thne 
only,  T{t) . 


x(y,t)  =  X(y)  T(t) 

Substituting  x(y,t)  into  Equation  12.41  yields 


or 


El 


A 


T  +  m 


d\ 

dt.. 


X 


1  d^T  ^  El  d^X  1 
T  dt2  m  dy^  X 


0 


Since  y  and  t  are  independent,  they  can  be  equated  to  a  separation  constant. 

In  this  example  the  separation  constant  will  be  the  square  of  the  mode  shape 
2 

frequency,  o)  .  The  result  will  be  two  linear  differential  equations. 


2 


(D 


T 


0 
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0 


El  d^X 

m  ^ 


2 


(I) 


X  = 


or, 


vdiere 


0 


a=./ir 

Y  m 

The  solutions  to  Equations  12.42  and  12.43  are 

T  =  A  sin  cct  +  B  cos  ojt 
X  =  C  sinh  y  +  D  cosh  Y 

+  E  sin  ^ 


vAiere,  A,  B,  C,  D,  E,  and  F  are  constants,  implying  the  boundary  conditions 
for  a  laniform  cantilever  beam:  X(o)  =0,  X'  (o)  =0,  X"  (L)  -  0,  and  X'"  (L)  = 
0,  yields  the  following  transcendental  equation.  The  solution  of  vhich  gives 
the  mode  shape  frequenc;^'- 


v^iere 


(12.44) 


The  solution  to  Equation  12.44  can  be  found  hy  iterative  techniques  or  through 
graphic  depiction.  What  follows  is  the  graphical  solution  that  satisfies 
Equation  12.44. 
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Bie  abcissas  of  the  points  of  intersection  of  the  curves  yield  values  for 

Vi- 


L  =  0.597  Ti, 


1.49 


or 


(0.597)^ 


a 


(1.49)^ 


(i  sufficiently  large) . 


FIGURE  12.84. 


GRAPHICftL  SCLUTION  TO  THE  TRANSCENDEOTAL  EQUATION 
OF  A  IWIFORM  CANTILEVER  BEAM 
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Having  determined  the  model  frequencies,  the  specific  mode  shapes,  X^(y) , 
can  be  derived  by  substituting  the  frequency  values  into  the  general  mode 
shape  equation  for  the  cantilever. 


Where  D  is  a  constant  determined  by  the  boundary  conditions.  Do  not  be 
concerned  at  how  these  equations  vere  derived.  They  are  sinply  the  solutions 
of  the  fourth  order  differential  equation  for  the  vibrating  cantilever  beam. 
The  inportant  thing  to  note  is  that  there  are  an  infinite  number  of  solutions 
and  that  for  each  o)^  there  corresponds  an  X^(y).  Figure  12.85  shows  four 
of  the  lowest  frequency  mode  shapes. 

Ihe  significance  of  the  application  of  mode  shapes  and  frequencies  is 
that  any  \dhration  of  any  elastic  body  will  be  a  summation  of  individual  rrade 
shape  vibrations.  The  characteristics  of  mode  shapes  are  such  that  (1)  for 
each  natural  vibration  frequeno^  there  exists  one  and  only  one  natural  mode 
shape,  (2)  in  the  vibration  of  any  single  mode,  the  displacements  of  the 
structure  reach  their  zeros  and  their  extreme  simultaneously  [as  is  evident 
fron  the  solution  x(y,t)  =  X(y)  T(t)],  (3)  the  natural  mode  shapes  are 

linearly  independent  (i.e.,  no  shape  can  be  formed  by  any  linear  ccrnbination 
of  the  others) ,  and  (4)  in  any  free  vibration  of  a  beam,  wing,  aircraft 
structure,  bridge,  and  so  forth,  any  combination  of  natural  modes  can  exist 
simultaneously  without  mutual  interference.  In  the  case  of  aircraft  flutter, 
the  mode  shapes  beoone  coupled  due  to  the  action  of  the  external  aerodynamic 
driving  forces,  vMch  are  in  turn  functions  of  the  critical  mode  shapes. 
Using  the  general  discrete  mass  free  body  problem  discussed  earlier,  we  can 
relate  to  a  specific  aircraft  example. 


[c]  [m]  X  +  x  =  0 
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vd-th  solution 


to]  [m]  [A} 


Ihe  "wing"  in  Figure  12.82  has  five  separate  stations,  each  with  mass  m^ 
and  each  having  two  degrees  of  freedom,  x^,  a  vertical  displacement,  and  an 
angular  rotation  about  the  elastic  axis  as  shown  in  Figure  12.86. 


X 


FUNDAMENTAL  MODE 


1  =  2 


(0.597)2  (^)*a 
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a 


FIGURE  12.85.  FREQUENCY  MODE  SHAPE  FOR  A  VIBRATING  BEAM 


The  model  vector  [A3  is  corprised  of  ten-degrees  of  freedom;  five 
vertical  displacements  and  five  angular  deflections. 
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FIGURE  12.86.  CKJJTILEVER  VITl-JG 

From  experimental  measurements  we  can  determine  the  mass  and  mcment  of  inertia 
of  each  section,  along  with  the  influence  coefficients  (D^j)  for  the  total 
wing  with  respect  to  the  ten  degrees  of  freedan.  Recall  that  €2^  represents 
the  vertical  deflection  X2  at  wing  Station  2  due  to  a  unit  mcnent  applied  at 
wing  Station  3  (a^  being  the  ei^th  degree  of  freedom  by  definition) . 

The  matrix  equation 


1  CA3  =  [c]  [m]  £A3 


thus,  represents  10  simultaneous  equations  in  the  11  unkncwns  to,  X^,  .  .  . 
a^.  Again,  the  mode  shapes  are  relative  displacements  in  the  A  matrix, 
typical  solution  might  be 


12.150 


0)  =  3.5,  A 


1.00 

1.21 

1.63 

2.23 

3.05 

0.74 

1.13 

1.74 

2.57 

3.49 


At  any  given  time,  therefore,  the  wing  station  deflections  for  the  given  3.5 
rad/ sec  mode  might  be  represented  by 


{x3  =  K  [A3  sin  (3.5t  +  (f.^) 


vAiere  the  values  K  and  (j)^  would  depend  on  the  initial  conditions.  Other  modes 
at  different  frequencies  will,  of  course,  exist  simultaneously.  Flutter  is 
the  result  of  the  aerodynamic  coupling  of  these  modes. 


Even  for  simple  structures  the  mode  shapes  and  frequencies  may  be 
difficult  to  determine  analytically.  As  a  result,  many  simplifying  assumptions 
are  normally  made.  For  the  wing  just  described,  the  assumption  that  is  can  be 
represented  by  five  rigid  sections  seems  a  bit  gross,  but  for  ttie  lo.v 
frequency  mode  shapes  the  simulation  is  quite  reasonable.  In  actual  aircraft 
problems  the  critical  modes  are  usually  the  lower  frequency  ones  for  very 
fundamental  reasons.  The  large  bending  and  structural  deformations  are  the 
ones  which  more  efficiently  extract  energy  from  the  airstream.  Ihese  are 
usually  the  lower  frequency  modes. 


A  useful  technique  for  the  flutter  analyst  is  to  assume  that  a  structure 
will  taJce  on  certain  predetermined  mode  shapes  based  on  previous  observations 
of  similar  type  structures.  Ihe  assuined  mode  shape  can  then  be  used  in  the 
energy  equations  to  arrive  at  the  particular  mode  shape  frequency.  Ihis 
technique  aids  in  the  solution  of  complex  structures  but  there  is  obviously  an 
acoonpanying  reduction  in  accuracy.  The  subject  of  mode  shapes  will  be 
broui^t  up  again  in  the  discussion  of  ground  vibration  tests. 
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12.5.9  Wind  Tunnel  Modeling 

The  use  of  aircraft  models  has  produced  solutions  to  practical  problems 
in  the  areas  vAiere  existing  theory  is  not  yet  dependable.  Particularly  in 
dealing  with  flutter,  the  testing  of  wind  tunnel  models  with  properly  scaled 
mass  and  stiffiiess  properties  has  often  yielded  better  results  than  equivalent 
analytical  efforts. 

12.5.10  Buckingham  ir  Theorem 

Model  theory  must  be  based  on  a  clear  understanding  of  the  principles  of 
dimensional  analysis.  The  writing  of  equations  in  dimensionless  form  with  a 
reduction  in  the  number  of  variables  was  generalized  by  Buckingham  in  his  "it 
Theorem"  vMch  states  that  if  a  pl^sical  situation  can  be  represented  by  the 
equation 

ij;  (S^,  S2,  S^,  .  .  .  ,  S^)  =  0 

vdiere  the  n  arguments  include  all  the  primary  quantities  (mass,  length, 
time,  etc.),  the  secondary  quantities,  and  dimensional  constants  which  must  be 
considered  in  the  problem,  the  equation  can  be  rewritten  in  the  form 

♦  <”1'  '2'  ’3 . W  -  0 

in  vdiich  ir, ,  7r_,  . .  ir  are  the  (n-m)  independent  products  of  the 

i  z  n-iri 

argunents  S^,...,  S^,  vdiich  are  dimensionless  in  primary  quantities.  The  form 
of  these  dimensionless  ir's  can  be  found  by  a  formal  procedure  but  they  can 
usually  be  constructed  by  inspection.  Typical  ir's  in  general  use  are  aspect 
ratio  AR,  reduced  frequency  k  =  ub/V,  Mach  M,  and  Reynolds  number  Re.  There 
are  two  main  advantages  in  using  dimensionless  variables.  First,  since  the 
dimensionless  equation  of  motion  is  cortpletely  unaffected  by  scale  effects  the 
values  of  the  dimensionless  variables  should  be  the  same  for  both  the  original 
problem  and  its  model. 

As  a  sanple  problem  we  can  use  the  wing  twist  relationship  derived 
earlier  to  determine  the  wing  divergence  speed. 

.  ^  q  s  fe  t  C 

"^6  -  9  Se 
a 
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6,  Cj^  and  are  alreac^  dimensionless.  The  other  quantities  have  the 


dimensions  as  indicated  (n  +  5) . 


Kj  -n. 


S  -  L' 


6  •"  L 


C  “  L 


Force  (F)  and  length  (L)  v^e  chosen  as  the  primary  quantities  (m  =  2) . 

By  inspection  we  can  form  independent  dimensionless  ratios  e/c  and  qSc/Kg 
so  that  the  equation  may  be  rearranged  to  read 


r)  [(f) 


Uie  number  of  arguments  has  been  reduced  by  two,  as  predicted  by  the 
Theorem  (since  there  are  the  two  primary  quantities  F  and  L  in  the  equation) . 
The  equation  really  inplies  that  any  model  of  the  semi  rigid  wing  having  the 
same  shape  and  0^  )  must  have  the  same  dimensionless  location  of  the 

elastic  axis  (e/c) ,  the  same  rigid  angle  of  attack  (a_) ,  and  the  same  ratio  of 
aerodynamic  to  elastic  forces  (qSc/Ke) .  It  is  also  desirable  to  have  similar 
Mach  and  Reynolds  numbers  if  the  effects  are  significant. 


12.5.11  Aeroelastic  Model 

All  aeroelastic  models  are  designed  along  three  fuixiamental  airplane 
properties:  (1)  structural  stiffness  distribution,  (2)  mass  distribution,  and 
(3)  the  external  shape.  The  model  designer  has  decide  befordiand  vhich  of 
the  properties  will  require  a  more  exact  simulation,  depending  on  the  nature 
of  the  planned  tests.  Property  (1)  is  important  for  accurate  loading 
measurements  and  property  (3)  is  important  for  aerodynamic  force  evaliiations. 
A  flutter  model  obviously  requires  an  accurate  reproduction  of  all  three 
properties  and  is  therefore  difficult  to  design. 
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The  usual  first  approach  is  to  attenpt  a  scaled  replica.  However,  for 
low  speed  models  the  required  skin  thickness  (using  the  same  material)  may  be 
prohibitively  small.  For  a  more  workable  thickness  a  softer  material  must  be 
used  but  this  in  tram  tends  to  reduce  the  accuracy  of  the  stiffness 
reproduction.  It  also  limits  the  amount  of  structural  detail  vbich  can  be 
obtained.  The  designer  most  often  finds  that  the  best  plan  is  to  lay  out  a 
simplified  structure,  making  sure  that  it  does  not  use  up  too  much  of  the 
available  mass.  A  shell  is  then  formed  to  enclose  the  structure  with  the 
proper  external  shape.  Finally,  the  remaining  mass  is  distributed  over  the 
sections.  Ihe  type  of  construction  used  for  both  the  structure  and  the 
external  shell  depends  upon  the  size,  speed  range,  and  ratio  of  the  air 
density  between  aircraft  fli^t  altitude  and  test  chamber  conditions. 

The  types  of  aeroelastic  testing  fall  into  three  general  areas,  vbether 
with  models  or  with  full  scale  airplanes.  The  first  area  requires  no 
airstream,  such  as  in  fatigue,  static  loading,  and  ground  vibration  tests. 
The  second  and  third  test  areas  involve  airflow,  in  either  the  tunnel  or  full 
scale  flight.  The  second  area  enconpasses  the  static  phenanena  divergence, 
(control  effectiveness)  and  the  third  area  includes  the  unsteady  phenanena 
such  as  flutter,  gust  loading  and  (^amic  stability. 

Ihe  test  programs  for  nearly  all  prototypes  include  "shake  testing",  to 
determine  the  normal  modes  of  vibration.  These  derived  mode  shapes  and 
frequencies  can  be  coipared  with  the  analytical  calculations  to  verify  the 
mass  and  stiffeess  properties  of  the  math  model.  Sonetiraes  the  experimental 
data  serve  as  the  basis  for  a  new  set  of  calculations. 

Once  an  adequate  model  is  built  the  testing  takes  place.  There  are 
basically  different  approaches  to  the  problem.  In  the  first,  the  testing 
is  designed  to  evaluate  the  coefficients  in  the  differential  equation 
governing  the  problem.  In  the  second  approach  the  model  is  designed  and 
tested  as  an  analog.  The  model  can  simulate  parts  of  airplanes,  such  as  wings 
or  tails,  or  the  vbole  airplane.  Each  schene  has  peculiar  advantages  and 
drawbacks  as  well  as  the  usual  problems  of  excitation,  moimting  and 
measurenent.  Ihe  particular  details  of  tunnel  testing  and  measurement  is  a 
total  subject  in  itself. 


12.154 


12.5.12  Wind  Tunnel  Model  Flutter  Prediction  Methods 

Four  methods  used  to  measure  the  siibcritical  (below  the  actual  flutter 
speed)  re^nse  characteristics  are  co/quad,  randcradec,  power  spectral 
density  (PSD),  and  peak-hold  spectrum  methods.  These  methods  are  used  to 
measure  the  frequency  and  damping  (or  an  inverse  response  amplitude 
proportional  to  the  damping  in  the  peak-hold  spectrum  case)  in  the 
predondnant  or  critical  vibration  modes.  Ety  suitably  plotting  and 
extrapolating  the  siibcritical  damping  in  the  vibration  mode  or  modes  of 
interest,  the  flutter  point  can  usually  be  established.  With  each  method,  the 
response  can  be  approximated  by  that  of  a  single-degree-of-freedcm  system. 
All  of  these  methods  can  be  used  real-time,  that  is,  used  to  translate  the 
response  time  history  samples  into  quantitative  information  for  the  test 
engineer  vhile  the  test  is  in  progress. 

Briefly,  the  co/quad  method  measures  the  in-jhase  and  out-of-phase 
components  of  the  forces  re^xjnse  generated  by  the  sinusoidal  frequency  sweep 
technique.  The  randondec  method,  a  relatively  new  method  used  in  the  F-16 
flight  flutter  testing,  makes  use  of  ensemble  averaging  of  transient  response 
to  random  excitation.  The  PSD  method  is  a  well-known  procedure  for  the 
analysis  of  random  refuse  data.  It  is  obtained  directly  from  an  ensemble 
average  of  the  square  of  the  magnitude  of  the  Fourier  transform  of  a  number  of 
segments  of  the  time  history.  In  the  peak-hold  spectrum  method,  Fourier 
components  of  a  number  of  time  history  segments  are  determined  and  the 
envelope  of  the  peak  values  of  these  components  is  obtained  as  a  function  of 
frequency. 

Co/Quad  Method.  The  co/quad  method  involved  measuring  the  forced 
response  of  a  model  to  an  input  force  such  as  that  generated  by  a 
trailing-edge  control  surface  as  illustrated  schamatically  in  Figure  12.87. 
If  a  transfer  function  relating  the  response  to  the  input  force  is  determined 
as  a  function  of  frequenoy,  then  the  damping  on  each  mode  can  be  cistained. 
Cross  spectrum  between  the  control  surface  command  signal  and  the  model 
(tynamic  response  can  be  determined  with  a  co/quad  analyzer.  This  "type 
analyzer  presents  two  outputs  in  terras  of  in-phase  (called  co  for  coincident) 
and  out-of-phase  (called  quad  for  quadrature)  components  between  signals. 
Several  means  of  calculating  the  damping  are  available  directly  from  a  co  and 
quad  type  of  presentation.  As  indicated  in  Figure  12.87,  the  damping  of  a 
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mode  can  be  estimated  from  the  out-of-phase  conponent  by  the  frequencies 
labeled  f^  and  fg.  These  are  the  frequencies  at  the  half-power  points  and  the 
structiaral  daitping  g  can  be  espressed  in  terms  of  these  frequencies  (Figure 
12.87). 


SWEEP 

CONTROL 

SURFACE 

r 

MODEL 

1 

OSCILLATOR 

ACTUATOR 

DYNAMICS 

1  DYNAMICS  1 

COMMAND  SIGNAL  6^ 

h 

_ _ \!^ 

1 

r- 

11 

“1 

TRACKING 

FILTERS 

-  1  ih 

I'PHASE 

1  : 

r 

. 

1 

h 

6c 

.  OUT-OF- 

a  _  IV*b1^  - 

ANALYZER 

1 

r 

1 

1  .  .  . 

X-Y 

s 

10  15  20 

25 

PLOTTER 

FREQUENCY, 

Hz 

FIGURE  12.87.  IMPLEMENTATICN  OF  CO/QUAD  METHOD 
Bandoitidec  Method.  The  Randan  Decren:ent  method  is  basically  an  ensenible 
averaging  of  the  turbulence-induced  random  vibrations  of  the  test  article.  As 
is  illustrated  in  Figure  12.88,  triggering  each  data  sample  at  a  constant 
level,  Y^.  Assuming  linear  superposition,  the  time  history  of  each  sample  can 
be  regarded  as  the  combined  solution  from  (1)  an  initial  step  displacement, 
(2)  an  initial  velocity  and  (3)  a  randan  forcing  function.  Note  that  the 
Figure  12.88C  sample  represents  the  reponse  to  the  same  initial  displacement 
as  Figure  12.88B,  a  different  initial  velocity  with  the  opposite  sign,  and  a 
different  randan  forcing  function.  It  can  be  reasoned  intuitively  that  vAien  a 
large  number  of  samples  are  averaged,  only  the  re^nse  to  the  constant 
initial  displacement  will  remain  because  the  average  of  responses  due  to  the 


12.156 


alternating  initial  velocities  and  the  randan  forcing  functions  vdll  tend  to 
zero.  Ihus,  it  is  seen  that  the  ensertble  average  converges  toward  the 
transient  re^onse  to  an  initial  step.  R>r  a  constant  trigger  level,  the 
ensenble  average  (Randondec  Signature)  will  be  constant  even  if  the  airplitude 
of  the  forcing  function  varies.  If  the  ensemble  average  is  made  rp  of  sairples 
with  initial  positive  slopes  only,  then  the  resulting  trace  r^resents  the 
transient  re^x>nse  to  a  conbined  step  and  initial  velocity.  Under  these 
conditions  the  Randcmdec  Signature  would  vary  with  the  intensity  of  the  forcing 
function,  thus  minimizing  the  use  of  the  signature  trace  as  a  failure 
detector.  However,  the  dairping  as  determined  fron  the  decay  rate  of  the 
signatuire  trace  would  be  valid. 


(•)  (c) 


(RAN  DOM  DEC  SIGNATURE) 
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FIGURE  12.88.  RANDOM  DEa®y!E3SlT  CmCEPT 
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Structiaral  danping,  g,  may  be  detennined  directly  by 


PSD  and  Peak-Hold  Spectrum  Methods.  Ihe  PSD  and  the  peak-hold  spectrum 
methods  are  inplemented  as  shown  in  Figure  12.89.  Both  methods  are 
implemented  using  a  spectroscope.  This  analyzer  etiploys  time  ccttpression 
techniques  to  achieve  minimum  analysis  time  for  the  frequency-tuned  band-pass 
filter  to  convert  the  irput  signal  from  the  time  domain  to  the  frequency 
domain.  Following  conpression,  the  input  signal  is  frequency  analyzed.  Shown 
on  the  left  of  Figure  12.89  is  a  typical  PSD  obtained  from  the  model  dynamic 
response  h.  The  resulting  signature  has  a  peak  for  each  structural  mcxie  and, 
for  well-separated  peaks,  the  damping  ratio  may  be  obtained.  As  indicrated  in 
Figure  12.89,  the  structural  damping  is  equal  to  the  frequency  bandwidth, 
taken  at  the  half-power  point,  and  divided  by  the  mode  frequency. 

An  additional  mode  of  operation  of  a  spec:troscope  allc^ws  for  detection 
and  storage  of  the  peak  values  of  frequency  windows.  In  this  mode  of 
operation,  an  ensemble  specrhrum  composed  of  frequency  windows  is  obtained. 

Upon  receipt  of  each  subsequent  spectrum,  peak  filter  response  at  each 
lcx:ation  is  updated  in  a  positive  direction.  That  is,  only  an  increase  in 
value  clauses  an  update  to  the  new  hi^er  value.  On  the  right  of  Figure  12.89, 
a  typical  peak-hold  spectrum  is  shown.  With  this  method  the  damping  parameter 
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FIGURE  12.89.  IMPLEMEOTATICW  OF  SPECTRUM  METHODS 

is  not  obtained.  However,  the  reciprocal  of  the  peak  spectrum  aitplitude  1/P 

is  proportional  to  the  darrping  ratio  and  is  used  as  a  measure  of  system 

stability.  The  peak-hold  method  can  be  applied  using  two  forms  of  excitation, 

model  response  to  tunnel  turbulence  and  model  response  to  sinusoidal  force. 

Typical  results  obtained  from  the  four  subcritical  response  methods  are 

presented  and  coipared  in  Figures  12.90  and  12.91.  Figure  12.90  presents  the 

variation  of  structxiral  dairping  coefficient  of  a  delta  wing  flutter  model  with 

dynamic  pressure.  The  damping  results  obtained  with  co/quad,  randomdec,  and 

PSD  are  indicated  with  open  symbols.  The  model  fluttered  at  a  dynamic 

2 

pressure  of  5.89  kPa  (1231bf/ft  )  as  indicated  with  the  closed  symbol.  A  plot 
of  the  inverse  amplitvide  of  the  peak  spectrum  (vised  as  the  stability  criteria) 
is  presented  in  Figure  12.90B  as  a  function  of  c^amic  pressure.  Shown  are 
results  from  forced  excitation  and  random  excitation  (turbulence) . 
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(b)  PEAK-HOLD  SPECTRUM  RESULTS 

FIGURE  12.90.  COI'JPARISON  OF  SUBCRITICAL  METHODS,  DELTA-T-JING  MODEL  (M  =  O.SO) 

Further  illustration  of  the  type  of  data  generated  with  the  use  of  the 

four  suhcritical  methods  is  presented  in  Figure  12.91.  Shewn  are  the  data 

plots  frctTi  vhich  the  danping  levels  presented  in  Figure  12.90  were  obtained. 

Ihe  wind-tunnel  conditions  were  the  same  for  each  method  (Mach  M  =  0.90; 

2 

dynamic  pressure  q  =  5.42  kPa  (113  Ibf/ft  ) . 

12.5.13  Ground  Vibration  Testing  (GVT) 

As  mentioned  earlier,  the  primary  purpose  of  ground  vibration  testing 
(shake  testing)  is  to  measure  the  structural  mode  shapes  and  frequencies. 
Sometimes  the  data  serve  to  generate  a  new  set  of  calculation.  At  other 
tiroes,  the  measured  mode  shapes  serve  as  the  model  for  flutter  predictions, 
especially  vhen  the  effects  of  configuratian  changes  on  flutter  speeds  are 
being  investigated. 
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FIGURE  12.91.  niiUSTRATION  OF  SUBCRITICMi  METHODS 
M  =  0.90;  g  =  5.42  kPa  (113  Ibf/ft^) 


The  X-24B  "lifting  body"  is  an  excellent  example  of  a  reconfigured 
aircraft  (frcm  the  X-24A) .  The  structural  modifications  include  the  addition 
of  ri^t  and  left  strake/aileron  cor±>inations  and  the  modification  of  the  two 
outboard  vertical  fins.  Ihe  fin,  strake  and  aileron  form  a  single 
sub-assembly  mated  to  the  airframe.  From  analyses  and  previous  test 
experience  with  the  X-24A,  the  only  possible  flutter  problems  would  occur  with 
the  strake/aileron  and  fin/rudder  carponents.  Ftor  further  definition  of  the 
model  characteristics  of  these  carponents  a  ground  vilaraticai  test  program  was 
conducted. 

Uie  objective  in  the  program  (as  in  all  GVT's)  was  to  measure  the 
frequency,  mode  shape  and  stojctural  danping  for  each  significant  mode  of  the 
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new  tail  section  confLgtiration.  E^ecially  inportant  vas  the  degree  of 
interaction  between  the  strake/aileron  and  the  fin/rudder  in  the  vibration 
ncdes. 

For  GVT,  there  is  always  the  requirement  to  isolate  the  rigid  boc^  inodes 
so  that  the  hipest  frequency  of  rigid  bcxty  motion  is  well  below  the  lowest 
frequency  of  aiy  strucrtural  vibration  or  mode  shape.  For  exainple,  if  the 
natural  frequency  of  rotation  of  the  elevators  were  10  Hz,  the  aircraft  as  it 
is  supported  should  not  have  a  pitch  frequency  close  to  the  10  Hz  value. 
Otherwise  there  may  be  a  problem  in  properly  isolating  and  in  determining  the 
existence  of  that  particular  elevator  rotation  mode. 

Ideally,  the  way  to  isolate  the  rigid  bo^  modes  is  to  suspend  the 
aircraft  throu^  the  c.g.  by  a  long  cable  but  the  methods  cannot  conveniently 
be  aiplcyed.  Two  practical  techniques  are  (1)  to  support  the  vehicle  on  air 
cushion  stands,  or  (2)  to  suspend  the  aircraft  with  soft-spring  ceiling 
mounted  su^)ension  systems.  Ihe  approach  used  on  the  X-24B  was  to  deflate  the 
tires  to  half  pressure. 

The  way  to  excite  the  structural  modes  on  ground  vibration  tests  is 
through  seme  variable  frequency  vibrating  device  attached  to  the  structinre. 
Typically  used  are  constant  force  electromagnetic  or  hydraulic  shakers  on 
vhich  there  is  precise  frequency  tuning  and  an  automatic  frequency  sweep 
feature. 

For  the  X-24B,  elec±romagnetic  shakers  with  current  feedback  and 
individual  gain  and  phase  controls  were  used  to  input  the  sinusoidal  forcing 
function  to  the  structure.  Seven  excitation  conditions  involving  six 
different  shaker  locations  vere  required  to  carpletely  investigate  the  model 
characteristics.  Seven  accelercmeters  were  used  to  obtain  data  during  the 
test.  Acceleration  measurements  v^re  recorded  at  predetermined  grid  points 
marked  on  the  structure. 

The  test  procedure  was  to  first  conduct  a  frequency  sweep  from  10  to  200 
Hz  at  each  of  six  different  shaker  locations.  Possible  modes  were  identified 
by  observing  Lissajous  patterns  on  an  oscilloscope  of  input  force  versus 
reference  acceleration,  and  of  other  acceleration  ratios.  Essentially,  at 
resonant  condition  the  force  would  be  in  phase  with  the  velocity  and  this 
situation  would  be  indicated  by  the  scope  patterns.  Once  a  mode  was  spotted 
the  shaker  covild  be  retuned  more  precisely  to  the  resonant  frequency. 
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Acceleration  and  phase  measurements  were  than  taken  along  the  predetermined 
grid  locations  labeled  on  the  structure.  Structural  damping  values  were 
i  Tipd  by  dumping  the  armature  current  to  the  shakers  and  measuring  the 
vibration  decay  traces. 

The  shape  and  frequency  data  would  be  used  in  a  flutter  analysis 

vterdDy  the  aerodynamic  forces  (through  a  Theodorsen  type  analysis)  woiiLd  be 
evaluated  in  their  interaction  with  the  stjnacture.  Any  conditions  of 
instabilil^  could  be  predicted  and  the  results  used  to  refine  the  V-g  diagrams 
or  equivalent. 

12.5.14  Flight  Test 

Most  of  the  preflight  functions  have  been  accounted  for.  Vfe  can  new 
approach  the  subject  at  the  user's  level  and  discuss  the  subject  of  flight 
test.  Is  the  aircraft  free  of  aeroelastic  problems  and  if  so,  to  v^at  extent? 

The  criteria  for  strength  and  dynamic  instabilities  are  listed  in 
Reference  12.18.  V^,  the  limit  speed  for  the  basic  and  hi^  drag 

configurations,  is  the  maximum  attainable  speed  cenmensurate  with  the 
operational  use  of  the  airplane,  considering  shallow  and  steep  dive  angles, 
thrust,  operation,  and  nonoperation  of  speed  brakes,  and  inadvertent  upsets 
from  gusts,  or  as  specified  in  the  contract  documents.  The  airplane  or  its 
oonponents  shoxald  not  exhibit  flutter,  buzz,  divergence  or  other  aeroelastic, 
aerothermoelastic  or  aeroservoelastic  instabilities.  The  fifteen  percent 
safety  margin  shall  be  shown  by  analytical  or  experimental  data  (including 
fli^t  test  up  to  Vj).  In  addition,  the  damping  coefficient  g  for  any 
critical  flutter  mode  or  any  significant  dynamic  response  mode  shall  be  at 
least  three  percent  (0.03)  for  all  altitudes  and  flight  speeds  up  to 

Reference  12.18  also  furnishes  guidelines  on  how  the  math  analyses  are  to 
be  performed,  specifying  the  use  of  (1)  corpressible  aeroc^mamics  in  high 
subsonic  flight,  (2)  finite  span  assumptions,  and  (3)  three-dimensional  flow 
effects,  if  significant.  Model  and  ground  vibration  tests  are  also  specified 
vAien  necessary. 

For  the  project  manager  or  project  pilot,  then,  what  should  the  important 
considerations  be  in  the  formulation  of  the  test  program.  The  driving  force 
should  be  the  probability  of  the  flutter  condition  and  its  anticipated 
seriousness.  For  example,  if  the  airplane  were  apt  to  destrey  itself  well 
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inside  its  performance  envelope  the  program  would  be  handled  differently  from 
that  of  a  flutter  clearance  demonstration  vdiere  the  predicted  safety  margin 
was  50%. 

If  the  program  were  solely  for  flutter,  the  problem  would  probably  be 
very  specific.  In  a  larger  protoi^^  eval  program,  flutter  investigations  are 
usually  clearance  demonstrations,  althou^  there  mi^t  be  occasions  vAien 
specific  problems  are  anticipated. 

At  a  critical  flutter  condition  the  dartping  ratio  of  one  of  the  modes 
beccmes  zero.  The  object  of  fli^t  flutter  testing  is,  then,  to  obtain  a 
measure  of  the  dairping  values  associated  with  the  flight  modes  and  their 
trends  as  the  airspeed  is  progressively  increased  so  that,  from  the  danping 
trends  as  sub-critical  conditions,  the  approach  of  a  critical  condition  may  be 
indicated  and  its  speed  determined  by  extxapolation. 

Prototype  airplanes  represent  tremendous  investments  in  time  and  money. 
Extensive  analytical,  iiKxiel  and  shake  tests  are  performed  in  the  design  and 
developnent  phase  so  that  for  most  of  the  newer  aircraft  the  flutter  margins 
will  have  been  assured.  If  not,  the  predictions  on  speed  and  danping  for  the 
prc±)lem  areas  will  be  reasonably  accurate.  On  the  initial  flints  the  tests 
will  be  carried  out  to  increasingly  hi^er  speeds  and  an  actual  occurrence  of 
flutter  can  be  spectacularly  destructive.  Even  though  many  preventive 
measures  m^  have  been  taken  by  the  designers  the  possibility  of  such  an 
occurrence  cannot  be  overlcxjked. 

For  most  filter  type  aircraft  the  problon  areas  will  not  be  with  the 
basic  airplane  but  instead  with  the  various  wing  stores  configxorations.  On 
large  bcrtiber  or  transport  aircnraft  an  aeroelastic  instability  m^  exist  on 
the  basic  airplane.  Regardless  of  whether  the  aircraft  is  new  or  old,  large 
or  small,  the  primary  consideration  in  the  program  fonnulation  should  be  the 
seriousness  and  probability  of  the  flutter  condition. 

Other  considerations  in  the  program  development  are  time  and  monetary 
constraints.  Is  the  problem  serious  enough  that  the  airplane  may  face  either 
a  grounding  or  a  severe  restriction  on  the  envelope  until  a  solution  is  found? 
A  grounding  of  the  entire  T-37  fleet  due  to  a  potential  rudder  flutter  problem 
would  place  an  unacceptable  restriction  on  the  USAF  UPT  program.  An  airspeed 
restriction  on  a  front  line  attack  aircxaft  due  to  a  possible  wing  flutter 
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FIGUEE  12.92.  HYPOTHETICAL  V-g  DIAGRAM 


prc±!lem  would  be  intolerable  during  a  crisis  situation.  In  either  case  there 
is  an  iniTiediate  need  to  find  a  solution  -  as  was  the  case  for  the  w-ing 
divergence  prcblem  of  the  Fokter  D-8  during  WW  I. 

Along  with  the  time  constraint  is  the  problem  of  money.  Are  there  enough 
funds  to  conduct  a  thorou*^,  ground  vibration  test  (tine  permitting)?  What 
about  money  to  contract  for  an  in-depth  mathematical  analysis?  These  factors 
are  discussed  separately  siitply  for  clarity. 

In  preparing  the  test  plan,  of  prime  iitportance  is  attitude  and 
Tmderstanding  of  the  details  of  the  testing  to  be  performed.  Are  all  test 
^5ecifications  and  the  pertinent  aircraft  characteristics  totally  familiar? 
Can  critical  speeds  be  achieved  in  level  fli^t?  Will  there  be  time  at  the 
test  points  for  adequate  frequency  sweeps  by  an  on-board  excitation  systen? 
Do  the  handling  qualities  change  significantly  in  the  transonic  range?  Will 
the  on-board  method  of  excitation  will  be  adequate? 

During  the  late  30 's  the  Germans  encountered  destructive  flutter  on 
several  flight  tests  due  to  their  inability  to  detect  its  onset.  The 
recording  equipment  was  not  satisfactory  in  obtaining  the  required  data  at  the 
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time  needed.  The  excitation  equipnent  was  inadequate  in  properly  exciting  and 
controlling  the  critical  modes  at  svibcritical  speeds.  It  is  extremely 
inportant  that  the  approach  to  a  critical  flutter  condition  be  recognized 
observing  the  subcritioal  airplane  response.  Damping  trends  beccme  very 
important  as  the  ^speeds  are  increased  toward  the  potential  flutter 
condition,  especially  if  the  V-g  diagram  is  of  the  type  shown  in  Figure  12.92 
where  the  slope  steepens  quickly  at  the  critical  condition. 


12.5.15.  Ih-flight  Excitation 

To  assure  that  the  critical  mode  is  adequately  excited  at  subcritical 
speeds,  a  special  excitation  system  may  have  to  be  manufactured.  The 
following  methods  for  inrfli^t  excitation  presently  exist: 


(1)  Manual  Pulses  (stick  raps)  -  can  e3«:ite  frequencies  up  to  six  or 

seven  HZ  or  with  pcwered  control  systems,  up  to  10  to  12  Hz 
depending  on  aircraft  size  and  servo  capabilities.  The  main 

advantages  of  this  technique  are  convenience  and  cost.  A  quick, 
solid  rap  with  the  open  palm  to  the  side  of  the  stick  is  the  method 
usually  ettployed.  It  was  used  on  the  F-8  supercriticcil  wing  flutter 
investigation.  The  same  technique  was  used  with  the  TACT-F-111 
except  that  a  wooden  mallet  was  used  instead  of  the  open  hand.  The 
main  disadvantages  are  lack  of  frequency  selection  and  the  inabili'^ 
to  excite  the  hi^er  frequency  modes.  It  is  difficult  to  manually 
simulate  an  impulse  input.  The  energy  content  in  the  high  frequency 
modes  is  usually  deficient. 

(2)  Ballistic  Charges  "Bonkers"  -  pulse  charges  have  been  used  to  excite 
flutter  modes.  These  mini-esq)losions  more  closely  simulate  the 
ideal  impulse  with  its  large  frequency  content.  The  disadvantage  is 
that  these  charges  are  usually  one  shot  devices,  if  pyrotechnic. 

(3)  Sinusoidal  Shakers  (inertial  Exciters)  —  these  devices  consist  of 
either  a  rotating  out-of-balance  vheel  or  a  mass  wand  oscillating  at 
a  desired  frequency.  The  obvious  advantage  is  frequency 
selectivity.  The  disadvantage  of  inertia  exciters  becomes  apparent 
at  low  frequencies  (i.e.,  below  3  c.p.s.)  when  an  extremely  large 
out-of-balance,  or  equivalent,  is  required  to  produce  the  desired 
force  and  the  wei^t  of  the  exciter  may  becone  prohibitive.  Also, 
in  the  case  of  the  vheel  there  may  be  an  excess  of  shaking  force  at 
the  hi^  frequencies.  The  F— 5,  YF— 17  and  B— 1  prototypes  used  mass 
shaters  for  flutter  excitation. 

(4)  Oscillating  Vanes  -  In  this  method  an  axKiliary  aerofoil,  of 
synmetrical  section,  is  attached  externally  to  the  aircraft 
structTire  and  is  made  to  oscillate  with  a  sinusoidal  change  of 
incidence  by  a  variable-frequency  driving  mechanism.  The  method  is 
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best  suited  to  the  excitation  of  low-frequency  inodes.  In  order  to 
control  the  force  exerted  by  the  aerofoil  it  is  necessary  to  provide 
sane  type  of  force-^neasuring  like  between  the  aerofoil  and  the 
aircraft  structure  and  also  to  have  independent  control  of  incidence 
as  well  as  of  frequency.  Nevertheless,  it  may  be  difficult  to 
control  the  force  with  accuracy  in  the  transonic  region  vten  large 
changes  of  lift  m^  occur  with  small  changes  in  Mach.  Care  must  be 
taten  to  ensure  that  the  installation  of  the  aerofoil  does  not 
significantly  affect  the  flutter  characteristics  of  tiie  aircraft. 
Oscillating  Vanes  vrere  used  in  the  C-141,  C-5A,  B-52,  and  A-10  test 
programs. 

(5)  Autopilot/Autostabilizer  Excitation  -  As  the  state-of-the-art  is 

advanced  and  more  sophisticated  control  augmentation  systems  are 
developed,  increasing  dependence  will  be  place  on 
Autopilot/Autostabilizer  systems  for  flutter  excitation.  This 
method  m^  be  adopted  for  tests  on  aircraft  having  pcwered  flying 
controls.  A  sinusoidal  voltage  signal  of  variable  frequency  is  fed, 
as  an  error  signal,  into  the  autopilot  and  this  produces  appropriate 
oscillations  of  the  main  control  surfaces.  Alternatively,  the 
signal  may  be  fed  to  the  autostabilizer  unit  to  produce  the  same 
result.  Both  methods  will  be  effective  down  to  zero  frequency  but 
it  m^  be  difficult  to  excite  the  higher-frequency  modes  of  the 
aircraft  because  of  attenuation  \4iich  may  be  inherent  in  the 
transmissibilil^  characteristics  of  the  pcwer  control  system.  The 
advantages  are  obvious  but  the  control  system  response  plus  freeplay 
and  stiffness  must  be  very  accurately  known.  The  F-15  flight  test 
program  involved  use  of  the  CAS  for  flutter  excitation.  The 

Mvanced  Aerial  Refueling  Boon  (AARB)  flight  test  program  used  its 
fly-by-wire  control  system  to  program  frequency  sweeps  and  control 
inputs  to  excite  various  structural  modes. 

(6)  Turbulence  -  with  the  growth  of  statistical  analysis  applications, 
much  more  use  will  also  be  made  of  turbulence  and  other  atmosj^ieric 
phenomena  as  sources  of  excitation  for  flutter  testing.  A  spectral 
density  analysis  of  a  randan  input  can  be  analyzed  against  the 
output  to  find  the  correlation  function.  The  main  disadvantages  are 
that  there  is  no  control  over  the  amplitude  and  frequency 
characteristics  of  the  input,  and  that  turbulence  is  hard  to  find. 
Another  disadvantage  is  that  an  on-line  conputer  is  required  to  make 
^e  responses  meaningful,  real  time.  The  YF-16  test  program 
incorporated  the  use  of  turbulence  as  a  flutter  excitation  method. 


The  selection  of  the  excitation  method  is  very  important.  The  critical 
mode  must  have  sufficient  excitation  at  subcritical  speeds  so  that  the  damping 
behavior  can  be  observed.  This  point  cannot  be  overenphasized. 

Experience  has  shown  that  it  is  preferable  to  apply  the  excitation  as 
close  as  possible  to  the  surfaces  v^ose  flutter  characteristics  are  under 
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particular  investigation/  as  it  is  not  always  possible  to  obtain  adequate 
forced  airplitudes  of  the  surfaces  under  observation/  by  excitation  at  a  remote 
point  in  the  structxure.  This  will  mean  that  in  a  general  flutter 
investigation  is  may  be  necessary  to  provide  separate  excitation  equipnnent 
for  the  wings  and  tail  end  of  the  aircraft.  For  example/  if  inertial 
exciters  were  being  used/  it  mi^it  be  necessary  to  install  an  exciter 
outboard  in  each  wing  with  the  ability  to  operate  in  and  out  of  phase  and/  in 
addition/  independent  vertical  and  lateral  excitation  in  the  rear  fuselage. 
The  detail  positions  of  the  exciters  would  be  dictated  by  local  space  and 
strer^th  consideration,  always  bearing  in  mind  that  the  exciters  must  be 
p3^ced  as  far  as  possible  away  from  nodal  lines  for  the  modes  of  vibration  of 
interest  in  the  flutter  problem. 

Other  areas  which  mi^t  need  seme  attention  in  the  test  plan  fonrailation 
are  pilot— static  calibrations,  use  of  a  back-up  pilot  in  the  safety/photo 
chase  role,  and  so  forth.  The  test  plan  must  be  thorou^  and  explicit  with 
regard  to  all  responsibilities  and  functions.  Yet  it  must  be  realistic  enough 
that  no  one  has  any  problems  sticking  to  the  rules. 

During  the  early  50' s  at  Wright-Patterson  AFB  flutter  tests  were 

conducted  on  the  P-80  to  obtain  data  on  the  effects  of  tip  tank  fuel  c.g. 
travel.  weights  were  used  to  control  the  variables.  The  wing  modes  were 

excited  by  elevator  raps  for  the  symmetrical  case  and  aileron  r^  for  the 
anti-syrtmetrical  case. 

The  test  program  was  planned  to  cover  a  predetermined  ^)eed  range  for 
each  flight.  At  each  data  point  the  pilot  would  excite  the  wings  and  take 
oscillograph  records.  After  each  fli^t  the  data  would  be  analyzed  and  the 
speed  range  established  for  the  next  sortie. 

After  two  uneventful  fUc^ts  the  pilot  felt  that  too  much  time  was  being 
wasted  in  flying  and  collecting  the  data  so  he  decided  to  take  the  initiative 
to  speed  up  the  program.  A  few  more  airbed  increments  were  flown  that  the 
test  card  had  called  for  and  a  wing  bending-torsion  flutter  condition 
developed.  Fortunately,  he  was  able  to  jettison  the  tip  tanks  before  a 
catastrophic’  failure  occurred.  However,  the  wings  were  so  badly  ripped  that 
they  could  not  be  repaired.  Analysis  of  the  subcritical  response  on  the 
oscillograph  showed  that  there  was  sufficient  flutter  onset  warning  even 
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though  the  pilot  could  feel  no  change  in  the  aircraft  prior  to  the  flutter 
occurrence.  The  moral  is  -  develop  a  sound  test  plan  and  stick  to  it. 


12.5.16  Flight  Test  Execution 

There  are  also- do's  and  don'ts  in  the  flying  and  test  plan  execution. 


(1)  Precision  flying  is  a  must  -  do  not  exceed  the  maximum  intended 
aerodynamic  pressure.  Practice  build-up  runs  should  be  flown  so  that  the 
procedxires  becotie  rote  and  the  pilot  knows  exactly  how  the  critical 
speeds  will  be  approached. 

(2)  There  must  be  a  minimum  craii7.  For  multi-place  aircraft  each  crewmember 
must  be  thoroughly  familiar  with  the  egress  procedures  following  an 
emergency.  In  past  cases  of  flutter  testing  of  cargo  type  aircraft  where 
the  risks  v^e  high,  knotted  ropes  were  secured  along  the  floor  in 
converging  lines  toward  the  exits.  If  the  aircraft  became  ixncontrollable 
in  flight  the  crewmembers  would  be  able  to  pull  themselves  toward  the 
exits  for  bailout. 

(3)  The  pilot  must  be  conpletely  familiar  with  the  recovery  techniques  if  a 
flutter  condition  were  to  be  encountered.  There  would  be  an  irrmediate 
need  to  reduce  the  q  and  the  pilot  must  fully  understand  the  effects  on 
the  airplane  of  throttle  chops,  iimiediate  g  loadings,  use  of  S/B,  stores 
release,  and  so  forth.  In  other  words,  \4iat  is  the  exact  response  for 
the  pilot  if  a  flutter  condition  were  encountered. 

(4)  An  investigation  should  start  at  the  hi^er  altitudes,  working  first 
toward  the  Mach  limits.  In  subsequent  flights,  the  lower  altitude,  high  q 
limits  could  be  approached. 

(5)  If  at  any  time  the  danping  becomes  less  than  predicted,  or  if  they  reach 
minimum  planned  for  levels,  recover  the  aircraft  and  regroup. 

(6)  The  weather  on  test  day  should  be  near  perfect.  Turbulence  (unless  used 
as  the  excitation  source)  might  prematurely  excite  the  critical 
condition.  Obstructions  to  vision  such  as  clo;:id  layers  or  heavy  haze  are 
not  desirable  during  high  speed  dives.  Strong  wind  shears  can  lead  to 
dangerous  dive  attitude  changes. 

(7)  A  decision  must  be  made  on  vhether  to  test  over  land  or  over  water. 
There  is  usually  less  turbulence  and  more  altitude  for  dive  recoveries 
over  the  latter.  The  disadvantages  are  distance  to  the  recovery  base  and 
the  difficully  of  ^inreckage  retrieval  should  an  accident  occur. 

(8)  Real  time  IM  with  coiputer  hook  up  is  highly  recommended.  This  set-up 
would  allcw  the  engineers  to  monitor  energy  distribution  of  the  excited 
modes  and  determine  the  adequacy  of  excitation  at  a  subcritical  speed. 
Also,^  the  damping  trends  will  be  better  defined  so  that  engineering 
decisions  based  sinply  on  successi^^e  oscillograph  traces  can  be  avoided. 
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(9)  The  aircrew  shoiald  monitor  the  data  analysis  to  ensure  that  they  clearly 
understand  hew  the  tests  are  proceeding.  Their  inputs  relative  to 

mission  progression  and  foreseeable  problems  are  often  times  invaluable. 
Many  of  the  "rules”  just  presented  are  ccramon  sense  in  natixre. 

Nevertheless,  they -are  enphasized  sinply  because  of  the  many  occasions  in 
previous  flight  tests  v^iere  catastrophes  happen  sirtply  because  of  a  failure 
in  understanding  or  of  comnunication  by  the  aircrew.  Back  in  1941,  flutter 
tests  were  being  conducted  on  the  twin  engined  AT-8.  Ground  vibration  tests 
in  conjunction  with  a  wing  analysis  predicted  a  bending-torsion  flutter  at 
217K.  During  the  preliminary  flight  tests,  however,  no  large  wing 
oscillations  were  encountered  even  at  speeds  up  to  221  K  (only  small 
persistent  oscillations  v^e  recorded) . 

Consequently,  a  3/8  in-lb  rotational  unbalance  shaker  was  placed  at  the 
42%  soni-span  on  the  left  wing  rear  spar.  The  procedure  for  the  pilot  and 
flutter  engineer  era?  was  to  climb  to  15  thousand  feet,  start:  a  gentle  pewered 
dive  and  tune  the  shaker  to  resonance  at  each  incremental  speed  build-up. 
ipproaching  200  K,  the  vying  oscillations  became  so  large  that  the  pilot  became 
sonevbat  concerned.  However,  since  he  felt  that  the  engineer  in  back  was 
watching  the  same  oscillations  without  getting  excited  the  situation  was  imder 
control.  It  turned  out  that  the  flutter  engineer  was  on  his  knees  controlling 
the  knobs  on  the  excitation  and  recording  equipment  and  vyas  not  in  fact 
observing  the  large  oscillations. 

This  exaiiple  illustrates  the  importance  of  each  person  involved  knowing 
vdiat  his  and  everyone  else's  responsibility  is  in  a  hazardous  test  program. 
Another  misconception  was  noted  then  and  at  other  times  prior  and  subsequent; 
If  the  flutter  engineer  is  aboard  things  must  be  okay. 

A  successful  test  program  requires  that  the  crew  thoroughly  understand 
the  problan,  the  aircraft  capabilities  and  characteristics,  the  whys  and 
viberefores  in  the  data  accumulation  and  analyses,  and  the  individual 
responsibilities  in  the  flight  test  execution.  Like  in  all  group  efforts  (and 
so  vrell  exeirplif ied  in  the  game  of  football) ,  a  sound  game  plan ,  followed  by 
total  concentration  in  individual  execution  is  a  must.  Happy  landing. 
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12.5.17  Brief  Exaitple 

We  \d.ll  very  briefly  discuss  the  Northrop  T-38  as  an  exanple  of  how 
flutter  characteristics  of  a  prototype  aircraft  are  detennined  through  design, 
analysis,  tunnel  and  ground  vibration  tests.  Hie  enphasis  is  not  on  detail 
but  on  overall  approach.  Seme  of  the  measures  taken  to  identiify  and  eliminate 
possible  problem  areas  vdll  be  discussed. 

NOn-steatfy  aera^amic  theories  are  usually  conservative  in  that  the 
Ccilculated  speeds  at  vdiich  instability  will  occur  are  usually  lower  than 
actual.  Because  of  this,  and  due  to  the  limitations  of  beam  theory  analyses 
the  basic  flutter  inputs  into  the  T-38  design  philosophy  were  to  be  derived 
frem  wind  tunnel  model  tests.  Analytic  calculations  were  made  mainly  to 
define  speed  incranents  resulting  from  structural  changes  and  variations  in 
the  surface  boundary  conditions. 

Ground  vibrations  tests  on  the  airplane  were  run 'primarily  to  verify  the 
calculated  modes  and  to  illuminate  deficiencies  in  the  structural  boundaries 
and  the  safety  margins  on  the  final  wing  design.  Flutter  flight  test 
confirmed  the  safety  margins  on  the  production  aircraft. 

Hie  vibration  and  flutter  analyses  of  each  of  the  T-38  aerofynaraic 
surfaces  were  based  on  the  following  assunptions:  (1)  Spanwise  distributions 
of  bending  and  torsional  stiffness  were  represented  along  a  straight  swept 
elastic  axis;  (2)  Surface  mass  distributions  were  reparesented  by  inertially 
equivalent  strips  perpendicular  to  the  elastic  axis  as  shown  in  Figure  12.93; 
(3)  Oscillating  surface  aerodynamic  forces  were  approocimated  by 
two-dimensional  aerodynamic  forces  calculated  for  streanwise  strips  as  shotm 
in  Figure  12.94,  and  (4)  A  flutter  mode  vras  r^resented  by  a  si^ierposition  of 
a  finite  number  of  natural  vibration  modes. 
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first  two  assunptions  lead  to  a  static  beam  problan  with  ineirtia 
loadings  induced  by  the  mass  section  oscillations.  The  strip  theory 
aerxx^mamic  terms  used  in  the  flutter  analysis  were  initially  based  on  the 
inocnpressible  flx3w  oscillating  wing  section  theory  of  Bieodorsen.  In 
foUcw-up  analyses  the  lift  curve  slopes  were  adjusted  frcsn  the  infinite  AR  2v 
value  to  wind  tunnel  derived  ones. 

Ihe  three  low  speed  tunnel  models  used  the  spanwise  mass  distribution 
shown  in  Figure  12.93.  Ihe  hi^  qpeed  models  simulated  the  actual  aircraft 
structure  with  scaled  down  skins  and  spars.  All  models  were  vibrated  prior  to 
the  tunnel  tests  to  verify  the  dynamic  siinnla<-ion. 

Bie  general  results  of  the  wing  tests  are  as  follows.  Flutter  boundary 
determination  was  difficult  because  of  the  high  stability  of  the  wing.  No 
flutter  occurred  during  the  transonic  tests  at  CJomell  vAiere  115%  limit 
velocities  vere  achieved  at  Mach  between  0.78  and  0.94.  A  flutter  may 
have  occurred  with  the  addition  of  aft  ballast  at  M  =  0.86.  Based  on  an 
analytical  edxapolation  ccrtparing  the  bare  wing  with  the  ballasted  one  a 
safety  margin  of  32%  was  predicted. 

Tunnel  tests  at  AEDC  revealed  a  possible  low  damping  region  at  M  =  0.85 
and  q  =  2340  (V^  =  850  K) .  The  wing  was  lost  at  q  =  3480.  Even  a 

conservative  estimate  yielded  a  safefy  margin  of  38%. 

Ihe  hi^  speed  wing  tests  at  Cornell  covered  a  range  of  aileron  rotation 
to  second  wing  bending  (t0g/a)2)  from  0.92  to  1.48.  NO  aileron  flutter  or  buzz 
instability  was  observed.  Very  briefly,  the  overall  conclusions  of  the  T-38 
flutter  characteristics  were  as  follows: 


Wrng  and  Aileron  —  no  special  problefiis  were  found  in  this  area. 

Vertical  Fin  and  Rudder  -  except  for  the  destabilizing  effect  of  a  three 
pound  extended  fuel  vent,  the  vertical  fin  did  not  present  any  major 
stability  problems. 

Iforizontal  Stabi  1  izer  -  this  was  initially  the  most  critical  of  the  T-38 
aerodynamic  surfaces.  Flutter  problem  areas  were  associated  with 
actuator  stiffiiess,  torque  tube  bending  stiffness,  freeplay  in  the 
actuator  system,  and  the  ^sibility  of  fuselage  vertical 
bending/horizontal  stabilizer  bendii^  flutter.  TSie  presence  of  freepl^ 
in  the  actuator  system  resulted  in  a  significant  decrease  in  model 
flutter  speeds.  As  a  result,  rigging  procedures  were  adopted  to  elim¬ 
inate  feeeplay  on  the  airplane. 
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The  level  of  pitch  stiffness  on  the  stabilizer  vras  found  to  be  deficient 
during  the  ground  vibration  tests.  Increases  in  horn  size  and  hydraxilic 
cylinder  stiffiiess  \i«ere  required  to  alleviate  the  coupled  stabilizer  pitch  - 
torsion  problem.  In  fact,  one  of  the  more  significant  recaimendations  made  in 
the  report  was  that  more  rigid  controls  ware  required  in  the  future  cn  the 
specific  problem  of  actuator  stiffness  and  fre^lay. 

It  is  appropriate  at  this  time  to  present  a  quote  from  the  "DISCUSSION 
AND  CCXJCUJSIONS"  section  of  the  report. 


"Identification  and  elimination  of  the  1^-38  flutter  prdaLerr  areas  has 
required  a  coordinated  effort  in  the  fields  of  theoretical  analysis, 
ground  vibration  testing,  and  wind  tunnel  flutter  model  testing.  No  one 
of  these  investigations  woiald  have  been  completely  effective  without  the 
additional  information  received  from  the  other  studies.  The  final 
confirmation  of  aircraft  flutter  freedom  has  come  from  the  best  possible 
source  -  fli^t  flutter  tests  of  the  production  aircraft." 
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